Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,987,280 A→G R78R (CGT→CGC yecH ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,2800AG99.7% 997.6 / ‑4.2 307coding (234/240 nt)yecHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (155/151);  minor base . (1/0);  total (156/151)
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.44e-01

Only 100 of 313 total aligned reads displayed.

TGTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCGC  >  NC_000913/1987190‑1987369
                                                                                          |                                                                                         
tGTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCg                                                                                         >  2:3545506/1‑93 (MQ=255)
tGTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCg                                                                                         <  1:664713/93‑1 (MQ=255)
 gTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCgg                                                                                        >  1:5048031/1‑93 (MQ=255)
 gTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCgg                                                                                        >  2:7679725/1‑93 (MQ=255)
     tttATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCaaa                                                                                    <  1:9199022/93‑1 (MQ=255)
      ttATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAAt                                                                                   <  2:1366530/93‑1 (MQ=255)
       tATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATc                                                                                  >  1:9005155/1‑93 (MQ=255)
       tATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATc                                                                                  <  1:4439895/93‑1 (MQ=255)
       tATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATc                                                                                  >  1:4797158/1‑93 (MQ=255)
        aTCGACAAAGAACTCTCTCCCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCt                                                                                 <  1:3865162/93‑1 (MQ=255)
            aCAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTAct                                                                             >  1:6969575/1‑93 (MQ=255)
             cAAAAACCTCTCTCCCCCGGACACAAAAACCTAATGCCCCCGGGCGGCGGTTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTActc                                                                            <  2:7169635/93‑1 (MQ=255)
                 gAACTCTGTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGAt                                                                        >  2:6167460/1‑93 (MQ=255)
                  aaCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGTGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATc                                                                       <  1:5446005/93‑1 (MQ=255)
                  aaCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATc                                                                       <  1:4544329/93‑1 (MQ=255)
                   aCTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCg                                                                      >  1:3216543/1‑93 (MQ=255)
                     tctctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGt                                                                    >  1:6036248/1‑93 (MQ=255)
                     tctctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGt                                                                    >  2:401909/1‑93 (MQ=255)
                     tctctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGt                                                                    >  1:5874603/1‑93 (MQ=255)
                       tctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCg                                                                  <  1:3912064/93‑1 (MQ=255)
                        ctACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGa                                                                 >  1:2596455/1‑93 (MQ=255)
                         tACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaa                                                                >  1:4070825/1‑93 (MQ=255)
                          aCCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCCAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaca                                                               >  1:7648418/1‑91 (MQ=255)
                          aCCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaa                                                               >  1:5177795/1‑93 (MQ=255)
                            ccTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaaaa                                                             <  2:4538721/93‑1 (MQ=255)
                            ccTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGaaaaa                                                             >  1:3308597/1‑93 (MQ=255)
                             cTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAAc                                                            <  2:6088930/93‑1 (MQ=255)
                               cGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCg                                                          <  1:5926158/93‑1 (MQ=255)
                               cGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCg                                                          <  2:5094476/93‑1 (MQ=255)
                                 acacacAAAACTAATGCTCTCCTGCGGCAGGATTACTGCCGCCGCCTCGCATCAGCGGCGGCAAATCTTACTTTTATCGCTCGAAACACCgac                                                        >  2:207400/1‑91 (MQ=255)
                                   acacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCtt                                                      >  1:1013711/1‑93 (MQ=255)
                                   acacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCtt                                                      <  2:1210916/93‑1 (MQ=255)
                                   acacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCtt                                                      <  2:1796382/93‑1 (MQ=255)
                                    cacaAACCTAACGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCttt                                                     <  2:2515949/93‑1 (MQ=255)
                                    cacaAAACTAATGCTCTCGGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCttt                                                     <  1:3760084/93‑1 (MQ=255)
                                    cacaAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCttt                                                     <  2:7062287/93‑1 (MQ=255)
                                      caAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCg                                                   >  1:5730088/1‑93 (MQ=255)
                                       aaaaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGa                                                  >  1:7812424/1‑93 (MQ=255)
                                       aaaaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGa                                                  >  1:1468252/1‑93 (MQ=255)
                                        aaaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGaa                                                 <  1:3715779/93‑1 (MQ=255)
                                         aaCTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAg                                                >  2:824941/1‑93 (MQ=255)
                                            tAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGta                                             >  2:9156594/1‑93 (MQ=255)
                                            tAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGta                                             <  1:7382826/93‑1 (MQ=255)
                                             aaTGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGtat                                            >  1:7135892/1‑93 (MQ=255)
                                               tGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATaa                                          >  2:641896/1‑93 (MQ=255)
                                               tGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTATGATCGGTCGAAAAACCGTCTTTCGAAGGTATaa                                          >  2:1134898/1‑93 (MQ=255)
                                                gccctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATaaa                                         <  1:9135944/90‑1 (MQ=255)
                                                gCTCTGCAGCGGCAGGTTCACTGCCACCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATaaa                                         <  2:3447681/93‑1 (MQ=255)
                                                 ctctCCGGCGGCAGGTTTGCTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAt                                        <  1:8149035/93‑1 (MQ=255)
                                                  tctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAtt                                       <  1:2399055/93‑1 (MQ=255)
                                                  tctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAtt                                       <  1:6374121/93‑1 (MQ=255)
                                                   ctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAttt                                      >  2:9049812/1‑93 (MQ=255)
                                                   ctcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAAttt                                      >  2:9677440/1‑93 (MQ=255)
                                                    tcCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTg                                     >  1:3542198/1‑93 (MQ=255)
                                                     ccGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGc                                    <  2:6557985/93‑1 (MQ=255)
                                                     ccGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGc                                    <  2:5922709/93‑1 (MQ=255)
                                                     ccGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGc                                    >  2:4288549/1‑93 (MQ=255)
                                                       ggcggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGCAAAACCGTCTTTCGAAGGGATAAAATTTCCt                                  >  2:9111047/1‑93 (MQ=255)
                                                       ggcggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCt                                  <  1:9507587/93‑1 (MQ=255)
                                                         cggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCttt                                <  2:1073543/93‑1 (MQ=255)
                                                         cggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCttt                                <  1:2791187/93‑1 (MQ=255)
                                                         cggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCttt                                <  2:5527191/93‑1 (MQ=255)
                                                          ggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCtttt                               >  2:8706573/1‑93 (MQ=255)
                                                          ggcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCtttt                               <  2:9550207/93‑1 (MQ=255)
                                                           gcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtg                              >  1:557631/1‑93 (MQ=255)
                                                           gcAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtg                              >  2:4432414/1‑93 (MQ=255)
                                                            gagttttACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgc                             <  1:1886548/89‑1 (MQ=255)
                                                            cAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgc                             >  1:2113204/1‑93 (MQ=255)
                                                            cAGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTAAAAATTTGCCTTTtgc                             <  1:6012830/93‑1 (MQ=255)
                                                             aGGTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgct                            <  2:2805943/93‑1 (MQ=255)
                                                              ggTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctg                           >  2:8404259/1‑93 (MQ=255)
                                                              ggTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctg                           <  2:1680061/93‑1 (MQ=255)
                                                               gTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctgc                          >  2:4929952/1‑93 (MQ=255)
                                                               gTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctgc                          <  1:9765888/93‑1 (MQ=255)
                                                               gTTTACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTtgctgc                          >  1:6737932/1‑93 (MQ=255)
                                                                 ttACTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCa                        >  2:3970826/1‑93 (MQ=255)
                                                                   aCTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTGCGCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGa                      <  2:5242441/93‑1 (MQ=255)
                                                                   aCTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGa                      <  1:851569/93‑1 (MQ=255)
                                                                    cTGCCGCCGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaa                     >  2:1967797/1‑93 (MQ=255)
                                                                     tGCCGCGGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGaaa                    <  1:5458718/93‑1 (MQ=255)
                                                                      gccgccGTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAAc                   >  2:5672671/1‑93 (MQ=255)
                                                                         gccgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCt                >  1:8335251/1‑93 (MQ=255)
                                                                           cgTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTAc              <  2:5028297/93‑1 (MQ=255)
                                                                            gTTTCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTAcc             <  2:3157045/93‑1 (MQ=255)
                                                                             tttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCa            >  2:9380431/1‑93 (MQ=255)
                                                                             tttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCa            <  2:7018029/93‑1 (MQ=255)
                                                                             tttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCa            >  2:6537798/1‑93 (MQ=255)
                                                                             tttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTACCAGAAACGCTACCa            >  2:471844/1‑93 (MQ=255)
                                                                              ttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAg           <  2:2952019/93‑1 (MQ=255)
                                                                              ttCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAg           <  1:4746599/93‑1 (MQ=255)
                                                                               tCAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGc          <  2:6261995/93‑1 (MQ=255)
                                                                                cAAATCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGct         >  1:684967/1‑93 (MQ=255)
                                                                                  aaTCAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGctct       <  2:3778206/93‑1 (MQ=255)
                                                                                    acAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTcc     <  2:10263/92‑1 (MQ=255)
                                                                                     cAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTccc    <  1:1323159/93‑1 (MQ=255)
                                                                                     cAGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTccc    >  1:8840339/1‑93 (MQ=255)
                                                                                      aGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACTGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCg   >  2:3032533/1‑93 (MQ=255)
                                                                                      aGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCg   <  2:7131854/93‑1 (MQ=255)
                                                                                      aGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCg   >  2:6837242/1‑93 (MQ=255)
                                                                                      aGTGGCGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCg   >  2:5218555/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTATTTTATCGACAAAGAACTCTCTACCCTCGACACACAAAACTAATGCTCTCCGGCGGCAGGTTTACTGCCGCCGTTTCAAATCAGTGACGGCAAATCTTACTCTGATCGGTCGAAAAACCGTCTTTCGAAGGTATAAATTTGCCTTTTGCTGCCAGAAACGCTACCAGCTCTCCCGC  >  NC_000913/1987190‑1987369

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: