Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,485,384 A→T 6.1% intergenic (+25/+30) prkB → / ← yhfA putative phosphoribulokinase/OsmC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,485,3840AT 6.1% 475.1 / 13.8 147intergenic (+25/+30)prkB/yhfAputative phosphoribulokinase/OsmC family protein
Reads supporting (aligned to +/- strand):  ref base A (74/64);  new base T (6/3);  total (80/67)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 150 total aligned reads displayed.

TGCAAGGCTCATTCATTTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCTGCGTCTTTCAG  >  NC_000913/3485235‑3485533
                                                                                                                                                     |                                                                                                                                                     
     ggCTCATTCATTTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATc                                                                                                                                                  >  2:1084405/1‑150 (MQ=255)
     ggCTCATTCATTTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATc                                                                                                                                                  >  2:2240175/1‑150 (MQ=255)
      gCTCATTCATTTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATcc                                                                                                                                                 >  2:1824917/1‑150 (MQ=255)
       cTCATTCATTTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCg                                                                                                                                                <  1:876450/150‑1 (MQ=255)
       cTCATTCATTTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCg                                                                                                                                                >  2:1483121/1‑150 (MQ=255)
        acattcataTCCCACATCAATACCTTAGTGGTACCGGGCGGCAAAATGGGACTCGCAAGCGAAATAATAATGCAGCCCCTGCTGCAACCATTGATCCAACGAAACAAAATCGACAAACTCTGCTATTACGCCCGATAAAATACTAACCgg                                                                                                                                               >  2:969391/2‑150 (MQ=255)
          attcattTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAGTATTTTATCCgacg                                                                                                                                             <  2:2309772/150‑4 (MQ=255)
          attcattTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAGTATTTTATCCGGcg                                                                                                                                             <  1:2105358/150‑2 (MQ=255)
          attcattTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCt                                                                                                                                             >  1:1265506/1‑150 (MQ=255)
           ttcattTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCtt                                                                                                                                            >  2:312840/1‑150 (MQ=255)
            tcattTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATTGTTCTGGCAATGGAATTAATTATGCTTCCGCTGGTTCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCAGGATAAAATACTATCCGGCTTc                                                                                                                                           <  1:337034/150‑1 (MQ=255)
            tcattTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTc                                                                                                                                           <  1:1685678/150‑1 (MQ=255)
             cattTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTca                                                                                                                                          >  1:1463126/1‑150 (MQ=255)
                 tCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAAc                                                                                                                                      >  2:451549/1‑150 (MQ=255)
                  cccACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGATGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACg                                                                                                                                     >  2:1438430/1‑150 (MQ=255)
                    cacaTCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACggg                                                                                                                                   >  1:285574/1‑150 (MQ=255)
                         cAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGt                                                                                                                              <  2:1525132/150‑1 (MQ=255)
                             aCGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAg                                                                                                                          <  1:2282091/150‑1 (MQ=255)
                               gTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTc                                                                                                                        <  1:2042036/150‑1 (MQ=255)
                                  aGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCAcc                                                                                                                     >  2:841065/1‑149 (MQ=255)
                                     ggTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCgg                                                                                                                  >  1:813948/1‑150 (MQ=255)
                                     ggTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCgg                                                                                                                  >  1:443471/1‑150 (MQ=255)
                                          cGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAAcc                                                                                                             <  1:1824917/150‑1 (MQ=255)
                                          cGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAAcc                                                                                                             <  2:410589/150‑1 (MQ=255)
                                            ggcggcAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCTGCAACCAc                                                                                                           >  1:315088/1‑150 (MQ=255)
                                             gcggcAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACt                                                                                                          >  2:1193673/1‑150 (MQ=255)
                                              cggcAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACtt                                                                                                         >  1:100902/1‑150 (MQ=255)
                                                gcAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCa                                                                                                       <  2:2012384/150‑1 (MQ=255)
                                                 cAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCAt                                                                                                      <  2:1923788/150‑1 (MQ=255)
                                                  aaaaTGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATa                                                                                                     >  2:1111841/1‑150 (MQ=255)
                                                   aaaTGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATAc                                                                                                    >  2:347517/1‑150 (MQ=255)
                                                   aaaTGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATAc                                                                                                    >  1:551417/1‑150 (MQ=255)
                                                   aaaTGGGGCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGAATGATGGAAGGACAGAAACTCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATAc                                                                                                    <  2:1588713/150‑1 (MQ=255)
                                                    aaTGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACg                                                                                                   <  2:1578734/150‑1 (MQ=255)
                                                    aaTGGGTCTGGCAATGGAATTAATTATGCAGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACg                                                                                                   <  2:2145256/150‑1 (MQ=255)
                                                        ggTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTGACTCTGCTATTACGCCGGATAAAATACTATCCTGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACgagtg                                                                                               <  2:1741697/150‑1 (MQ=255)
                                                         gTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACgagtga                                                                                              >  2:179512/1‑150 (MQ=255)
                                                           cTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACgagtgagt                                                                                            <  2:1514015/150‑1 (MQ=255)
                                                           cTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACgagtgagt                                                                                            >  2:183262/1‑150 (MQ=255)
                                                            tGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTa                                                                                           >  2:2026081/1‑150 (MQ=255)
                                                             ggCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATGCGGCTTCACAACGGGATTGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTaa                                                                                          >  1:2089052/1‑150 (MQ=255)
                                                                aaTGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAatat                                                                                       <  1:1111841/150‑1 (MQ=255)
                                                                    gAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCAc                                                                                   >  2:136848/1‑150 (MQ=255)
                                                                     aattaattATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCAcc                                                                                  >  2:80188/1‑150 (MQ=255)
                                                                     aattaattATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCAcc                                                                                  <  2:263582/150‑1 (MQ=255)
                                                                      attaattaTGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTACCTCTGCTATTACGCCGGATAAAATAATATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATAAGAGTGAGTAATATTCAccc                                                                                 >  2:96647/1‑149 (MQ=255)
                                                                       ttaattaTGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGc                                                                                <  1:1375753/150‑1 (MQ=255)
                                                                       ttaattaTGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGc                                                                                <  1:432197/150‑1 (MQ=255)
                                                                              tgctgcCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTcc                                                                         <  2:1463126/150‑1 (MQ=255)
                                                                                ctgcCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAg                                                                       <  1:347517/150‑1 (MQ=255)
                                                                                ctgcCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAg                                                                       <  1:1387539/150‑1 (MQ=255)
                                                                                ctgcCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAg                                                                       <  1:627683/150‑1 (MQ=255)
                                                                                  gcCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGATTTTTCCAGCa                                                                     <  2:1912991/150‑1 (MQ=255)
                                                                                     gCTGGTGCAACGATTGATGGAAGGAAAGATAATCGAGTAACTCTGGTATTACGCCGGATATAATAGTATCCGGCTTATCAACGGGATAGTTAAGTCACGCGGCAACCACTTCATTCGAGTGAGTAATATTCACGGCTTTTTCCAGCATCa                                                                  <  2:1624263/150‑1 (MQ=255)
                                                                                        ggTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACg                                                               >  1:1667060/1‑150 (MQ=255)
                                                                                                 aTTGATGGAAGGAAAGAAAATCGAGTAACTCTTCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGATTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAAc                                                      <  1:142757/150‑1 (MQ=255)
                                                                                                    gATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGGATCAACGCCACTGAACAat                                                   >  1:1035769/1‑150 (MQ=255)
                                                                                                      tGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAatat                                                 >  1:1119236/1‑150 (MQ=255)
                                                                                                      tGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAatat                                                 <  2:816049/150‑1 (MQ=255)
                                                                                                           ggaaAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTtctc                                            >  1:1847206/1‑150 (MQ=255)
                                                                                                           ggaaAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTtctc                                            <  1:979950/150‑1 (MQ=255)
                                                                                                            gaaagaaaATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCg                                           >  2:32363/1‑150 (MQ=255)
                                                                                                             aaagaaaATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCgg                                          >  1:1506788/1‑150 (MQ=255)
                                                                                                             aaagaaaATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCgg                                          <  1:1941893/150‑1 (MQ=255)
                                                                                                              aagaaaATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGc                                         <  2:643768/150‑1 (MQ=255)
                                                                                                               agaaaATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCa                                        >  2:2424519/1‑150 (MQ=255)
                                                                                                                 aaaaTCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGCACAATATTTCTCGGCaga                                      >  2:1140170/1‑150 (MQ=255)
                                                                                                                  aaaTCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCagag                                     <  2:773158/150‑1 (MQ=255)
                                                                                                                  aaaTCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCagag                                     >  2:1649402/1‑150 (MQ=255)
                                                                                                                     tCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAgagaga                                  <  2:1425759/150‑1 (MQ=255)
                                                                                                                      cGAGTAACTCTGCTATTACGCCTGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGAt                                 <  1:1084405/150‑1 (MQ=255)
                                                                                                                       gAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATc                                <  1:26069/150‑1 (MQ=255)
                                                                                                                       gAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATc                                >  1:597486/1‑150 (MQ=255)
                                                                                                                        aGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCa                               <  1:748248/150‑1 (MQ=255)
                                                                                                                           aaCTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAAcc                            <  1:2158245/150‑1 (MQ=255)
                                                                                                                           aaCTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAAcc                            <  2:285574/150‑1 (MQ=255)
                                                                                                                           aaCTCTGCTATTACGCCGGATAAAATACTATCAGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAAcc                            >  2:2355522/1‑150 (MQ=255)
                                                                                                                             ctctGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGc                          >  1:181365/1‑150 (MQ=255)
                                                                                                                              tctGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCa                         <  1:251698/150‑1 (MQ=255)
                                                                                                                                 gCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCAcgc                      >  2:1886384/1‑150 (MQ=255)
                                                                                                                                 gCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCAcgc                      <  1:1362762/150‑1 (MQ=255)
                                                                                                                                   tATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACgcgc                    >  1:1923163/1‑150 (MQ=255)
                                                                                                                                   tATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACgcgc                    >  2:2386384/1‑150 (MQ=255)
                                                                                                                                     tgaagCCGGATAGTATTTTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCaa                  >  1:83085/5‑150 (MQ=255)
                                                                                                                                      gaagCCGGATAGTATTTTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAAc                 >  1:79702/4‑150 (MQ=255)
                                                                                                                                       aCGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAAcc                >  1:959758/1‑150 (MQ=255)
                                                                                                                                       aCGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCAAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAAcc                >  2:2392283/1‑150 (MQ=255)
                                                                                                                                        cGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCg               <  1:487286/150‑1 (MQ=255)
                                                                                                                                        agCCGGATAGTATTTTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCg               >  1:1765213/2‑150 (MQ=255)
                                                                                                                                        agCCGGATAGTATTTTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCg               >  1:298807/2‑150 (MQ=255)
                                                                                                                                         gCCGGATAGTATTTTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGc              >  2:2066753/1‑150 (MQ=255)
                                                                                                                                         gCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGc              >  1:2304155/1‑150 (MQ=255)
                                                                                                                                          ccGGATAAAATACTATCCGGCTTCCCAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCt             >  2:929815/1‑150 (MQ=255)
                                                                                                                                          ccGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATTCGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCt             <  2:1265506/150‑1 (MQ=255)
                                                                                                                                            ggATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCTGc           <  2:1119236/150‑1 (MQ=255)
                                                                                                                                                 aaaTACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCACCCGCTGCGTCtt      >  1:2204037/1‑150 (MQ=255)
                                                                                                                                                  aaTACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCTGCGTCttt     >  1:482491/1‑150 (MQ=255)
                                                                                                                                                   aTACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCTGCGTCTTTc    <  1:2026081/150‑1 (MQ=255)
                                                                                                                                                    tACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCTGCGTCTTTCa   >  2:901422/1‑150 (MQ=255)
                                                                                                                                                     aCTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCTGCGTCTTTCAg  <  2:171096/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCAAGGCTCATTCATTTCCCACATCAATACGTTAGTGGTACCGGGCGGCAAAATGGGTCTGGCAATGGAATTAATTATGCTGCCGCTGGTGCAACGATTGATGGAAGGAAAGAAAATCGAGTAACTCTGCTATTACGCCGGATAAAATACTATCCGGCTTCACAACGGGATAGTTAAGTCACGCGGCAACCACTTCATACGAGTGAGTAATATTCACCGCTTTTTCCAGCATCAACGCCACTGAACAATATTTCTCGGCAGAGAGATCAACCGCACGCGCAACCGCTGCGTCTTTCAG  >  NC_000913/3485235‑3485533

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