Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,086,641 C→A 9.1% intergenic (+239/‑65) yqgC → / → metK uncharacterized protein/S‑adenosylmethionine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,086,6410CA 9.1% 298.2 / 2.9 123intergenic (+239/‑65)yqgC/metKuncharacterized protein/S‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base C (45/61);  new base A (13/3);  total (59/64)
Fisher's exact test for biased strand distribution p-value = 5.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTA  >  NC_000913/3086492‑3086789
                                                                                                                                                     |                                                                                                                                                    
ccACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTCGTAGGATc                                                                                                                                                      >  1:105504/1‑150 (MQ=255)
   cGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTTAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGAAGGAACCGc                                                                                                                                                   >  1:277726/1‑150 (MQ=255)
    gTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCt                                                                                                                                                  <  1:590256/150‑1 (MQ=255)
        ttAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGAACTAAAATAGCCATCCAGATGTTAATACATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGAGGTAAGATCCGCTAccc                                                                                                                                              >  1:629479/1‑149 (MQ=255)
          aaCAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGCTAAGCTAAcaca                                                                                                                                            >  1:1302946/1‑150 (MQ=255)
              gAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGaaa                                                                                                                                        >  1:1521510/1‑150 (MQ=255)
              gAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGaaa                                                                                                                                        >  2:448876/1‑150 (MQ=255)
              gAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGaaa                                                                                                                                        <  1:1363306/150‑1 (MQ=255)
              gAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATACGCTACCACAGaaa                                                                                                                                        >  1:1471965/1‑150 (MQ=255)
                 tGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATc                                                                                                                                     <  1:880473/150‑1 (MQ=255)
                 tGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATc                                                                                                                                     >  1:1292469/1‑150 (MQ=255)
                  gagaCACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATcc                                                                                                                                    <  2:816380/150‑1 (MQ=255)
                    gaCACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCcac                                                                                                                                  >  2:913097/1‑150 (MQ=255)
                     acacGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCcaca                                                                                                                                 <  2:1691007/150‑1 (MQ=255)
                          aTTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACacaaca                                                                                                                            <  2:78524/150‑1 (MQ=255)
                          aTTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATACGCTTCCACAGAAAATCCACacaaca                                                                                                                            >  1:192436/1‑150 (MQ=255)
                             cAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGtt                                                                                                                         >  2:223848/1‑150 (MQ=255)
                             cAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGGAGGATACGCTACCACAGAAAATCCACACAACAGtt                                                                                                                         >  1:1318259/1‑150 (MQ=255)
                              aaaaaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTACCACTCAGACTGCCAGTGTTTTTAACATGCAGAGTCGTGGTAGGATCCCCTACCACAGAAAAACCACacaacacgtt                                                                                                                        >  2:1628665/1‑146 (MQ=255)
                              aaaaaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGttt                                                                                                                        >  1:1104045/1‑150 (MQ=255)
                              aaaaaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGttt                                                                                                                        >  2:546910/1‑150 (MQ=255)
                               aaaaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTg                                                                                                                       >  2:45064/1‑150 (MQ=255)
                               aaaaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTg                                                                                                                       <  2:379768/150‑1 (MQ=255)
                                 aaaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAg                                                                                                                     <  2:1572717/150‑1 (MQ=255)
                                  aaaaGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGc                                                                                                                    <  1:321940/150‑1 (MQ=255)
                                      gTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAAc                                                                                                                <  2:202248/150‑1 (MQ=255)
                                       tGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCAACCACACAAAATCCACACAACAGTTTGAGCaaaaa                                                                                                               >  1:1547630/1‑148 (MQ=255)
                                           aaTAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATTTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTAGCACAGAAAATCCACACAACAGTTTGAGCTAACCAAAt                                                                                                           <  1:575134/150‑1 (MQ=255)
                                                  tGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCttt                                                                                                    <  2:1464967/150‑1 (MQ=255)
                                                  tGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCttt                                                                                                    <  2:1604472/150‑1 (MQ=255)
                                                   gaagaaTTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTa                                                                                                   <  2:212145/150‑1 (MQ=255)
                                                   gaagaaTTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTACGATCCGCTACCACAGAAAATCCACCAAAAAGTTTAAGCAAACCAAATTCTCTTTa                                                                                                   >  1:446125/1‑150 (MQ=255)
                                                   gaagaaTTGAACTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTTTTTTTAACCTGCAGAGTCGTGGTAGGATCATCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTa                                                                                                   <  1:316799/150‑1 (MQ=255)
                                                       aaTTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTg                                                                                               <  1:155465/150‑1 (MQ=255)
                                                       aaTTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTg                                                                                               <  2:1479424/150‑1 (MQ=255)
                                                           gACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGCCTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGatat                                                                                           >  2:664409/1‑150 (MQ=255)
                                                            aCCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATAtt                                                                                          >  2:39130/1‑150 (MQ=255)
                                                             ccTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTa                                                                                         >  1:918278/1‑150 (MQ=255)
                                                             ccTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTa                                                                                         <  1:290771/150‑1 (MQ=255)
                                                               tAAAATAGCCATCCAGATGTTAATCCACCCATACCGATTAACACTCAGACTGCAAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTaaa                                                                                       <  1:1336001/150‑1 (MQ=255)
                                                                     cgCCATCCATATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGc                                                                                 <  1:674945/149‑1 (MQ=255)
                                                                     aGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGTATATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGc                                                                                 <  2:1292469/150‑1 (MQ=255)
                                                                     aGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGc                                                                                 <  2:1568774/150‑1 (MQ=255)
                                                                     aGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGc                                                                                 <  2:1480168/150‑1 (MQ=255)
                                                                     aGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGc                                                                                 <  1:1028976/150‑1 (MQ=255)
                                                                      gCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCa                                                                                >  2:1013615/1‑150 (MQ=255)
                                                                      gCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATACACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCa                                                                                >  1:670511/1‑150 (MQ=255)
                                                                       ccATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCaa                                                                               >  1:1031615/1‑150 (MQ=255)
                                                                        cATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATACGCTACCACACAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCaaa                                                                              >  1:1526331/1‑150 (MQ=255)
                                                                          tCCAGATGTTAATCAATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGAACCGCTACCACATAAAAAACACAAAACAGTTTGAGATAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAac                                                                            >  1:927175/1‑150 (MQ=255)
                                                                           ccAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGAAATTAAATATGGCAAAaca                                                                           >  1:67437/1‑150 (MQ=255)
                                                                            cAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAacac                                                                          <  1:351203/150‑1 (MQ=255)
                                                                             aGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACAcc                                                                         <  2:105504/150‑1 (MQ=255)
                                                                             aGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACAcc                                                                         <  1:97971/150‑1 (MQ=255)
                                                                                tGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGGAGGATACGCGACCACACAAAATACAACCAAAACTTTTAGCTACACAAAATATATTTAAGTGATATTAAATATTGCCAAACAccaat                                                                      >  2:1407142/1‑147 (MQ=255)
                                                                                 gTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGAACCGCTACCAAATAAAATACCAAAAAAAGTTTGATCTAACCAAAATATATTTAAGTGATTTTAAATATTGCAAAACAACattt                                                                     >  2:358509/1‑150 (MQ=255)
                                                                                   tAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTa                                                                   <  2:1088202/150‑1 (MQ=255)
                                                                                    aatcaatcCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGCAAATCCACACAACCGGTTGACCTATCCAAATTCTCTTTAGGTGATATTATATATGGGAAAACACCTTTTTaa                                                                  >  2:449135/1‑149 (MQ=255)
                                                                                    aaTCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATACGCTACCACAGAAAATCCACACAACAGTTTGAGCAAACCAAATTCTCTTTAGGTGATATTAAAAATGGCAAAACAACTTTTTaa                                                                  >  1:426475/1‑149 (MQ=255)
                                                                                    aaTCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGAACCGCTACCACAGCAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTATAc                                                                  >  1:63670/1‑150 (MQ=255)
                                                                                       ccatccatACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTTCATAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTc                                                               <  2:714642/150‑1 (MQ=255)
                                                                                       ccatccatACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGt                                                                <  2:931029/149‑1 (MQ=255)
                                                                                       ccatccatACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTc                                                               <  1:913097/150‑1 (MQ=255)
                                                                                       ccatccatACCGATTAACACTCAGACTGACAGTGTTTTTAACGTTATGATTTTTTTTATGATACGGTTACAAAGTATATCTACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTc                                                               <  2:1547630/150‑1 (MQ=255)
                                                                                        catccatACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATTTCAAAACACCTTTTTACGTcc                                                              >  1:460430/1‑150 (MQ=255)
                                                                                           ccatACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAg                                                           <  2:1471965/150‑1 (MQ=255)
                                                                                             atACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTc                                                         <  2:277726/150‑1 (MQ=255)
                                                                                               aCCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCg                                                       <  2:538508/150‑1 (MQ=255)
                                                                                                ccGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGt                                                      >  1:82448/1‑150 (MQ=255)
                                                                                                  gATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCAAAGCAAATACACACAACAGTTTGAGCTAACCAAATTGTCTTTAGGTGAAATTAAAAAAGGCAAAACACCTTTTTAAGTACGAGTCCGtct                                                    >  2:716439/1‑150 (MQ=255)
                                                                                                     tAACACTCAGACTGCCAGTGTTTTTAACCTGCAGATTTTTTTTAGGATACGGTTACAAAGTATATCTACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTg                                                 <  2:1318259/150‑1 (MQ=255)
                                                                                                     tAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACACAAAAAACACACAACAGTTTTAGCTAACCAAATTAGCTTTAGGTGATATTAAATATTGCAAAACCACATTTTAAGTCCGAGTACGTATCTg                                                 >  2:1332072/1‑150 (MQ=255)
                                                                                                      aaCACTCCGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGAGCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGa                                                <  2:475967/150‑1 (MQ=255)
                                                                                                      aaCACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCACTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGa                                                <  1:1635804/150‑1 (MQ=255)
                                                                                                        cacTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAg                                              >  1:9016/1‑150 (MQ=255)
                                                                                                        cacTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAg                                              >  2:701611/1‑150 (MQ=255)
                                                                                                        aacTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTTGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAg                                              <  2:132074/149‑1 (MQ=255)
                                                                                                          cTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAggg                                            <  2:907279/150‑1 (MQ=255)
                                                                                                              gACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATc                                        <  1:1344786/150‑1 (MQ=255)
                                                                                                               aCTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATcc                                       <  1:913194/150‑1 (MQ=255)
                                                                                                                 tGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTg                                     <  2:1261713/150‑1 (MQ=255)
                                                                                                                  gCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCGTTAGGTGATATTAAATATGGCAAAACACCTTTTTAAGTCCGAGTCCGTCGATGAAGGGCATACTGa                                    >  1:379911/1‑150 (MQ=255)
                                                                                                                        ttttttAACCTGCAGAGTCGTTGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAAt                              <  2:782035/149‑1 (MQ=255)
                                                                                                                        gtTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGGCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAAt                              <  1:662309/150‑1 (MQ=255)
                                                                                                                        gtTTTTAACCTGCAGAGTCGTGGTAGGATACGCTACCACAGAAAATCCACAAAAAAATTTGAGATAACCAAATTCTCTTTTGGTGATATTAAAAATTGCAAAACACCATTTTACGTCCGAGTCGGTCTATGAAGGGCATCCTGACaaaaa                              >  1:840122/1‑149 (MQ=255)
                                                                                                                            ttAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCt                          <  2:1104045/150‑1 (MQ=255)
                                                                                                                            ttAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCt                          <  2:918278/150‑1 (MQ=255)
                                                                                                                               aCCTGCAGAGTCGTTGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGAc                       <  2:542490/150‑1 (MQ=255)
                                                                                                                                ccTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGAcc                      <  1:394995/150‑1 (MQ=255)
                                                                                                                                ccTGCAGAGTCGTGGTAGGATACCCTAACACAGAAAATCCACACAACAGTTTGAGCGAAACAAAATCTCATTAGGTGATATTAAATATGGAAAAAAACCATTTTACGTCCGACTCCGTCGAGGAACGGCATCCTGACAAAATTTCTGaac                      >  2:1627713/1‑148 (MQ=255)
                                                                                                                                 cTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGCGCTAACCAAATTATCTTTAGGTGATATTAAAAATGGGAAAAAACAGTTTTTCTTCCGCGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCa                     >  2:223537/1‑150 (MQ=255)
                                                                                                                                        gTCGTGGTAGGATCCGCTACCACAGATAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCt              >  2:703010/1‑150 (MQ=255)
                                                                                                                                        gTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCt              >  2:1696760/1‑150 (MQ=255)
                                                                                                                                          cGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGa            <  1:45064/150‑1 (MQ=255)
                                                                                                                                               tAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCg       <  1:448876/150‑1 (MQ=255)
                                                                                                                                                aGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCAGAGTCCGTCTCTGCAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGt      >  2:75620/1‑150 (MQ=255)
                                                                                                                                                 ggATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGtt     >  2:1162441/1‑150 (MQ=255)
                                                                                                                                                  gATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGttt    <  1:223848/150‑1 (MQ=255)
                                                                                                                                                  gATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGttt    >  1:1436331/1‑150 (MQ=255)
                                                                                                                                                    tCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGGCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTa  <  1:39130/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCACGTTTTTAACAGAATGAGACACGATTCAAAAAAAAGTGGAAATAGGGTGAAGAATTGACCTAAAATAGCCATCCAGATGTTAATCCATCCATACCGATTAACACTCAGACTGCCAGTGTTTTTAACCTGCAGAGTCGTGGTAGGATCCGCTACCACAGAAAATCCACACAACAGTTTGAGCTAACCAAATTCTCTTTAGGTGATATTAAATATGGCAAAACACCTTTTTACGTCCGAGTCCGTCTCTGAAGGGCATCCTGACAAAATTGCTGACCAAATTTCTGATGCCGTTTTA  >  NC_000913/3086492‑3086789

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: