Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,446 C→T 6.1% A793V (GCG→GTG)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,4460CT 6.1% 345.0 / 7.0 113A793V (GCG→GTG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (55/51);  new base T (2/5);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 113 total aligned reads displayed.

GGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGA  >  NC_000913/1660299‑1660593
                                                                                                                                                   |                                                                                                                                                    
ggCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGa                                                                                                                                                    <  2:1209223/150‑1 (MQ=255)
  cTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑g                                                                                                                                                 >  2:1619887/1‑150 (MQ=255)
   tACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑gg                                                                                                                                                <  1:1048587/150‑1 (MQ=255)
    aacAGCCAAAGTAACACCACGTTTTCTCCCTGGAGTTAGCGCTATGGTCCAGGGAGCCTGGCAGGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCTCTCTGCGCCCTTCACCACTGGCGAA‑ggg                                                                                                                                               <  2:1429500/148‑1 (MQ=255)
    aCCATCCAAAGTCCCACCACGTATTATCCCTGGAGTTCGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCATTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑ggg                                                                                                                                               <  1:250888/150‑1 (MQ=255)
    aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑ggg                                                                                                                                               >  1:1316581/1‑150 (MQ=255)
     ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAA‑gggc                                                                                                                                              <  1:1041791/150‑2 (MQ=255)
     ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAA‑gggc                                                                                                                                              <  2:1701739/150‑2 (MQ=255)
     ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAA‑gggc                                                                                                                                              <  2:28520/150‑2 (MQ=255)
      cAGCCAAAGTAACACGACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGGACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAA‑gggcg                                                                                                                                             <  2:201649/150‑3 (MQ=255)
       aGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGaaa                                                                                                                                            >  1:141352/1‑150 (MQ=255)
         ccAAAGTAACACCTCGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAAcc                                                                                                                                          <  1:1191290/150‑1 (MQ=255)
         ccAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATACGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAccc                                                                                                                                          >  1:1083119/1‑150 (MQ=255)
          cAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAAccc                                                                                                                                         >  2:428204/1‑150 (MQ=255)
            aaGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCCCCACTGGCGAA‑GGGAAGCCCgc                                                                                                                                       >  1:1266339/1‑150 (MQ=255)
             aGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCCCGAGGCCAATATGTCTGGCGCTAAAATCGACCATGGCTGCTGAGTGAATACTCTAACAACTCTGCGCCCATCACCACTGGCGAA‑GGTAAACCCgca                                                                                                                                      >  2:311660/1‑150 (MQ=255)
             aGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCgca                                                                                                                                      <  1:848921/150‑1 (MQ=255)
             aGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGCCTGTGTGAATACGCTACCCACTCTGCGCCCTTCCCCCCTGGCGAA‑GGGAAACCCgca                                                                                                                                      >  2:502792/1‑150 (MQ=255)
                aaCACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCgcagcc                                                                                                                                   >  1:1302442/1‑149 (MQ=255)
                 acacCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGcac                                                                                                                                  <  1:512612/150‑1 (MQ=255)
                  caccacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGcaca                                                                                                                                 >  1:1572322/1‑150 (MQ=255)
                     cacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTa                                                                                                                              >  1:568652/1‑150 (MQ=255)
                          aTTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTg                                                                                                                         >  2:864608/1‑150 (MQ=255)
                             ctcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTc                                                                                                                      <  1:1000729/150‑1 (MQ=255)
                                ccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGTCCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCgag                                                                                                                   <  1:39051/150‑1 (MQ=255)
                                   ggAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAATATCGACGATGGCGGCTGTGTGAATACGCTTACCACTCTGCGCCCTTCACCGCTGGCGAA‑GGGATACCCGCAGCTCACTAATCTGGTCGAGATc                                                                                                                >  1:1530084/1‑150 (MQ=255)
                                          gcgcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaaa                                                                                                         <  1:1077786/150‑1 (MQ=255)
                                          gcgcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaaa                                                                                                         >  2:317777/1‑150 (MQ=255)
                                             cTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACAACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAatat                                                                                                      >  2:1012595/1‑150 (MQ=255)
                                              tATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAAtata                                                                                                     <  2:933449/150‑1 (MQ=255)
                                                 gggCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAattataa                                                                                                   >  2:712917/1‑145 (MQ=255)
                                                 gggCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAcc                                                                                                  >  2:773931/1‑150 (MQ=255)
                                                 gggCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAcc                                                                                                  >  1:733237/1‑150 (MQ=255)
                                                    ccAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCAc                                                                                               <  1:1432326/150‑1 (MQ=255)
                                                       gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAACTATAACCCACGAc                                                                                            >  2:444779/1‑150 (MQ=255)
                                                          aGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAAc                                                                                         >  1:61874/1‑150 (MQ=255)
                                                          aGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAAc                                                                                         <  1:1709951/150‑1 (MQ=255)
                                                            ccTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTACCCACTCTGCGCCCTTCACCCCTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGCCAACCa                                                                                       >  2:1581235/1‑150 (MQ=255)
                                                            ccTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCa                                                                                       <  1:499269/150‑1 (MQ=255)
                                                            ccTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCa                                                                                       >  1:551907/1‑150 (MQ=255)
                                                            ccTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCT‑GTCGAGATCGAAAAAATATAACCCACGACAACCAt                                                                                      >  2:1510522/1‑150 (MQ=255)
                                                                 cACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAAt                                                                                  <  2:1079672/150‑1 (MQ=255)
                                                                      ggCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGCTCGAAAAAATATAACCCACGACAACCATAAATTCTgg                                                                             >  2:690146/1‑150 (MQ=255)
                                                                       gCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGc                                                                            <  2:1692740/150‑1 (MQ=255)
                                                                        ccAATATGTCTGGCGATAAATTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGTA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGTTCGAAAAATTATTTCCGACGACACCCTTAATTTCTGGCa                                                                           >  2:250591/1‑150 (MQ=255)
                                                                        ccAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCCCTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCCCGACAACCATAAATTCTGGCa                                                                           >  1:94611/1‑150 (MQ=255)
                                                                        ccAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTTGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCa                                                                           <  1:297350/150‑1 (MQ=255)
                                                                        ccAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATACCCCACGACAACCATAAATTcagtca                                                                           >  1:236373/1‑145 (MQ=255)
                                                                         aaaTATGTCTGTCGATAATATCGACCATGTCGGATGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTTGCGAA‑GGGATACCCGCAGCACACTAATCTGGTCGAGATCGATAAAATATAACCCACGACAACCATAAATTCTGGCAt                                                                          <  2:882945/149‑1 (MQ=255)
                                                                          aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATg                                                                         >  2:1513262/1‑150 (MQ=255)
                                                                              tGTCTGGCTATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAt                                                                     <  1:169582/150‑1 (MQ=255)
                                                                               gTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTAACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtt                                                                    >  2:564595/1‑150 (MQ=255)
                                                                                tCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAttt                                                                   >  2:1700278/1‑150 (MQ=255)
                                                                                 cTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                                  <  2:733237/150‑1 (MQ=255)
                                                                                   ggCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAATCCCGCAGCACACTAATCTGGTCGAGATCGAAAAACTATTTCCCACGACAACCTTAAATTCTGGCTTGACTGtttgt                                                                >  1:1049849/1‑150 (MQ=255)
                                                                                   ggCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                                >  1:1023999/1‑150 (MQ=255)
                                                                                   ggCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                                >  2:771933/1‑150 (MQ=255)
                                                                                   ggCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                                >  2:738515/1‑150 (MQ=255)
                                                                                    gCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtt                                                               <  1:1275143/150‑1 (MQ=255)
                                                                                     cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgttt                                                              <  2:737223/150‑1 (MQ=255)
                                                                                     cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCT‑GTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtttg                                                             <  1:1510522/150‑1 (MQ=255)
                                                                                      gATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtttg                                                             <  2:233110/150‑1 (MQ=255)
                                                                                      gATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgtttg                                                             <  2:1654035/150‑1 (MQ=255)
                                                                                       aTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                            >  1:179814/1‑150 (MQ=255)
                                                                                       aTAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                            >  1:1314720/1‑150 (MQ=255)
                                                                                        tAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaa                                                           <  1:28339/150‑1 (MQ=255)
                                                                                        tAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaa                                                           <  1:104886/150‑1 (MQ=255)
                                                                                               gACCATGGCGGCTGTGTGAATACTCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCTAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAAt                                                    <  1:325112/150‑1 (MQ=255)
                                                                                                    tGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                               <  1:502792/150‑1 (MQ=255)
                                                                                                    tGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                               <  1:600859/150‑1 (MQ=255)
                                                                                                     ggcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGa                                              >  2:1147632/1‑150 (MQ=255)
                                                                                                        ggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTa                                           <  2:705131/150‑1 (MQ=255)
                                                                                                         gcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTaa                                          >  1:1237471/1‑150 (MQ=255)
                                                                                                         gcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTaa                                          >  1:418616/1‑150 (MQ=255)
                                                                                                           tgtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatg                                        <  2:111012/150‑1 (MQ=255)
                                                                                                               tgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAAtgatga                                    >  2:1545093/1‑150 (MQ=255)
                                                                                                                  aTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAC‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATACCCCACGACCACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGataa                                 >  1:1166523/1‑150 (MQ=255)
                                                                                                                    aCGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATc                               <  1:212233/150‑1 (MQ=255)
                                                                                                                      gCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCa                             <  2:852088/150‑1 (MQ=255)
                                                                                                                       cTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCaa                            >  2:1306897/1‑149 (MQ=255)
                                                                                                                       cTAACCACTCTGCGCCCTTCACCACTGGCGAA‑CGGAAACCAGCAGCAAACTAATCTTGTCTAGATCTAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCAt                            <  2:1133227/150‑1 (MQ=255)
                                                                                                                         aaCCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAc                          <  2:962921/150‑1 (MQ=255)
                                                                                                                         aaCCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATAc                          <  1:708868/150‑1 (MQ=255)
                                                                                                                               gggtttcccTTCGCCAGTGGTGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACaga                    >  2:525454/7‑150 (MQ=255)
                                                                                                                                ctGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACagag                   <  1:565271/150‑1 (MQ=255)
                                                                                                                                  gcgcCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgc                 >  1:1617341/1‑150 (MQ=255)
                                                                                                                                    tcccTTCGCCAGTGGTGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGgcgc               >  2:1294689/2‑150 (MQ=255)
                                                                                                                                       cTTCACCACTGGCGAA‑GGGAAACCCGCCGCTCACTAATCTGGTCGAGATCGAAAAAATATAACCCACGCCAACCATAAATTCTGGCATGACATTTTGTTTGAATAGCAATAAGTGAGTAATGATGAAAATCCATACCCCAGAGGCGCTCa            >  1:95639/1‑150 (MQ=255)
                                                                                                                                       cTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCa            <  2:1061496/150‑1 (MQ=255)
                                                                                                                                       cTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCa            >  1:1211438/1‑150 (MQ=255)
                                                                                                                                        ttCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCAt           <  1:1669742/150‑1 (MQ=255)
                                                                                                                                        ttCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCAt           <  2:916375/150‑1 (MQ=255)
                                                                                                                                          caccacTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATATGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGa         >  2:927241/1‑150 (MQ=255)
                                                                                                                                            ccacTGTCGAA‑GTGAAACCCGCAGCACACTAATCTGGTCGAGATCGAATAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAATAGCTATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAg       <  1:336653/150‑1 (MQ=255)
                                                                                                                                            ccacTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAg       >  2:1475426/1‑150 (MQ=255)
                                                                                                                                            ccacTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAg       >  1:1241428/1‑150 (MQ=255)
                                                                                                                                              acTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGc     <  2:1083119/150‑1 (MQ=255)
                                                                                                                                              acTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAACCCATACCACAGAGGCGCTCATGAAGGc     <  1:304831/150‑1 (MQ=255)
                                                                                                                                              acTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGAGAGTAATGATGAAAATCCATACCACAGAGGCGATCATGAAGGc     >  2:1408626/1‑150 (MQ=255)
                                                                                                                                                 ggCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGa  >  1:1265260/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAA‑GGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGA  >  NC_000913/1660299‑1660593

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: ‑