Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,617 A→T 5.9% V13V (GTT→GTA) ‡ sodC ← superoxide dismutase, Cu, Zn, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,724,6170AT 5.9% 287.5 / 6.3 102V13V (GTT→GTA) ‡sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (51/44);  new base T (4/2);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

Only 100 of 104 total aligned reads displayed.

TCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGCTTTA  >  NC_000913/1724468‑1724766
                                                                                                                                                     |                                                                                                                                                     
tCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGca                                                                                                                                                       >  2:1609850/1‑150 (MQ=255)
      gggggTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAAcc                                                                                                                                                 <  2:1358924/150‑1 (MQ=255)
          gTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAgcg                                                                                                                                             >  2:943509/1‑150 (MQ=255)
          gTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAgcg                                                                                                                                             <  2:398872/150‑1 (MQ=255)
          gTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAgcg                                                                                                                                             >  2:1533333/1‑150 (MQ=255)
              tGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCTGGTTGTTGCAACCGgcgcaca                                                                                                                                         <  2:748933/150‑4 (MQ=255)
              tGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAg                                                                                                                                         <  2:1487188/150‑1 (MQ=255)
               gCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGa                                                                                                                                        >  1:1048056/1‑150 (MQ=255)
               gCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGa                                                                                                                                        >  1:368820/1‑150 (MQ=255)
                 tttCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCTGGTTGTTGCAACCGgcgcacaagc                                                                                                                                      <  2:358343/150‑7 (MQ=255)
                  ttCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCACTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACAAGCGCCAGAATa                                                                                                                                     >  1:1399480/1‑150 (MQ=255)
                    cataTCGCGCGAACAATCCAGACCTTTATCGCTTTCAGTAATGTTGACGCGACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGc                                                                                                                                   <  1:302507/147‑1 (MQ=255)
                    cAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGc                                                                                                                                   >  1:1319871/1‑150 (MQ=255)
                     aGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGATTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGcc                                                                                                                                  >  1:1065077/1‑150 (MQ=255)
                       aTCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAg                                                                                                                                >  1:1435447/1‑150 (MQ=255)
                         cGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGAc                                                                                                                              <  2:1624910/150‑1 (MQ=255)
                           ggCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTa                                                                                                                            >  2:1378323/1‑150 (MQ=255)
                           ggCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTa                                                                                                                            >  2:1646322/1‑150 (MQ=255)
                            gCGAAAACTCCAGACCTTTATCGGTTTCAGTAATTGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTaa                                                                                                                           <  1:22469/150‑1 (MQ=255)
                            gCGAAAACTCCAGACCTTTATCGGTTGCAGTTATGCTGACGCTACCAATTGACTCCCCTACCCCTTGCGACGTGACGAGGATCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTaa                                                                                                                           <  1:1152427/150‑1 (MQ=255)
                              gAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTaaaa                                                                                                                         >  2:861347/1‑150 (MQ=255)
                               aaaaCTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTTCGCCGGTTGCAACAACCCGCGCCAGAATAGCCAGACTAAAAc                                                                                                                        >  1:302218/1‑150 (MQ=255)
                               aaaaCTCCAGACCTTATTCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTGGCAACAACCAGCGCCAGAATAGCCAGACTAAAAc                                                                                                                        <  2:1085114/150‑1 (MQ=255)
                                aaaCTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACg                                                                                                                       >  1:1618762/1‑150 (MQ=255)
                                    tCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTc                                                                                                                   <  2:32824/150‑1 (MQ=255)
                                     ccAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCa                                                                                                                  >  2:1609766/1‑150 (MQ=255)
                                       aGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATa                                                                                                                <  1:1150007/150‑1 (MQ=255)
                                         aCCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAgg                                                                                                              <  2:546319/150‑1 (MQ=255)
                                          ccTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGa                                                                                                             <  1:1718581/150‑1 (MQ=255)
                                          ccTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCACTTGAGACGTGACGAGGTTCATCTAGAATTTTTCACTTGCAGATTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGa                                                                                                             >  1:699438/1‑150 (MQ=255)
                                             ttATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGAGACGAGAAGAGGTTCATCTCGAATTTTTAACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCt                                                                                                          >  1:1638392/1‑150 (MQ=255)
                                                  ggTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGtt                                                                                                     >  2:1548227/1‑150 (MQ=255)
                                                  ggTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGtt                                                                                                     >  2:1354117/1‑150 (MQ=255)
                                                  ggTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGAGCCGGTTTCAACAACCAGCTCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCAGtt                                                                                                     >  1:995763/1‑150 (MQ=255)
                                                      tCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCAtt                                                                                                 >  2:312244/1‑150 (MQ=255)
                                                       cAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                <  2:1513138/150‑1 (MQ=255)
                                                       cAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                <  1:1533333/150‑1 (MQ=255)
                                                       cAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                <  2:1319871/150‑1 (MQ=255)
                                                       cAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                >  1:1259455/1‑150 (MQ=255)
                                                       cAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                <  1:943509/150‑1 (MQ=255)
                                                         gTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtgtg                                                                                              >  2:1356720/1‑150 (MQ=255)
                                                          tAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGTCAGCTTGTGCGCCGGTTGCAACTACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGc                                                                                             <  1:1064974/150‑1 (MQ=255)
                                                            aTGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCAt                                                                                           <  1:1428659/150‑1 (MQ=255)
                                                             tGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGATTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCAtt                                                                                          <  2:1065077/150‑1 (MQ=255)
                                                                    gCTCCCAATTGACTTCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCGTGTGCGCCGGTTGCAACTACCAGCGCCATAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAAt                                                                                   <  2:925875/150‑1 (MQ=255)
                                                                    gCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAAt                                                                                   <  2:1002717/150‑1 (MQ=255)
                                                                      tACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATaa                                                                                 >  2:191440/1‑150 (MQ=255)
                                                                         cAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAgtg                                                                              <  1:1646322/150‑1 (MQ=255)
                                                                            ttGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTaa                                                                           <  2:1451068/150‑1 (MQ=255)
                                                                              gACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACCACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAAc                                                                         >  2:227919/1‑150 (MQ=255)
                                                                              gACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAAc                                                                         <  1:861347/150‑1 (MQ=255)
                                                                              gACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAAc                                                                         >  1:374294/1‑150 (MQ=255)
                                                                                 tGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAg                                                                      >  2:308459/1‑150 (MQ=255)
                                                                                 tGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAg                                                                      <  1:1378323/150‑1 (MQ=255)
                                                                                   cccTACCCCTTGCGACGTGACGAGGTTCAGCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCTCCAGAATAGCCATACTAAAACGTTTCATAGGACCTCCATTGATTGTGCATTCCTGAATAAGTGTAAACCAGTg                                                                    <  1:960426/150‑1 (MQ=255)
                                                                                      tACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGaca                                                                 <  1:312244/150‑1 (MQ=255)
                                                                                          ccTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACAtt                                                             >  2:1344996/1‑150 (MQ=255)
                                                                                            ttGCGATGTGTCGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACAtttt                                                           >  1:237748/1‑150 (MQ=255)
                                                                                            ttGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACAtttt                                                           <  2:249210/150‑1 (MQ=255)
                                                                                             tGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACAttttt                                                          >  2:1551000/1‑150 (MQ=255)
                                                                                              tcGACGTGACGAGGGTCATCTCGACTGTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTg                                                         <  2:72118/149‑1 (MQ=255)
                                                                                                gACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGaa                                                       >  2:414910/1‑150 (MQ=255)
                                                                                                gACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGaa                                                       >  1:194204/1‑150 (MQ=255)
                                                                                                   gTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAAcgg                                                    >  1:1488852/1‑150 (MQ=255)
                                                                                                     gACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAAcggcg                                                  >  1:29251/1‑150 (MQ=255)
                                                                                                       cGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGGAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTACGTTCATAGTGCATACCTGAATAAGAGTAAACCAGTGACAAATTTTTGAAAggctgc                                                >  1:1113552/1‑150 (MQ=255)
                                                                                                          ggTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTAt                                             >  1:1478155/1‑150 (MQ=255)
                                                                                                            ttCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTc                                           >  1:1223548/1‑150 (MQ=255)
                                                                                                                tctcGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAACCGTTTCATAGGACCTCCGTACATTGTGCATTCCTGAATAAGTGTAACCCAGTGACACATTTTTGAACGGCGGCTATTCCTaa                                       >  1:1027006/1‑150 (MQ=255)
                                                                                                                 ctcGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTaaa                                      >  1:299322/1‑150 (MQ=255)
                                                                                                                 ctcGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCAAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTaaa                                      >  1:264705/1‑150 (MQ=255)
                                                                                                                     aCTTTTTCTCTGGCAGCTTGATCGCAGGTTGCAACAACCAGCGACAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTTACACATTTTTGAACGGCGGCTATTCCTAAAAGTg                                  <  2:88658/150‑1 (MQ=255)
                                                                                                                     aCTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTg                                  <  1:166666/150‑1 (MQ=255)
                                                                                                                     aCTTTTTCACTGGCAGCTCGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTg                                  <  1:534715/150‑1 (MQ=255)
                                                                                                                      cTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGc                                 >  1:1190129/1‑150 (MQ=255)
                                                                                                                      cTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGc                                 >  1:1346001/1‑150 (MQ=255)
                                                                                                                        ttttCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCtt                               <  2:368820/150‑1 (MQ=255)
                                                                                                                          ttCACTGGCAGCTTGTGCGCCGGTTGCAACCACCAGCGCCAGAATAGCCAAACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTAc                             >  2:816489/1‑150 (MQ=255)
                                                                                                                          ttCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTAc                             <  1:1609766/150‑1 (MQ=255)
                                                                                                                            cACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTTAACGGCGGCTATTCCTAAAAGTGCTTACgg                           <  2:1399480/150‑1 (MQ=255)
                                                                                                                            cACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTTAAAGTGCTTACgg                           >  2:62783/1‑150 (MQ=255)
                                                                                                                              cTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTa                         >  2:670185/1‑150 (MQ=255)
                                                                                                                              cTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTa                         >  1:1198463/1‑150 (MQ=255)
                                                                                                                                ggCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTAcg                       <  2:1618762/150‑1 (MQ=255)
                                                                                                                                ggCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTAcg                       >  1:1588979/1‑150 (MQ=255)
                                                                                                                                    gCTTGTGCGCCGGTTGCCACAACCAGCGCCAGAATAGCCAGACTAAACCGTTTCATAGGACCTGCGTTCATTGTGCTTTCCTGCATAAGTGTAAACCAGTGCCCCCGTTTTGAACGGCGGCTATTCCTAAAAGTTCTTACGGTaagaagt                   >  2:814179/1‑144 (MQ=255)
                                                                                                                                       tgtgCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAct                <  1:308459/150‑2 (MQ=255)
                                                                                                                                        tggcgcTGGTTGTTGCAACCGGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTTCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAccc               >  1:980834/3‑150 (MQ=255)
                                                                                                                                        gtgCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAccc               >  1:1361822/1‑150 (MQ=255)
                                                                                                                                        gtgCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAccc               <  2:932714/150‑1 (MQ=255)
                                                                                                                                         tgCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAcccc              >  1:46679/1‑150 (MQ=255)
                                                                                                                                         ggcgcTGGTTGTTGCAACCGGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAcccc              >  2:1016319/2‑150 (MQ=255)
                                                                                                                                         ggcgcTGGTTGTTGCAACCGGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAcccc              >  1:752621/2‑150 (MQ=255)
                                                                                                                                          gcgcTGGTTGTTGCAACCGGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCa             >  2:498795/1‑150 (MQ=255)
                                                                                                                                          gcgaCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTACGTTCATTGTGCATTCCTGAATAAATGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACACCa             >  1:178648/1‑150 (MQ=255)
                                                                                                                                           cgcCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAg            >  2:49944/1‑150 (MQ=255)
                                                                                                                                             ccGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTg          <  1:1609850/150‑1 (MQ=255)
                                                                                                                                              cGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGc         <  1:473202/150‑1 (MQ=255)
                                                                                                                                                gTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCg       <  2:1259455/150‑1 (MQ=255)
                                                                                                                                                     aacaacCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGCTTTa  >  2:360677/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGCTTTA  >  NC_000913/1724468‑1724766

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