Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,400,393 C→T 6.1% P73L (CCA→CTA)  USA300HOU_1299 → ImpB/MucB/SamB family DNA‑damage repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,400,3930CT 6.1% 185.9 / 2.8 66P73L (CCA→CTA) USA300HOU_1299ImpB/MucB/SamB family DNA‑damage repair protein
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base T (2/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

TTGAAAAGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTATGGAT  >  CP000730/1400255‑1400532
                                                                                                                                          |                                                                                                                                           
ttGAAAAGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACcac                                                                                                                                           >  1:199186/1‑141 (MQ=255)
   aaaaGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACATAAATCACAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATa                                                                                                                                        >  2:617167/1‑141 (MQ=255)
   aaaaGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGAGTCAGGGTTCTGTAATATTGGCTGCGCCACCTAAATTAAAAGAATTAGGCATCAAGTCCGGGTCGCGATTGTTTGAAATACCACATa                                                                                                                                        <  2:435815/141‑1 (MQ=255)
    aaaGGGGCTAGATGCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAg                                                                                                                                       >  2:340682/1‑141 (MQ=255)
    aaaGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAg                                                                                                                                       >  1:3467/1‑141 (MQ=255)
        gggCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAAt                                                                                                                                   >  1:496080/1‑141 (MQ=255)
        gggCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAAt                                                                                                                                   >  1:632195/1‑141 (MQ=255)
        gggCTAAATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAATAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAAt                                                                                                                                   >  2:592349/1‑141 (MQ=255)
         ggCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATg                                                                                                                                  <  2:57118/141‑1 (MQ=255)
           cTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGat                                                                                                                                <  1:945079/141‑1 (MQ=255)
           cTAGATCCATTAGAAACAAAGCCAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGACTTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGat                                                                                                                                <  1:265882/141‑1 (MQ=255)
            tAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACAAAGAAATGata                                                                                                                               >  2:741916/1‑141 (MQ=255)
                tCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTa                                                                                                                           <  2:760803/141‑1 (MQ=255)
                 ccATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTAc                                                                                                                          <  1:9043/141‑1 (MQ=255)
                     tAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGTGATTGTTTGAAATACCACATAGAAATGATATTTACATTa                                                                                                                      <  1:346189/141‑1 (MQ=255)
                       gAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATc                                                                                                                    >  1:700839/1‑141 (MQ=255)
                       gAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGAAACTGTAATATTGGCTGCGACACCTAAATTAAAAGAAGTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATc                                                                                                                    >  2:417530/1‑141 (MQ=255)
                       gAAACAAAGCTAGCTGTTGTTGACAAGACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATc                                                                                                                    >  1:457110/1‑141 (MQ=255)
                            aaaGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATcc                                                                                                               >  2:836642/1‑141 (MQ=255)
                             aaGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGGTTGTTTGAAATACCACATAGAAATGATGTTTACATGATAAATCCa                                                                                                              >  2:15671/1‑141 (MQ=255)
                              aacttgctGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGAGATTGTTTGAAATACAACATAGAAATGATATTAACATTATCAATCCaa                                                                                                             >  2:112480/6‑141 (MQ=255)
                                  agctGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTAt                                                                                                         >  2:127405/1‑141 (MQ=255)
                                      gttgttGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGt                                                                                                     >  1:776514/1‑141 (MQ=255)
                                           tgCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAata                                                                                                <  2:690857/141‑1 (MQ=255)
                                              aGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCt                                                                                             <  2:640468/141‑1 (MQ=255)
                                                  aCTAATCGTCAGGGTTCTGTAATATTGGCTGTGACACCTAGATTAAAAGAATTAGGCATCAAGACAAGGTCGCGATTGTTTGAAGTACCACAGAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAAt                                                                                         >  1:694399/1‑141 (MQ=255)
                                                  aCTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATAACTTAAt                                                                                         >  2:88846/1‑141 (MQ=255)
                                                   cTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGTGTCGCTATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATg                                                                                        <  2:315881/141‑1 (MQ=255)
                                                   cTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATg                                                                                        >  2:922680/1‑141 (MQ=255)
                                                   cTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATg                                                                                        <  2:363800/141‑1 (MQ=255)
                                                      aGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGttt                                                                                     >  1:69587/1‑141 (MQ=255)
                                                       gCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTc                                                                                    >  2:869572/1‑141 (MQ=255)
                                                           cAGGGTTCTGTAATATTGGCTGCAACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGtt                                                                                >  2:674811/1‑141 (MQ=255)
                                                              ggTTCTGTACAATTGGCTGCGACAACTAAATTAAACGCTGGATACATACAGTCAGTGCCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCt                                                                             >  1:690286/1‑141 (MQ=255)
                                                                 tCTGTAATATTGGCTGCGAACGCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTAtt                                                                          >  1:226838/1‑141 (MQ=255)
                                                                       agatggacTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAg                                                                    <  1:440556/134‑1 (MQ=255)
                                                                         atgtgCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTATATTATCAACCCTAGTATGCGTAATTATCTTAATGTTTCAGTTGCTATTTCTAAGat                                                                  <  2:444141/137‑1 (MQ=255)
                                                                             gCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTACCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgca                                                              >  1:297519/1‑141 (MQ=255)
                                                                              cTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcat                                                             >  1:533477/1‑141 (MQ=255)
                                                                              cTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcat                                                             >  2:861732/1‑141 (MQ=255)
                                                                               tgggACACCGAAATTAAAAGAATTAGGGATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcatt                                                            <  1:382235/138‑1 (MQ=255)
                                                                               tGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcatt                                                            <  1:245004/141‑1 (MQ=255)
                                                                                 cGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcattgc                                                          <  2:457110/141‑1 (MQ=255)
                                                                                 cGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcattgc                                                          <  2:700839/141‑1 (MQ=255)
                                                                                   acacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCCTTGCGt                                                        >  2:461839/1‑141 (MQ=255)
                                                                                   acacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGt                                                        >  2:644428/1‑141 (MQ=255)
                                                                                   acacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGt                                                        >  2:285826/1‑141 (MQ=255)
                                                                                    cacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGtt                                                       <  1:741916/141‑1 (MQ=255)
                                                                                         aaaTTAAAAGAATTAGGCATCAAGACAGGTTCGCGATTGTTTGAAATACCACGTAGAAAGGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATAtt                                                  <  2:382858/141‑1 (MQ=255)
                                                                                         aaaTTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTACATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATAtt                                                  <  1:414003/141‑1 (MQ=255)
                                                                                         aaaTTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATAtt                                                  <  1:236526/141‑1 (MQ=255)
                                                                                           aTTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTcc                                                <  1:127405/141‑1 (MQ=255)
                                                                                                 aGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTAAGTTGCTATTTCTAAGATTGCATTGCGTTATTTTCAACCTGa                                          >  2:367562/1‑141 (MQ=255)
                                                                                                  gAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGaa                                         >  1:743421/1‑141 (MQ=255)
                                                                                                   aaTTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTACACCTGAAg                                        >  1:677805/1‑141 (MQ=255)
                                                                                                       acacATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAATTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGAttt                                    <  1:1920/138‑1 (MQ=255)
                                                                                                       aGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGAttt                                    >  1:39576/1‑141 (MQ=255)
                                                                                                           aTCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTacac                                >  2:789464/1‑141 (MQ=255)
                                                                                                            tCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACAcc                               <  2:406068/141‑1 (MQ=255)
                                                                                                                 aCAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAatat                          >  1:153098/1‑141 (MQ=255)
                                                                                                                  cAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGCAGATTTACACCAAtata                         >  2:630257/1‑141 (MQ=255)
                                                                                                                          cgATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGAc                 <  1:303035/141‑1 (MQ=255)
                                                                                                                           gATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACg                <  1:340682/141‑1 (MQ=255)
                                                                                                                           gATTGTTTAAAATACCACATAGAAATGATATTTACATGATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACg                <  1:153428/141‑1 (MQ=255)
                                                                                                                                tttGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAAttt           >  1:108688/1‑141 (MQ=255)
                                                                                                                                   gAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTGATATTCCACCTGAAGATTTACACCAATATGGTATTGAGGAAtttttt        >  1:21422/1‑141 (MQ=255)
                                                                                                                                    aaaTACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTGCTAAGATTGCAGTGGGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTa       <  2:690286/141‑1 (MQ=255)
                                                                                                                                         ccACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTATGGAt  >  2:285383/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TTGAAAAGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTATGGAT  >  CP000730/1400255‑1400532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑