Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,610 C→T 5.8% G16S (GGC→AGC)  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,6100CT 5.8% 326.6 / 4.4 103G16S (GGC→AGC) sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (56/40);  new base T (4/2);  total (60/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 105 total aligned reads displayed.

CCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGC  >  NC_000913/1724462‑1724759
                                                                                                                                                    |                                                                                                                                                     
ccATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCg                                                                                                                                                      >  1:1513138/1‑150 (MQ=255)
  aTGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGt                                                                                                                                                    >  2:486825/1‑150 (MQ=255)
     ttCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCACGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGc                                                                                                                                                 >  1:1451068/1‑150 (MQ=255)
      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:1646322/1‑150 (MQ=255)
                                 ggCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTa                                                                                                                     >  2:1378323/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)
                                      aaaCTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACg                                                                                                                <  2:1653466/150‑1 (MQ=255)
                                          tCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTc                                                                                                            <  1:1622422/150‑1 (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)
                                                        ggTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGAGCCGGTTTCAACAACCAGCTCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCAGtt                                                                                              >  1:995763/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:374294/1‑150 (MQ=255)
                                                                                    gACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAAc                                                                  <  1:861347/150‑1 (MQ=255)
                                                                                       tGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAg                                                               <  1:1378323/150‑1 (MQ=255)
                                                                                       tGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAg                                                               >  2:308459/1‑150 (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)
                                                                                                           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)
                                                                                                                     aTCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAATGGCGGCTATTCCTa                                 <  1:969365/150‑1 (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)
                                                                                                                                              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)
                                                                                                                                                    |                                                                                                                                                     
CCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGC  >  NC_000913/1724462‑1724759

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: