Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,624 A→G 6.7% L11P (CTG→CCG)  sodC ← superoxide dismutase (Cu‑Zn)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,724,6240AG 6.7% 270.0 / 7.1 91L11P (CTG→CCG) sodCsuperoxide dismutase (Cu‑Zn)
Reads supporting (aligned to +/- strand):  ref base A (41/43);  new base G (3/3);  total (44/47)
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.25e-01

GGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGCTTTACG  >  NC_000913/1724476‑1724768
                                                                                                                                                    |                                                                                                                                                
gggTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGCCGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAg                                                                                                                                                 <  2:1426808/150‑1 (MQ=255)
  gTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCTGGTTGTTGCAACCGgcg                                                                                                                                               >  2:1718121/1‑150 (MQ=255)
   tAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCTGGTTGTTGCAACCGgcgc                                                                                                                                              <  2:716996/150‑1 (MQ=255)
     aTGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCa                                                                                                                                            >  1:1374616/1‑150 (MQ=255)
      tGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAg                                                                                                                                           >  1:659685/1‑150 (MQ=255)
      tGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAg                                                                                                                                           >  2:787572/1‑150 (MQ=255)
        cTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCTGGTTGTTGCAACCGgcgcacaag                                                                                                                                         <  1:1718121/150‑6 (MQ=255)
        cTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGaa                                                                                                                                         <  2:1513816/150‑1 (MQ=255)
         tttCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAAt                                                                                                                                        <  2:1586799/150‑1 (MQ=255)
           tCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAg                                                                                                                                      >  2:236358/1‑150 (MQ=255)
              gATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCa                                                                                                                                   >  2:1282955/1‑150 (MQ=255)
              gATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCa                                                                                                                                   <  1:1624904/150‑1 (MQ=255)
                    gCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTaa                                                                                                                             <  2:56893/150‑1 (MQ=255)
                      gAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTaaaa                                                                                                                           <  2:1715913/150‑1 (MQ=255)
                      gAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTaaaa                                                                                                                           <  1:1389421/150‑1 (MQ=255)
                       aaaaCTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAAc                                                                                                                          <  1:593388/150‑1 (MQ=255)
                        aaaCTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACg                                                                                                                         >  1:701917/1‑150 (MQ=255)
                              cAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCAGACTAAAACGTTTCATa                                                                                                                  <  2:1968705/150‑1 (MQ=255)
                                gACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAg                                                                                                                 >  2:17157/1‑150 (MQ=255)
                                    tttATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGAcc                                                                                                             <  1:787572/150‑1 (MQ=255)
                                      tATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTc                                                                                                           <  2:1968106/150‑1 (MQ=255)
                                          ggTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGtt                                                                                                       >  1:434255/1‑150 (MQ=255)
                                           gTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTc                                                                                                      <  1:1885032/150‑1 (MQ=255)
                                           gTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCACAGGACCTCCGtta                                                                                                      >  2:638087/1‑149 (MQ=255)
                                               cAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGAc                                                                                                              <  1:1818121/138‑1 (MQ=255)
                                               cAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGAc                                                                                                              >  2:1818121/1‑138 (MQ=255)
                                               cAGTAATGGTGACGCGACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                  <  2:1052448/150‑1 (MQ=255)
                                               cAGAAATGTTGACGCTACAAATTGACAGCCATACCCCTTGCCACGTGACGAGGTTAATCTACAATTATTCACAGGCAGCTTGTGCGACCGTAGCAACAACAAGCGCCACAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCAttt                                                                                                  >  2:338097/1‑149 (MQ=255)
                                                aGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtgt                                                                                                 >  1:170715/1‑150 (MQ=255)
                                                aGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                  <  2:686681/149‑1 (MQ=255)
                                                aGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtg                                                                                                  >  1:686681/1‑149 (MQ=255)
                                                 gTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATtgtg                                                                                                <  2:659685/150‑1 (MQ=255)
                                                        tGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCt                                                                                         >  1:1152540/1‑150 (MQ=255)
                                                         gACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTTGTAGAACCTGAGCCAGAATAGCCAGAGTTATACTTTTCATTCGACCTCTGTTCATTGTGAATTCCTg                                                                                        <  2:1859369/150‑1 (MQ=255)
                                                         gACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTg                                                                                        >  2:83954/1‑150 (MQ=255)
                                                         gACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTg                                                                                        <  1:236358/150‑1 (MQ=255)
                                                         gACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAAACAGATACAGAATAGCCAGACTAAAACTTTTCATATGACATCCGTTCATTGTGCATTACTg                                                                                        <  2:737381/150‑1 (MQ=255)
                                                           cGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCTGGTTGTTGCAACCGgcgc                                                                                                                                              >  1:1382793/1‑94 (MQ=255)
                                                           cGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCTGGTTGTTGCAACCGgcgc                                                                                                                                              <  2:1382793/94‑1 (MQ=255)
                                                            gCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAAt                                                                                     <  1:844979/150‑1 (MQ=255)
                                                               aCCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACACCCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAg                                                                                  >  1:1306476/1‑150 (MQ=255)
                                                                  aaTTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGtt                                                                                                                       >  1:533143/1‑110 (MQ=255)
                                                                  aaTTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGtt                                                                                                                       <  2:533143/110‑1 (MQ=255)
                                                                      gACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAAc                                                                           >  1:1863745/1‑150 (MQ=255)
                                                                      gACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAAc                                                                           >  2:1701350/1‑150 (MQ=255)
                                                                         tGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAg                                                                        <  1:151831/150‑1 (MQ=255)
                                                                          gCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGt                                                                       >  2:174876/1‑150 (MQ=255)
                                                                          gCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGt                                                                       >  1:738473/1‑150 (MQ=255)
                                                                            ccTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGa                                                                     <  2:1636995/150‑1 (MQ=255)
                                                                                ccacTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGAcaca                                                                 >  2:508083/4‑150 (MQ=255)
                                                                                   cTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACAttt                                                              <  2:170715/150‑1 (MQ=255)
                                                                                    ttGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACAtttt                                                             <  2:1863745/150‑1 (MQ=255)
                                                                                       cGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGa                                                          >  2:737032/1‑150 (MQ=255)
                                                                                       cGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGa                                                          >  2:1546215/1‑150 (MQ=255)
                                                                                        gACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGaa                                                         >  2:1534557/1‑150 (MQ=255)
                                                                                            tGACGATGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCTCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTTCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAAcggc                                                     <  2:1602935/150‑1 (MQ=255)
                                                                                             gACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAAcggcg                                                    <  1:1368956/150‑1 (MQ=255)
                                                                                              aCGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAAcggcgg                                                   >  2:233233/1‑150 (MQ=255)
                                                                                               cGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACggcggc                                                  >  2:90619/1‑150 (MQ=255)
                                                                                               cGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACggcggc                                                  <  1:1701350/150‑1 (MQ=255)
                                                                                                      cATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCt                                           >  2:731523/1‑150 (MQ=255)
                                                                                                        tctcGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCAAGACTAAAACGTTTCATAGGACCTCCGTTCCTTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTaa                                         >  1:813859/1‑150 (MQ=255)
                                                                                                           cGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAg                                      <  2:434255/150‑1 (MQ=255)
                                                                                                           cGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAg                                      <  1:508083/150‑1 (MQ=255)
                                                                                                             aCTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTg                                    >  2:1152984/1‑150 (MQ=255)
                                                                                                             aCTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTg                                    <  1:90619/150‑1 (MQ=255)
                                                                                                              cTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGc                                   >  1:363559/1‑150 (MQ=255)
                                                                                                               tttttCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCt                                  <  1:1282955/150‑1 (MQ=255)
                                                                                                                  ttCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTAc                               <  2:250686/150‑1 (MQ=255)
                                                                                                                  ttCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTAc                               <  1:233233/150‑1 (MQ=255)
                                                                                                                    cACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACgg                             >  2:204891/1‑150 (MQ=255)
                                                                                                                     aCTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGt                            <  1:260351/150‑1 (MQ=255)
                                                                                                                      cTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCt                                           <  2:1885767/134‑1 (MQ=255)
                                                                                                                      cTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCt                                           >  1:1885767/1‑134 (MQ=255)
                                                                                                                       tGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTAc                          >  1:574367/1‑150 (MQ=255)
                                                                                                                        ggCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCt                                  <  2:543146/141‑1 (MQ=255)
                                                                                                                        ggCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCt                                  >  1:543146/1‑141 (MQ=255)
                                                                                                                        ggCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTAcg                         >  1:1813421/1‑150 (MQ=255)
                                                                                                                         gCAGCTTGTGCGCCGGTTGCAACAACCAGCGGCATAATATCGAGACTAATACGTTTAATAGTACGTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACTTTTTTGATCGGCGGCTATTCCTAAAAGTGCTTAGGGTAcgt                        <  1:338097/150‑1 (MQ=255)
                                                                                                                            gCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTAcgtcgt                     >  2:1409501/1‑150 (MQ=255)
                                                                                                                               tgtgCGCCGGTTGCAACAACCAGCGCCAGAATAGCAAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAcc                  >  1:705567/1‑150 (MQ=255)
                                                                                                                                 ggcgcTGGTTGTTGCAACCGGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAcccc                >  2:959538/2‑150 (MQ=255)
                                                                                                                                 ggcgcTGGTTGTTGCAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTAcccc                >  1:914001/2‑150 (MQ=255)
                                                                                                                                    gcCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGt             >  2:1298144/1‑150 (MQ=255)
                                                                                                                                     ccGGTTGCAACTACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATATGTGTTAAGCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTTCGGTACGTCGTAGCGCAGTg            <  2:429943/150‑1 (MQ=255)
                                                                                                                                     ccGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTg            <  1:204891/150‑1 (MQ=255)
                                                                                                                                        gttgttgcaacCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCg         <  2:914001/143‑1 (MQ=255)
                                                                                                                                        gTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCg         >  2:1193345/1‑150 (MQ=255)
                                                                                                                                         ttGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGc        <  1:17157/150‑1 (MQ=255)
                                                                                                                                               caacCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGCTTTACg  <  2:738473/150‑1 (MQ=255)
                                                                                                                                               caacCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGCTTTACg  <  1:174876/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGTGCGCCGGTTGCAACAACCAGCGCCAGAATAGCCAGACTAAAACGTTTCATAGGACCTCCGTTCATTGTGCATTCCTGAATAAGTGTAAACCAGTGACACATTTTTGAACGGCGGCTATTCCTAAAAGTGCTTACGGTACGTCGTACCCCAGTGCCGCTTTACG  >  NC_000913/1724476‑1724768

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: