Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,704 G→C 5.9% intergenic (+23/‑205) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7040GC 5.9% 467.7 / 16.0 138intergenic (+23/‑205)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (72/57);  new base C (4/4);  total (77/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 141 total aligned reads displayed.

GACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGA  >  NC_000913/3334558‑3334852
                                                                                                                                                  |                                                                                                                                                     
      tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCgg                                                                                                                                               >  1:222549/1‑148 (MQ=255)
         gCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATg                                                                                                                                           <  3:601139/149‑1 (MQ=255)
           gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATgg                                                                                                                                          <  1:273131/148‑1 (MQ=255)
            tGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGAt                                                                                                                                            <  3:757070/145‑1 (MQ=255)
            tGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGtt                                                                                                                                        >  4:101930/1‑149 (MQ=255)
             gCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGtt                                                                                                                                        <  2:63677/148‑1 (MQ=255)
             gCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGC‑GCTGCGCGGGATgaggg                                                                                                                                       <  4:513791/149‑5 (MQ=255)
              ccGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                       <  1:781490/148‑1 (MQ=255)
               cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCa                                                                                                                                     <  4:865292/149‑1 (MQ=255)
                  ggAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAATGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAAC‑GCGCGGGATGGTTCCAGt                                                                                                                                   >  1:684565/1‑148 (MQ=255)
                    aagaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAg                                                                                                                                >  3:651509/1‑149 (MQ=255)
                     agaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAg                                                                                                                                >  2:816497/1‑148 (MQ=255)
                       aagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGtt                                                                                                                              <  1:170066/148‑1 (MQ=255)
                       aagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGtt                                                                                                                              <  2:899050/148‑1 (MQ=255)
                         gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTGGCGCACATCGCTGTTCAACGCGATCGTTAAGCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTATGTTCCa                                                                                                                           <  4:719740/149‑1 (MQ=255)
                         gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTcc                                                                                                                            <  2:676962/148‑1 (MQ=255)
                         gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTcc                                                                                                                            <  1:183603/148‑1 (MQ=255)
                          cAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCAt                                                                                                                          <  3:895839/149‑1 (MQ=255)
                              cAAAGCCATCGCAGCGTTCCAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACCTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGTGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                      <  3:839716/149‑1 (MQ=255)
                               aaaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGATCCATaaaaa                                                                                                                     >  4:50066/1‑149 (MQ=255)
                                 aGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAac                                                                                                                    >  2:213869/1‑148 (MQ=255)
                                    cATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACtt                                                                                                                <  4:540721/149‑1 (MQ=255)
                                             gTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                       <  4:912642/149‑1 (MQ=255)
                                               tACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctct                                                                                                      >  1:821683/1‑148 (MQ=255)
                                                  aGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                  >  4:629312/1‑149 (MQ=255)
                                                       cTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACAcc                                                                                              <  1:213869/148‑1 (MQ=255)
                                                       cTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCt                                                                                             >  3:912456/1‑149 (MQ=255)
                                                        tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                            >  4:243934/1‑149 (MQ=255)
                                                          ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGa                                                                                           >  2:36729/1‑148 (MQ=255)
                                                            ggACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaaaa                                                                                        >  3:127899/1‑149 (MQ=255)
                                                               cacCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAAt                                                                                      <  2:610870/148‑1 (MQ=255)
                                                               cacCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                     <  3:926622/149‑1 (MQ=255)
                                                               cacCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                     <  4:184981/149‑1 (MQ=255)
                                                                acCCCTTGGCGAGAAGCACTCATTGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAATACACTTATTCAGCTCTATAACACCTGAAAAATat                                                                                    <  4:664587/149‑2 (MQ=255)
                                                                 ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAgg                                                                                   <  3:50066/149‑1 (MQ=255)
                                                                   ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAggg                                                                                  >  1:904758/1‑148 (MQ=255)
                                                                   ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                 >  3:678585/1‑149 (MQ=255)
                                                                    cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                 >  2:203131/1‑148 (MQ=255)
                                                                    cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                 >  2:54862/1‑148 (MQ=255)
                                                                    cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                >  3:725559/1‑149 (MQ=255)
                                                                     ttGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                     >  4:342319/1‑143 (MQ=255)
                                                                     ttGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                     <  3:342319/143‑1 (MQ=255)
                                                                        gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                    <  4:401965/141‑1 (MQ=255)
                                                                        gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                    >  3:401965/1‑141 (MQ=255)
                                                                          gaaaaGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCGATAACACCTGAAAAATAGGGATGTGTAt                                                                          >  3:843444/4‑149 (MQ=255)
                                                                               gCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGc‑gcgcg                                                                                                                                                >  1:164559/1‑74 (MQ=255)
                                                                               gCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGc‑gcgcg                                                                                                                                                <  2:164559/74‑1 (MQ=255)
                                                                                   tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                                 >  4:20076/1‑149 (MQ=255)
                                                                                    cATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCCGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                                 <  2:513272/148‑1 (MQ=255)
                                                                                     aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                               >  4:146756/1‑149 (MQ=255)
                                                                                     aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                               >  3:485623/1‑149 (MQ=255)
                                                                                      tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                               <  1:280820/148‑1 (MQ=255)
                                                                                      tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                               <  1:787844/148‑1 (MQ=255)
                                                                                       cGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                              >  2:37838/1‑148 (MQ=255)
                                                                                              gcgcACATCGCTGTTCAACGCGATCGTTAGTCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATAtt                                                      >  3:861600/1‑149 (MQ=255)
                                                                                              gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATAtt                                                      <  3:25120/149‑1 (MQ=255)
                                                                                                   caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGTGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcaca                                                  >  2:831036/1‑148 (MQ=255)
                                                                                                       gCTGTTCAACGCGATCGTTAGTCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCtt                                             <  4:861600/149‑1 (MQ=255)
                                                                                                           ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTATTTACa                                         <  3:496549/149‑1 (MQ=255)
                                                                                                            tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACa                                         <  2:839026/148‑1 (MQ=255)
                                                                                                               aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGaa                                     >  3:379997/1‑149 (MQ=255)
                                                                                                               aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGa                                      >  1:384044/1‑148 (MQ=255)
                                                                                                                   gATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGAAACttt                                 >  3:6220/1‑149 (MQ=255)
                                                                                                                    aTCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACtttt                                >  4:219550/1‑149 (MQ=255)
                                                                                                                      cGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                              <  4:651509/149‑1 (MQ=255)
                                                                                                                      cGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGCTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                               <  2:904758/148‑1 (MQ=255)
                                                                                                                        ttAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                             <  1:54862/148‑1 (MQ=255)
                                                                                                                          aaTCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGc                          >  4:782058/1‑149 (MQ=255)
                                                                                                                              ccctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                      >  4:731018/1‑149 (MQ=255)
                                                                                                                              ccctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                      >  4:732561/1‑149 (MQ=255)
                                                                                                                               cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                      >  1:87804/1‑148 (MQ=255)
                                                                                                                               cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                     >  3:203310/1‑149 (MQ=255)
                                                                                                                                ctcccctcCTCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACAcc                                                                                              >  1:832810/1‑75 (MQ=255)
                                                                                                                                ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACAcc                                                                                              <  2:832810/75‑1 (MQ=255)
                                                                                                                                 tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg                    >  1:848712/1‑148 (MQ=255)
                                                                                                                                    cctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCAt                >  3:773492/1‑149 (MQ=255)
                                                                                                                                     ctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATACACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATc               <  3:67457/149‑1 (MQ=255)
                                                                                                                                      tcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa              >  3:179029/1‑149 (MQ=255)
                                                                                                                                       cATCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg             >  4:790892/1‑149 (MQ=255)
                                                                                                                                       cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTc                                                     >  4:186840/1‑109 (MQ=255)
                                                                                                                                       cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa              <  1:36729/148‑1 (MQ=255)
                                                                                                                                       cATCCCGCGC‑GCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg             >  3:271739/1‑149 (MQ=255)
                                                                                                                                        aTCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTatatat                                                                >  3:76858/1‑97 (MQ=255)
                                                                                                                                        aTCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTatatat                                                                <  4:76858/97‑1 (MQ=255)
                                                                                                                                        aTCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg             >  1:403958/1‑148 (MQ=255)
                                                                                                                                        aTCCCGCGC‑GCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg            >  3:131779/1‑149 (MQ=255)
                                                                                                                                         tCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                    >  1:321862/1‑76 (MQ=39)
                                                                                                                                         tCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                    <  2:321862/76‑1 (MQ=39)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                       <  2:614439/57‑1 (MQ=255)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                       >  1:614439/1‑57 (MQ=255)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTTACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa           >  3:198604/1‑149 (MQ=255)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg            <  1:908458/148‑1 (MQ=255)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg            <  1:267750/148‑1 (MQ=255)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa           <  4:6220/149‑1 (MQ=255)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa           >  3:118387/1‑149 (MQ=255)
                                                                                                                                         tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa           <  4:912456/149‑1 (MQ=255)
                                                                                                                                           ccGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGt         >  4:140599/1‑149 (MQ=255)
                                                                                                                                                 aGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGt    <  1:816497/148‑1 (MQ=255)
                                                                                                                                                 aGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTg   <  3:170335/149‑1 (MQ=255)
                                                                                                                                                  gc‑gcgcGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGa  <  4:118387/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGA  >  NC_000913/3334558‑3334852

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑