Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,596,244 A→G pseudogene (1744/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,596,2440AG99.1% 405.9 / NA 108pseudogene (1744/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (47/60);  total (48/60)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 108 total aligned reads displayed.

CGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTC  >  NC_000913/1596101‑1596389
                                                                                                                                               |                                                                                                                                                 
cgtcgCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGcgtaccc                                                                                                                                            <  1:23388/151‑1 (MQ=255)
    gCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGcgtacccgtac                                                                                                                                        <  2:357700/151‑1 (MQ=255)
    gCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGcgtacccgta                                                                                                                                         <  2:938192/150‑1 (MQ=255)
       ttGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGcgtacccgtacc                                                                                                                                       >  2:434579/1‑149 (MQ=255)
        tGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGcgtacccgtacc                                                                                                                                       >  3:505693/1‑148 (MQ=255)
        tGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGcgtacccgtacc                                                                                                                                       >  2:905081/1‑148 (MQ=255)
         ggTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGAt                                                                                                                                   <  3:814579/151‑1 (MQ=255)
            gATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTc                                                                                                                                <  2:301794/151‑1 (MQ=255)
              tAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTccc                                                                                                                              >  2:1325316/1‑151 (MQ=255)
               aaaTCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCt                                                                                                                             >  1:147235/1‑151 (MQ=255)
                aaTCCGGCATCCGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCt                                                                                                                             >  3:443367/1‑150 (MQ=255)
                 aTCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTAc                                                                                                                           <  3:1092371/151‑1 (MQ=255)
                   ccGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACg                                                                                                                          <  4:239819/150‑1 (MQ=255)
                   ccGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGt                                                                                                                         <  2:6389/151‑1 (MQ=255)
                        aTCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACat                                                                                                                     >  2:1218205/1‑150 (MQ=255)
                          cTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATAtt                                                                                                                  >  2:82498/1‑151 (MQ=255)
                           tGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTg                                                                                                                 >  3:378740/1‑151 (MQ=255)
                            gccgccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATAtt                                                                                                                  <  4:52338/149‑1 (MQ=255)
                            gccgccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATAtt                                                                                                                  <  4:342132/149‑1 (MQ=255)
                            gccgccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATAtt                                                                                                                  <  4:371136/149‑1 (MQ=255)
                              cgccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCCCCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAAt                                                                                                              <  3:499408/151‑1 (MQ=255)
                              cgccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAAt                                                                                                              <  2:382598/151‑1 (MQ=255)
                                ccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATg                                                                                                             <  3:33001/150‑1 (MQ=255)
                                ccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGt                                                                                                            <  3:441986/151‑1 (MQ=255)
                                ccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGt                                                                                                            <  3:1075207/151‑1 (MQ=255)
                                ccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGt                                                                                                            <  2:964875/151‑1 (MQ=255)
                                ccgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGt                                                                                                            <  2:657220/151‑1 (MQ=255)
                                 cgccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGtt                                                                                                           <  2:1258982/151‑1 (MQ=255)
                                   ccTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTc                                                                                                         <  2:49131/151‑1 (MQ=255)
                                    cTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTc                                                                                                         >  3:1143393/1‑150 (MQ=255)
                                     tCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCa                                                                                                        >  2:1380895/1‑150 (MQ=255)
                                     tCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCa                                                                                                        >  2:277470/1‑150 (MQ=255)
                                       aCCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCaaga                                                                                                     <  3:434612/151‑1 (MQ=255)
                                       aCCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCaaga                                                                                                     <  3:68320/151‑1 (MQ=255)
                                        ccGTGAGCAGCGCCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCaaga                                                                                                     >  3:1177115/1‑150 (MQ=255)
                                        ccGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCaagaa                                                                                                    <  2:549824/151‑1 (MQ=255)
                                        ccGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCaaga                                                                                                     <  3:225385/150‑1 (MQ=255)
                                          gTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAagaaga                                                                                                  <  1:154088/151‑1 (MQ=255)
                                            gAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCaagaag                                                                                                   >  2:925198/1‑148 (MQ=255)
                                            gAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAagag                                                                                                 <  3:1113651/150‑1 (MQ=255)
                                            gAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAagag                                                                                                 <  3:934327/150‑1 (MQ=255)
                                              gcagcaCCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACa                                                                                                                      <  2:453589/127‑1 (MQ=255)
                                              gcagcaCCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACa                                                                                                                      >  3:453623/1‑127 (MQ=255)
                                               cagcaCCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGtt                                                                                             <  1:7927/151‑1 (MQ=255)
                                                 gcaCCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCt                                                                                           >  1:223174/1‑151 (MQ=255)
                                                 gcaCCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCt                                                                                           <  2:1350007/151‑1 (MQ=255)
                                                       aTGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCAt                                                                                     >  2:1296440/1‑151 (MQ=255)
                                                       aTGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCAt                                                                                     >  2:429339/1‑151 (MQ=255)
                                                         gTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATc                                                                                    <  1:476700/150‑1 (MQ=255)
                                                               cGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTc                                                                             <  3:825167/151‑1 (MQ=255)
                                                                gTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCa                                                                            <  3:492695/151‑1 (MQ=255)
                                                                  tCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCAc                                                                           <  3:925106/150‑1 (MQ=255)
                                                                   ccTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACg                                                                          <  1:234013/150‑1 (MQ=255)
                                                                       cagcagGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGcc                                                                     <  2:1255535/151‑1 (MQ=255)
                                                                           agGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCtt                                                                  >  4:139081/1‑150 (MQ=255)
                                                                             gCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACt                                                               >  3:696241/1‑151 (MQ=255)
                                                                             gCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACt                                                               >  2:654044/1‑151 (MQ=255)
                                                                             gCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACt                                                               <  3:942235/151‑1 (MQ=255)
                                                                                cAAAAACACCTAAGTAAACTTCTCCACCTTCGACATCCAGCGTACCGTTATTCGTCAAGGTTAGCTTACCCGTTCCAGATGTCCCTACGTACATATGGAATGTTTAAAGAAGAGAGTTCTGCCCAACCCCCCGCACGTCGCCCTTAc                                                                <  3:747213/147‑1 (MQ=255)
                                                                                 aaaaaCACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTg                                                           <  3:1325386/151‑1 (MQ=255)
                                                                                      cacCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACtttt                                                            <  1:226235/145‑1 (MQ=255)
                                                                                      cacCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACtttt                                                            >  4:226233/1‑145 (MQ=255)
                                                                                      cacCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCAtt                                                      <  4:442326/151‑1 (MQ=255)
                                                                                       acCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTg                                                     >  2:716246/1‑151 (MQ=255)
                                                                                        ccTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTga                                                    <  2:209022/151‑1 (MQ=255)
                                                                                           aaGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTgag                                                   <  2:304574/149‑1 (MQ=255)
                                                                                             gTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCt                                               <  1:114279/151‑1 (MQ=255)
                                                                                              tAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCt                                               >  3:909892/1‑150 (MQ=255)
                                                                                               aaaCTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTa                                              <  1:96910/150‑1 (MQ=255)
                                                                                                 aCTTCTCCCCCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCt                                           <  3:988516/151‑1 (MQ=255)
                                                                                                 aCTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCt                                           <  3:961460/151‑1 (MQ=255)
                                                                                                    tctcCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCtg                                          <  4:260090/149‑1 (MQ=255)
                                                                                                      tcCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTgtgt                                       >  1:226757/1‑150 (MQ=255)
                                                                                                       ccACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTgtgt                                       <  2:341501/149‑1 (MQ=255)
                                                                                                        cACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATc                                    >  2:1390160/1‑151 (MQ=255)
                                                                                                        cACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATc                                    <  3:1296509/151‑1 (MQ=255)
                                                                                                        cACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATc                                    >  2:1106392/1‑151 (MQ=255)
                                                                                                           cTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCAcg                                 >  1:86694/1‑151 (MQ=255)
                                                                                                            ttCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCAcg                                 >  1:226080/1‑150 (MQ=255)
                                                                                                             tCGACATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCAcgc                                <  3:1218266/150‑1 (MQ=255)
                                                                                                                aCATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCtt                            >  3:695963/1‑151 (MQ=255)
                                                                                                                aCATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCtt                            >  2:17718/1‑151 (MQ=255)
                                                                                                                aCATTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTCCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCtt                            >  2:1117213/1‑151 (MQ=255)
                                                                                                                  aTTCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTa                          <  3:925233/151‑1 (MQ=255)
                                                                                                                   ttCAGCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTa                          >  3:10010/1‑150 (MQ=255)
                                                                                                                       gCGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACtat                     >  4:40741/1‑151 (MQ=255)
                                                                                                                        cGTACCGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACtat                     >  3:429077/1‑150 (MQ=255)
                                                                                                                             cGTTATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTcc               >  2:470428/1‑151 (MQ=255)
                                                                                                                                tATTCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCggt            >  1:313001/1‑151 (MQ=255)
                                                                                                                                  ttCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCggtg           <  2:1140180/150‑1 (MQ=255)
                                                                                                                                   tCGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCggtggt         >  3:583975/1‑151 (MQ=255)
                                                                                                                                    cGTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCggtggt         >  4:446768/1‑150 (MQ=255)
                                                                                                                                     gTCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTaa       <  2:378705/151‑1 (MQ=255)
                                                                                                                                      tCAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAat      >  3:721301/1‑151 (MQ=255)
                                                                                                                                       cAGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAat      >  3:1096639/1‑150 (MQ=255)
                                                                                                                                        aGGGTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAatat    >  1:326869/1‑151 (MQ=255)
                                                                                                                                         gggTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATAtt   >  3:1277133/1‑151 (MQ=255)
                                                                                                                                          ggTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTc  >  3:466647/1‑151 (MQ=255)
                                                                                                                                           gTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTc  >  1:323483/1‑150 (MQ=255)
                                                                                                                                           gTTAGCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTc  >  3:334739/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACACCTAAGTAAACTTCTCCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTC  >  NC_000913/1596101‑1596389

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: