Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,685,073 C→T 5.2% P272L (CCG→CTG)  tus → DNA replication terminus site‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,0730CT 5.2% 249.3 / 2.1 77P272L (CCG→CTG) tusDNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base T (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GCCACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTG  >  NC_000913/1684924‑1685222
                                                                                                                                                     |                                                                                                                                                     
gccACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGcc                                                                                                                                                       <  1:929418/150‑1 (MQ=255)
gccACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGcc                                                                                                                                                       >  2:1808366/1‑150 (MQ=255)
gccACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGcc                                                                                                                                                       >  1:1923682/1‑150 (MQ=255)
gccACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGcc                                                                                                                                                       <  2:384967/150‑1 (MQ=255)
        aCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTg                                                                                                                                               <  2:121764/150‑1 (MQ=255)
        aCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTg                                                                                                                                               >  1:1860555/1‑150 (MQ=255)
             aaGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAg                                                                                                                                          <  2:1968668/150‑1 (MQ=255)
                 tAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAgttgt                                                                                                                                      >  1:1153677/1‑150 (MQ=255)
                 tAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAgttgt                                                                                                                                      >  1:954134/1‑150 (MQ=255)
                 tAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAgttgt                                                                                                                                      >  2:1810967/1‑150 (MQ=255)
                  aaaaaTCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAgttgtt                                                                                                                                     <  2:1254938/150‑1 (MQ=255)
                    aaaTCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCTGGACGTTGGTGAGTTGTTaa                                                                                                                                   >  2:1543894/1‑150 (MQ=255)
                    aaaTCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTaa                                                                                                                                   <  2:1967649/150‑1 (MQ=255)
                     aaTCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCTGGACGTTGGTGAGTTGTTaaa                                                                                                                                  >  1:128506/1‑150 (MQ=255)
                     aaTCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCTGGACGTTGGTGAGTTGTTaaa                                                                                                                                  <  2:128506/150‑1 (MQ=255)
                     aaTCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTaaa                                                                                                                                  >  1:1284568/1‑150 (MQ=255)
                         aaaCGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTAc                                                                                                                              >  1:1353319/1‑150 (MQ=255)
                           aCGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGa                                                                                                                            <  1:611907/150‑1 (MQ=255)
                                cGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCg                                                                                                                       <  2:285149/150‑1 (MQ=255)
                                   tGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACa                                                                                                                    >  2:1621626/1‑150 (MQ=255)
                                   tGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACa                                                                                                                    >  2:196847/1‑150 (MQ=255)
                                        gTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTg                                                                                                                                               <  2:644419/118‑1 (MQ=255)
                                        gTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTg                                                                                                                                               >  1:644419/1‑118 (MQ=255)
                                         tGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTgc                                                                                                              <  1:197882/150‑1 (MQ=255)
                                            aGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTgcagc                                                                                                           <  1:1796527/150‑1 (MQ=255)
                                                       cgcgTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATaaa                                                                                                <  1:1810967/150‑1 (MQ=255)
                                                        gcgTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGTGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAAc                                                                                               <  1:1635869/150‑1 (MQ=255)
                                                           tCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTc                                                                                            <  2:155345/150‑1 (MQ=255)
                                                           tCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTc                                                                                            >  2:1639997/1‑150 (MQ=255)
                                                              ggTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAgg                                                                                         <  2:1153677/150‑1 (MQ=255)
                                                              ggTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAgg                                                                                         >  2:1621717/1‑150 (MQ=255)
                                                               gTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGgc                                                                                        <  1:1621717/150‑1 (MQ=255)
                                                                tACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGgcg                                                                                       >  2:1205452/1‑150 (MQ=255)
                                                                tACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGgcg                                                                                       >  1:1807701/1‑150 (MQ=255)
                                                                 aaaaaGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACGGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGAGAATGTGCAGCACCGTTATAAACCTCAGgcgc                                                                                      >  2:1537059/3‑150 (MQ=255)
                                                                 aCAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGgcgc                                                                                      >  1:294496/1‑150 (MQ=255)
                                                                  cAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCa                                                                                     >  2:449358/1‑150 (MQ=255)
                                                                   aaaGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAg                                                                                    <  1:593281/150‑1 (MQ=255)
                                                                     aGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACCCCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGcc                                                                                  >  2:1162561/1‑150 (MQ=255)
                                                                     aGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGcc                                                                                  >  1:950189/1‑150 (MQ=255)
                                                                           tCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCg                                                                            >  1:1888974/1‑150 (MQ=255)
                                                                           tCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCg                                                                            >  2:564840/1‑150 (MQ=255)
                                                                            cAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGt                                                                           >  2:1574521/1‑150 (MQ=255)
                                                                              aaaaCAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCTGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGttt                                                                         <  1:1543894/150‑1 (MQ=255)
                                                                              aaaaCAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGttt                                                                         <  2:1807701/150‑1 (MQ=255)
                                                                              aaaaCAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGttt                                                                         >  1:1171543/1‑150 (MQ=255)
                                                                              aaaaCAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGttt                                                                         <  2:1991215/150‑1 (MQ=255)
                                                                                aaCAAGTCCAACCCGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGCCAATGTGCAGCACCGTTATAAACCTCAGGCGCCGCCGCGTCGTTTGa                                                                       >  2:1713301/1‑150 (MQ=255)
                                                                                 aCAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGAt                                                                      >  1:836942/1‑150 (MQ=255)
                                                                                  cAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATc                                                                     <  1:1574521/150‑1 (MQ=255)
                                                                                  cAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATc                                                                     <  2:1092306/150‑1 (MQ=255)
                                                                                     gTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCAtt                                                                  >  2:1936583/1‑150 (MQ=255)
                                                                                       ccAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTcc                                                                <  1:1808366/150‑1 (MQ=255)
                                                                                            acGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGc                                                           <  1:1205452/150‑1 (MQ=255)
                                                                                             cGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCt                                                          >  2:25601/1‑150 (MQ=255)
                                                                                               ccTGCCCTACACCACCGATTGCCCTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGc                                                        <  2:1605212/150‑1 (MQ=255)
                                                                                                  gCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCAcc                                                     <  2:698422/150‑1 (MQ=255)
                                                                                                     cTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGt                                                  >  2:1671146/1‑150 (MQ=255)
                                                                                                      tACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTa                                                 >  2:951515/1‑150 (MQ=255)
                                                                                                      tACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTa                                                 <  1:1810396/150‑1 (MQ=255)
                                                                                                      tACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTa                                                 <  1:615149/150‑1 (MQ=255)
                                                                                                      tACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTa                                                 <  1:1936583/150‑1 (MQ=255)
                                                                                                            ccTGATTGCAGTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATACCACGGCTGCACCTGTATGTTGc                                           <  2:607748/149‑1 (MQ=255)
                                                                                                                  tGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTa                                     >  2:194496/1‑150 (MQ=255)
                                                                                                                    cACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAAc                                   <  1:949597/150‑1 (MQ=255)
                                                                                                                      cTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGc                                 >  2:1997031/1‑150 (MQ=255)
                                                                                                                        gATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGccc                               >  1:1994214/1‑150 (MQ=255)
                                                                                                                            aaTCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCt                           <  1:1639997/150‑1 (MQ=255)
                                                                                                                             aTCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCtt                          <  2:294496/150‑1 (MQ=255)
                                                                                                                               cGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTc                        >  2:164318/1‑150 (MQ=255)
                                                                                                                               cGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTc                        <  2:954134/150‑1 (MQ=255)
                                                                                                                               cGGGATAATGGCGCGGGCGTGCCGGACCTTGGTGAGTTGTTAAATTACGATGCCCACAATGTGCAGCACCGTTATAAACCTCACGCGCAGCCGCTTCGTTAGATCATTCCACGGCTGCACCACTATGTAGCAGATTAACTCCCGGCTTTc                        >  2:1063741/1‑150 (MQ=255)
                                                                                                                                  gATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATa                     >  2:1958971/1‑150 (MQ=255)
                                                                                                                                   aTAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATAc                    <  1:25601/150‑1 (MQ=255)
                                                                                                                                       tGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGcc                >  1:552276/1‑150 (MQ=255)
                                                                                                                                       tGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGcc                >  2:1778417/1‑150 (MQ=255)
                                                                                                                                        ggCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCg               >  2:170383/1‑150 (MQ=255)
                                                                                                                                             gggCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCAt          >  2:113604/1‑150 (MQ=255)
                                                                                                                                              ggCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATc         <  2:950189/150‑1 (MQ=255)
                                                                                                                                               gCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCt        <  2:1860555/150‑1 (MQ=255)
                                                                                                                                                   gCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTc    <  2:836942/150‑1 (MQ=255)
                                                                                                                                                     cGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTg  >  2:557547/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTG  >  NC_000913/1684924‑1685222

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: