Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,417 T→C 6.5% intergenic (+335/‑37) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,4170TC 6.5% 355.3 / 11.3 124intergenic (+335/‑37)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base T (59/56);  new base C (6/2);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 127 total aligned reads displayed.

AGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGC  >  NC_000913/1426272‑1426566
                                                                                                                                                 |                                                                                                                                                     
aGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGGGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgc                                                                                                                                                   <  2:1428920/150‑1 (MQ=255)
     cAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTgg                                                                                                                                              <  1:1879846/150‑1 (MQ=255)
     cAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTgg                                                                                                                                              <  1:1066932/150‑1 (MQ=255)
     cAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAAATAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTgg                                                                                                                                              <  2:156419/150‑1 (MQ=255)
      aGTGTTTGTAGTTGTTAATTACGTACACTCTGTTACAGTGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTggg                                                                                                                                             <  2:1487034/150‑1 (MQ=255)
          tgtgtAGTTGTTAATTCCGTACACTCTGTTACAGAGATGTGATGGGGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttt                                                                                                                                         <  1:1716396/150‑1 (MQ=255)
          tgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttt                                                                                                                                         <  1:1724566/150‑1 (MQ=255)
           gtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGttttt                                                                                                                                        <  2:414998/150‑1 (MQ=255)
           gtgtAGTTGTTAATTACGTACAATCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCTATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGttttt                                                                                                                                        <  2:1486006/150‑1 (MQ=255)
            tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGATAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCAACTTCGCGGTGGGtttttt                                                                                                                                       >  2:476742/1‑150 (MQ=255)
            tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCTCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttttt                                                                                                                                       >  1:340179/1‑150 (MQ=255)
            tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttttt                                                                                                                                       <  1:898904/150‑1 (MQ=255)
             gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttg                                                                                                                                      <  2:507871/150‑1 (MQ=255)
             gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttg                                                                                                                                      <  2:1825038/150‑1 (MQ=255)
             gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                      <  1:982485/150‑5 (MQ=255)
             gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                      <  1:315059/150‑5 (MQ=255)
                gttgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgttg                                                                                                                                   >  1:965654/1‑150 (MQ=255)
                  tgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCTCTGTGGGTTTTTTGTtctc                                                                                                                                 >  2:334774/1‑150 (MQ=255)
                  tgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTc                                                                                                                                 >  1:428997/1‑150 (MQ=255)
                  tgttAAATACGTACACTCTGTTACAGAGATGTGATGGGGTCTTAAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTc                                                                                                                                 <  2:1649911/150‑1 (MQ=255)
                    ttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAg                                                                                                                               <  2:673239/150‑1 (MQ=255)
                        ttACGTACACTCTGTTACAGAGATGTGATTGTGTCTTTACTTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTAATATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGt                                                                                                                           <  2:899999/150‑1 (MQ=255)
                         tACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGtt                                                                                                                          >  1:1973077/1‑150 (MQ=255)
                         tACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGtt                                                                                                                          <  2:1397046/150‑1 (MQ=255)
                           cGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGtttt                                                                                                                        <  1:1144968/150‑1 (MQ=255)
                            gTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTta                                                                                                                       >  2:106597/1‑150 (MQ=255)
                             tACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaa                                                                                                                      <  1:799749/150‑1 (MQ=255)
                             tACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaa                                                                                                                      <  2:162583/150‑1 (MQ=255)
                                   ctGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATg                                                                                                                <  2:211844/150‑1 (MQ=255)
                                   ctGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATg                                                                                                                >  2:539564/1‑150 (MQ=255)
                                         cAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCagag                                                                                                          >  2:1520208/1‑150 (MQ=255)
                                           gagaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGTTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGATCa                                                                                                        >  2:1005770/1‑150 (MQ=255)
                                            agaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCaa                                                                                                       >  1:409471/1‑150 (MQ=255)
                                              aTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAc                                                                                                     <  2:544342/150‑1 (MQ=255)
                                               tgtgATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCATGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                                    <  1:987307/150‑1 (MQ=255)
                                               tgtgATGGTGTCTTTAAATAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                                    <  1:406107/150‑1 (MQ=255)
                                                 tgATGTTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACACACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTc                                                                                                  <  2:389119/150‑1 (MQ=255)
                                                         tCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGc                                                                                                         >  2:1537285/1‑135 (MQ=255)
                                                         tCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGc                                                                                                         <  1:1537285/135‑1 (MQ=255)
                                                         tCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCg                                                                                          <  2:332109/150‑1 (MQ=255)
                                                          cTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACATCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGc                                                                                         <  1:294729/150‑1 (MQ=255)
                                                           tttAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCa                                                                                        >  1:806491/1‑150 (MQ=255)
                                                           tttAATTAAATCAATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCa                                                                                        >  1:523859/1‑150 (MQ=255)
                                                              aattaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaata                                                                                     <  1:427967/150‑1 (MQ=255)
                                                                 taaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCTCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAatt                                                                                  <  1:1768156/150‑1 (MQ=255)
                                                                 taaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAatt                                                                                  >  2:1726516/1‑150 (MQ=255)
                                                                 taaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAatt                                                                                  >  1:1066825/1‑150 (MQ=255)
                                                                  aattCGATGATGCTCCTGGAGAAAATCATTGCTTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAatta                                                                                 <  1:1569929/147‑1 (MQ=255)
                                                                   aaTCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAattat                                                                                >  2:971558/1‑150 (MQ=255)
                                                                   aaTCGATGATGCACCTGGAGAAAAGCATTGCTTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCAACTTCGCGGTGGGTTTTTTGTTGACAGGAGTTATAATAAATCTCAGAGCAAACCTCGCGTCTCGAAATAattat                                                                                >  1:1233439/1‑150 (MQ=255)
                                                                     tCGATGAAGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAATTGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTg                                                                              <  1:245016/150‑1 (MQ=255)
                                                                       gatgatGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCTCTATATCTACTATTATGTCGCCTGAAACCCACTTCTCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGAt                                                                            <  2:2011788/150‑1 (MQ=255)
                                                                             gCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACt                                                                      >  2:1689815/1‑150 (MQ=255)
                                                                              cTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTg                                                                     >  2:1500376/1‑150 (MQ=255)
                                                                                 cTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGTAg                                                                  >  2:2001384/1‑150 (MQ=255)
                                                                                 cTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAg                                                                  >  1:1228069/1‑150 (MQ=255)
                                                                                 cTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAg                                                                  <  1:920647/150‑1 (MQ=255)
                                                                                  tGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGc                                                                 >  1:819502/1‑150 (MQ=255)
                                                                                    gagaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTGAGCGa                                                               >  2:850387/1‑150 (MQ=255)
                                                                                    gagaAAAGCATTGCGTGGCCTCGTAATCGCTATATATACTATTATGTCGCCTGAAACCCACTTCGCTGTGGGTTATTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGa                                                               >  2:1765248/1‑150 (MQ=255)
                                                                                      gaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                             >  2:692737/1‑150 (MQ=255)
                                                                                            cATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAaataat                                                       <  2:558416/150‑1 (MQ=255)
                                                                                             aTTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATg                                                      <  1:50396/150‑1 (MQ=255)
                                                                                                   tGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACa                                                <  1:1358637/150‑1 (MQ=255)
                                                                                                       cTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCt                                            >  2:1079841/1‑150 (MQ=255)
                                                                                                        tCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTg                                           >  1:1920781/1‑150 (MQ=255)
                                                                                                          gTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAc                                         >  1:1538742/1‑150 (MQ=255)
                                                                                                          gTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAc                                         <  1:405245/150‑1 (MQ=255)
                                                                                                            aaTCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATTAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcct                                       >  2:501050/1‑150 (MQ=255)
                                                                                                                gCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGCCctcctc                                   >  2:62447/1‑150 (MQ=255)
                                                                                                                 cTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCa                                  >  1:159854/1‑150 (MQ=255)
                                                                                                                   atatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCAtt                                <  2:1973077/150‑1 (MQ=255)
                                                                                                                    tatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCAttt                               <  2:1337528/150‑2 (MQ=255)
                                                                                                                       ctactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGta                            >  1:647106/1‑149 (MQ=255)
                                                                                                                       ctactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGta                            >  2:1644967/1‑149 (MQ=255)
                                                                                                                       ctactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTc                            <  2:1540142/150‑1 (MQ=255)
                                                                                                                       ctactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTc                            >  1:655587/1‑150 (MQ=255)
                                                                                                                       atactaTTATGACGCCTGAAACCCACTTCTCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAATTGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTc                            <  2:1617276/149‑1 (MQ=255)
                                                                                                                            attatGTCGCCTGACACCCACTTCGCGGGGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAAt                       <  2:1192830/150‑1 (MQ=255)
                                                                                                                                tGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGaag                   >  2:1949267/1‑148 (MQ=255)
                                                                                                                                tGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACg                   <  1:405269/150‑1 (MQ=255)
                                                                                                                                tGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACg                   >  2:276065/1‑150 (MQ=255)
                                                                                                                                   cGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATTGCACTGGAGCGAAAAATAATGCTGGCAATCTGACCTCCTGATTAGTCAAAATGACGCAg                >  1:1882053/1‑150 (MQ=255)
                                                                                                                                   cGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGCCCTCCTCATTAGTCAAAATGACGCAg                >  2:1927704/1‑150 (MQ=255)
                                                                                                                                    gCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTAGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgg               >  1:1073915/1‑150 (MQ=255)
                                                                                                                                    acaaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgg               >  1:545784/6‑150 (MQ=255)
                                                                                                                                    acaaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTATAAAAAATGGCAGAGCAAACCTCGCGTCACGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCATCATTAGACAAAATGACGCAgg               >  1:731914/6‑150 (MQ=255)
                                                                                                                                    acaaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGGTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgg               >  1:1948398/6‑150 (MQ=255)
                                                                                                                                    acaaaaaaCCCACCGCGAAGTGGGTTTGTTGTTGGCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCGCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGCCAATCTGACCTCCTCATTAGTCAAAATGACGCAgg               >  1:212176/6‑150 (MQ=255)
                                                                                                                                     caaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGc              >  1:847129/5‑150 (MQ=255)
                                                                                                                                      cTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt             >  1:554397/1‑150 (MQ=255)
                                                                                                                                      cTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt             >  1:453399/1‑150 (MQ=255)
                                                                                                                                         aaaCCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTggg          <  1:1868816/150‑1 (MQ=255)
                                                                                                                                          aaCCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAACTAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTgggg         <  2:1266018/150‑1 (MQ=255)
                                                                                                                                               aCTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCTCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAAct    <  1:106597/150‑1 (MQ=255)
                                                                                                                                               aCTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAAct    <  1:539564/150‑1 (MQ=255)
                                                                                                                                               aCTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCACGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAAct    >  2:1794814/1‑150 (MQ=255)
                                                                                                                                                 ttCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAActgc  <  2:965654/150‑1 (MQ=255)
                                                                                                                                                 ttCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAActgc  <  1:909610/150‑1 (MQ=255)
                                                                                                                                                 ttCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGGCAAAATGACGCAGGCTGGGGAAActgc  <  2:1061199/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGC  >  NC_000913/1426272‑1426566

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