Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,555,781 G→C 5.1% intergenic (‑112/+45) maeA ← / ← sra malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme/stationary‑phase‑induced ribosome‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,7810GC 5.1% 353.8 / 9.8 117intergenic (‑112/+45)maeA/sramalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme/stationary‑phase‑induced ribosome‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (63/48);  new base C (3/3);  total (66/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-02

Only 100 of 121 total aligned reads displayed.

TATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCAGTC  >  NC_000913/1555633‑1555929
                                                                                                                                                    |                                                                                                                                                    
   aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAgg                                                                                                                                                  >  2:238754/1‑150 (MQ=255)
    aaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATTAAAACAAAAACACCATTGCGCAGGc                                                                                                                                                 >  1:344338/1‑150 (MQ=255)
         aaCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATgg                                                                                                                                            >  1:822062/1‑150 (MQ=255)
          aCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGgt                                                                                                                                           <  1:4495/150‑1 (MQ=255)
            gCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGgtgt                                                                                                                                         >  1:831223/1‑150 (MQ=255)
            gCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGgtgt                                                                                                                                         <  2:1429613/150‑1 (MQ=255)
             cTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGtt                                                                                                                                        >  2:1631323/1‑150 (MQ=255)
              tGTTTTTTTGTGTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGTGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGttt                                                                                                                                       <  1:311113/150‑1 (MQ=255)
               gTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAAAAAAAACACCATTGCGCAGGCAATGGTGTTTa                                                                                                                                      >  2:123482/1‑150 (MQ=255)
               gTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGAAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTa                                                                                                                                      >  1:222361/1‑150 (MQ=255)
               gTTTTTTTGTATTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACAGCCAGGGCGTCGGCCACAAGTATAAATAAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGtata                                                                                                                                      >  1:418341/1‑150 (MQ=255)
                 ttttttGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGtttttg                                                                                                                                    <  2:533761/150‑4 (MQ=255)
                    tttgttgttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTTCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAACACACCATTGCGCAGGCAATGGTGTTTAATCGt                                                                                                                                 <  1:202918/150‑1 (MQ=255)
                       gtttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAt                                                                                                                              <  2:1885068/150‑1 (MQ=255)
                        tttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAtt                                                                                                                             >  1:1999273/1‑150 (MQ=255)
                        tttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAtt                                                                                                                             <  2:781355/150‑1 (MQ=255)
                        tttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAtt                                                                                                                             <  1:625586/150‑1 (MQ=255)
                                   aTGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATg                                                                                                                  <  2:734897/150‑1 (MQ=255)
                                       cactcactCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGTACTGATGGTTa                                                                                                              >  1:707213/1‑150 (MQ=255)
                                       cactcactCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTa                                                                                                              >  1:1620816/1‑150 (MQ=255)
                                           cactcTTTTTTGAATATCCATCCCTGGGGGGCTTATATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGAAAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTGATGaa                                                                                                          >  2:1004278/1‑150 (MQ=255)
                                             ctctTTTTTGAATATCCATCCCTGGGGGGATTTTATCGTCTTTGCTTTACCGCAAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGAAATGGTGTTTAATCTACATTGAGGACTGATTATTATGAAtt                                                                                                        >  1:1312356/1‑150 (MQ=255)
                                               ctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTaa                                                                                                      >  2:485539/1‑149 (MQ=255)
                                               ctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTAc                                                                                                      >  2:1511655/1‑150 (MQ=255)
                                               ctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTAc                                                                                                      <  2:212855/150‑1 (MQ=255)
                                               ctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGCACTGATGGTTATGAATTAc                                                                                                      >  1:1486969/1‑150 (MQ=255)
                                               ctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGGCATTGAGGACTGATGGTTATGAATTAc                                                                                                      >  2:973161/1‑150 (MQ=255)
                                                ttttttGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACt                                                                                                     <  2:101102/150‑1 (MQ=255)
                                                  ttttGAATATCCATCCCTGGGGGGCTTTTATCGGCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACttt                                                                                                   <  2:1411020/150‑1 (MQ=255)
                                                      gAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTTTTAATTACTTTTccg                                                                                               >  1:966235/1‑148 (MQ=255)
                                                       aaTATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCCGc                                                                                              >  1:1033218/1‑150 (MQ=255)
                                                            ccatccCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTTATGGTTATGAATTAATTTTAAGCTGGGc                                                                                         >  2:195419/1‑150 (MQ=255)
                                                                cccTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCCTTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGttt                                                                                     >  1:1814994/1‑150 (MQ=255)
                                                                  cTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTc                                                                                   >  1:122691/1‑150 (MQ=255)
                                                                   tGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCt                                                                                  <  2:674976/150‑1 (MQ=255)
                                                                        ggCGTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCgg                                                                             <  2:213495/150‑1 (MQ=255)
                                                                         gCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGt                                                                            >  1:1410268/1‑150 (MQ=255)
                                                                          cTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAAACGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGtt                                                                           >  2:1385884/1‑150 (MQ=255)
                                                                              tATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGt                                                                       <  1:1947301/150‑1 (MQ=255)
                                                                                   tCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTAttt                                                                  <  2:1930245/150‑1 (MQ=255)
                                                                                    cTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTta                                                                 >  1:1085823/1‑150 (MQ=255)
                                                                                    cTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTta                                                                 >  1:533601/1‑150 (MQ=255)
                                                                                        gCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCAATGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGGTTTCTGCCGGTTGGTTTATGtactac                                                             >  1:1777458/1‑150 (MQ=255)
                                                                                         cTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACg                                                            >  2:1625771/1‑150 (MQ=255)
                                                                                         cTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGAGGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACg                                                            >  2:1387027/1‑150 (MQ=255)
                                                                                          tttACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCCGATGGGTTATTAACTACGc                                                           >  2:516484/1‑150 (MQ=255)
                                                                                            tACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGGTATTTACTACGCTg                                                         <  2:708774/150‑1 (MQ=255)
                                                                                              ccGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGa                                                       <  2:516488/150‑1 (MQ=255)
                                                                                               cGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGAt                                                      <  1:485539/150‑1 (MQ=255)
                                                                                                gCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGAtt                                                     >  2:953196/1‑150 (MQ=255)
                                                                                                  cAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTg                                                   >  2:1757588/1‑150 (MQ=255)
                                                                                                   aGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGt                                                  >  2:509540/1‑150 (MQ=255)
                                                                                                   aGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGt                                                  >  2:1827899/1‑150 (MQ=255)
                                                                                                     ggCGTCGGCCTCAAGTATAAAGCAGACAAAAACAAAAACACCCTTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGTATTACTTTGCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCa                                                <  1:1384262/150‑1 (MQ=255)
                                                                                                     ggCGTCGGCCGCAAGTCTAAAGCAGATCAAACCAACAACACCCTTGCGCAGGCAATGGTGTTTAATCGGCATTGAGGACTGATGGGTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCa                                                <  1:266155/150‑1 (MQ=255)
                                                                                                           ggCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTATTCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTTTAACCTTCg                                          >  2:1039639/1‑150 (MQ=255)
                                                                                                            gCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCgg                                         >  1:1543284/1‑150 (MQ=255)
                                                                                                            gCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCgg                                         >  1:655047/1‑150 (MQ=255)
                                                                                                            gCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCgg                                         <  1:1385884/150‑1 (MQ=255)
                                                                                                             ccTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGt                                        <  2:1620816/150‑1 (MQ=255)
                                                                                                             ccTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTAAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGt                                        >  2:1341981/1‑150 (MQ=255)
                                                                                                              cTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTa                                       >  1:456719/1‑150 (MQ=255)
                                                                                                              cTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTa                                       >  2:1553236/1‑150 (MQ=255)
                                                                                                                  aGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTa                                   <  1:413161/150‑1 (MQ=255)
                                                                                                                  aGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTa                                   <  2:1999273/150‑1 (MQ=255)
                                                                                                                       aaaGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTg                              <  2:73206/150‑1 (MQ=255)
                                                                                                                       aaaGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGAGGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTg                              <  1:390700/150‑1 (MQ=255)
                                                                                                                       aaaGCAGATAAAAACAAAAACACCATTGAGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATCGTTATGAATTACTTATCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTg                              >  1:245398/1‑150 (MQ=255)
                                                                                                                        aaGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTgc                             <  1:1631323/150‑1 (MQ=255)
                                                                                                                         accAGATAAATAAAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGCCTGCTGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGGTATATACTACGCTGGATTTGTCACCTTCGGTAACTATTTTgcg                            <  2:22627/148‑1 (MQ=255)
                                                                                                                          tcAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGTATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTgcgc                           <  1:763279/149‑1 (MQ=255)
                                                                                                                          gCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTgcgc                           >  2:1454737/1‑150 (MQ=255)
                                                                                                                          gCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTgcgc                           >  1:1184841/1‑150 (MQ=255)
                                                                                                                           cAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCt                          <  1:1596681/150‑1 (MQ=255)
                                                                                                                             gATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCCTTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTgg                        <  2:1230298/150‑1 (MQ=255)
                                                                                                                             gATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGATGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCTGTTGGGTTATTTACTACGATGGATTTGTCACCTTCGGTAACTATTTTGCGCTgt                        >  1:1428201/1‑149 (MQ=255)
                                                                                                                             gATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTACGGTAACTATTTTGCGCTgg                        >  1:1806662/1‑150 (MQ=255)
                                                                                                                               tAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGtt                      <  2:1697790/150‑1 (MQ=255)
                                                                                                                                aaaaacaaaaacACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTa                     >  1:1431041/1‑150 (MQ=255)
                                                                                                                                 aaaacaaaaacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGGTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAg                    <  2:1810100/150‑1 (MQ=255)
                                                                                                                                   cgattaaacaCCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGAGCTGGTTAGaa                  >  2:1653821/6‑150 (MQ=255)
                                                                                                                                   aacaaaaacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCTGTAACTATTTTGCGCTGGTTAGaa                  >  2:1243209/1‑150 (MQ=255)
                                                                                                                                       aaaacaCCTTTGCGCAGGCAATGGTGTTTAATCGTCTTTGAGGACTGATTGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAAttt              <  1:110449/150‑1 (MQ=255)
                                                                                                                                       aaaacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAAttt              <  2:1481536/150‑1 (MQ=255)
                                                                                                                                         aacaCCATTGTGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTa            >  2:1925986/1‑150 (MQ=255)
                                                                                                                                         aacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTa            <  1:238754/150‑1 (MQ=255)
                                                                                                                                         aacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTAATGAATTACTTTTCATCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAAGTATTTTGCGCTGGTTAGAAATTTTa            >  1:1204339/1‑150 (MQ=255)
                                                                                                                                         aacaCCATTGCCTGCGCAATGGTGTTTAATCGTCATAGAGGACTGAAGGTTATGAATTACTTTACAGCGGGGCGTTTTCAGCCGGGTGGGTAATTAAATACGCTGGAATTGTCACCTTCGGTAACTATTTTGCGCTGGTAAGAAATTTTa            >  2:1046522/1‑150 (MQ=255)
                                                                                                                                           cacCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATg          <  2:1106601/150‑1 (MQ=255)
                                                                                                                                           cacCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATg          <  1:1060070/150‑1 (MQ=255)
                                                                                                                                           cacCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATg          >  2:890267/1‑150 (MQ=255)
                                                                                                                                           cacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATg          >  1:2011082/1‑150 (MQ=255)
                                                                                                                                            acCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATgg         <  2:1644639/150‑1 (MQ=255)
                                                                                                                                              cATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtc       >  2:1552977/1‑150 (MQ=255)
                                                                                                                                              cATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGCTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtc       <  1:1039639/150‑1 (MQ=255)
                                                                                                                                                 tGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccag    <  2:1951954/150‑1 (MQ=255)
                                                                                                                                                 tGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccag    >  1:247623/1‑150 (MQ=255)
                                                                                                                                                  gcgcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccagt   >  2:629104/1‑150 (MQ=255)
                                                                                                                                                   cgcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTATAAATTTTATGgtccagtc  >  1:1315949/1‑150 (MQ=255)
                                                                                                                                                   cgcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccagtc  >  1:744568/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCAGTC  >  NC_000913/1555633‑1555929

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: