Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,381,857 C→T 5.3% T372M (ACG→ATG) ‡ degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,381,8570CT 5.3% 357.9 / 4.3 111T372M (ACG→ATG) ‡degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (60/45);  new base T (4/2);  total (64/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 111 total aligned reads displayed.

ATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCG  >  NC_000913/3381710‑3382004
                                                                                                                                                   |                                                                                                                                                   
aTAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCCCCTCTTCGTCGGCCAGCGCTGAAATGATCACGc                                                                                                                                                   >  2:379774/1‑150 (MQ=255)
  aGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGAGACGCTCGATACCAGCACCTCTGCGTCGGCCAGCGCGGCAATGATCACGCCa                                                                                                                                                 >  1:910834/1‑150 (MQ=255)
    cTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAgc                                                                                                                                               >  2:194727/1‑150 (MQ=255)
      ttGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAgcgc                                                                                                                                             <  2:1505953/150‑1 (MQ=255)
      ttGCTGAGTTGCGCTCTCGTACCGCGCCCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAgcgc                                                                                                                                             >  2:633049/1‑150 (MQ=255)
       tGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCt                                                                                                                                            <  2:1716937/150‑1 (MQ=255)
            aGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCTCTGGAAg                                                                                                                                       <  1:1476852/150‑1 (MQ=255)
             gTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGACAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGGCGTGATCATTTCAGCGCTggccga                                                                                                                                      <  1:766493/150‑5 (MQ=255)
             gTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAgg                                                                                                                                      >  1:661500/1‑150 (MQ=255)
                   ctctCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAAc                                                                                                                                <  1:1158992/150‑1 (MQ=255)
                     ctcGTATCGCGACCACCGAGCCGGGAACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCCCGCCCGCGCTGGAAGGTGCAACGt                                                                                                                              >  2:81190/1‑150 (MQ=255)
                       cGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTg                                                                                                                            <  2:546613/150‑1 (MQ=255)
                        gTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGa                                                                                                                           >  1:475356/1‑150 (MQ=255)
                         tATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAg                                                                                                                          <  2:1379449/150‑1 (MQ=255)
                          aTCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGc                                                                                                                         >  1:434590/1‑150 (MQ=255)
                          aTCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGc                                                                                                                         <  2:545042/150‑1 (MQ=255)
                             gcgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCCGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGAt                                                                                                                      >  1:496778/1‑150 (MQ=255)
                              cgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACTGCAAACCACTGGAAGTAGAAGTGACGCACGATACCAGCACATCTTCGTCGGCCAGCGCTGAAATGATAACCCCTGCGCTGGAAGGTGCAACGTTGAGccaag                                                                                                                     >  2:1286415/1‑146 (MQ=255)
                              cgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATg                                                                                                                     >  2:1366872/1‑150 (MQ=255)
                               gACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACAGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCCCGCCACCGCTGGAAGGTGCAACGGTGAGCGATgg                                                                                                                    >  2:1328606/1‑150 (MQ=255)
                                accaccGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                                                                                   >  1:893563/1‑150 (MQ=255)
                                 cccccGAGCCGGGCAGGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAATCCACTGGATGTAGAAGTCCCGCTCGATACCAGCCCCTCTTCGTCGGCCAGCGCAGAAATGATCCCGCCAGCGCTGGAAGGTGCAACGTTGCGCGAggggc                                                                                                                  >  2:253440/4‑148 (MQ=255)
                                 ccaccGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                                                                                  >  2:636268/1‑150 (MQ=255)
                                  caccGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCa                                                                                                                 >  2:880097/1‑150 (MQ=255)
                                   accGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGATGTAGAAGTGACGCTCGATATCAGCCCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                                                                                >  1:918569/1‑150 (MQ=255)
                                   accGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                                                                                <  2:895217/150‑1 (MQ=255)
                                   accGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                                                                                >  1:1013554/1‑150 (MQ=255)
                                    ccGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGc                                                                                                               <  2:641046/150‑1 (MQ=255)
                                     cGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCACACCACTGGAAGTAGAAGTGCCGCTCGATACCAGCCCCTCTTCGTCGGCCAGCGCTGAAATGATCCCGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCt                                                                                                              >  1:1145097/1‑150 (MQ=255)
                                     cGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCCCTTCGTCGGCCAGCGCTGAAATGAGCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCt                                                                                                              <  2:683456/150‑1 (MQ=255)
                                          cGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAg                                                                                                         <  2:572309/150‑1 (MQ=255)
                                           gggCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGa                                                                                                        >  2:1572087/1‑150 (MQ=255)
                                             gCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGATGGTGCAACGGTGAGCGATGGTCAGCTAAAAGATg                                                                                                      >  2:38172/1‑150 (MQ=255)
                                              cACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATgg                                                                                                     >  2:1212236/1‑150 (MQ=255)
                                               aCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTACAACGTTGAGCGATGGTCAGCTAAAAGATggc                                                                                                    >  2:924010/1‑150 (MQ=255)
                                                cGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCCCGCCAGCGCTGGAAGGTGCAACGGTGAGCGAGGGTCAGCTAAAAGAGggcg                                                                                                   >  1:819605/1‑150 (MQ=255)
                                                     ttGAAGCTTGGCCTGCTTCGTAACGGCAAAGCACTGGAAGTAGAAGTGACCCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaaa                                                                                              <  1:1477241/149‑1 (MQ=255)
                                                      tGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                                                                                             >  1:793752/1‑150 (MQ=255)
                                                        aaGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGCAGAAGTGCCGCTCGATACCAGCCCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGt                                                                                           <  1:341491/150‑1 (MQ=255)
                                                         aGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTGAGCTAAAAGATGGCGGCAAAGGTa                                                                                          >  2:840597/1‑150 (MQ=255)
                                                         aGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTa                                                                                          >  2:1688235/1‑150 (MQ=255)
                                                         aGCTAGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGACGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAggga                                                                                          >  1:87531/1‑148 (MQ=255)
                                                                ccTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAAt                                                                                   <  1:984866/150‑1 (MQ=255)
                                                                 ctgctgCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATc                                                                                  >  1:1655619/1‑150 (MQ=255)
                                                                  tgctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCg                                                                                 >  2:557629/1‑150 (MQ=255)
                                                                        gTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGACAAAGGTATTAAAATCGATGAAg                                                                           >  1:392100/1‑150 (MQ=255)
                                                                           aCGGCAAACCTCTTGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTg                                                                        <  2:277937/150‑1 (MQ=255)
                                                                            cGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGt                                                                       >  2:1232272/1‑150 (MQ=255)
                                                                            cGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGt                                                                       >  2:67508/1‑150 (MQ=255)
                                                                             ggCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGta                                                                      >  2:1617259/1‑149 (MQ=255)
                                                                             ggCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTc                                                                      >  2:1608567/1‑150 (MQ=255)
                                                                                aaaCCACTGGAAGTAGAAGTGACGCTCGATACAAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaaa                                                                   >  1:1600014/1‑150 (MQ=255)
                                                                                 aaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAg                                                                  >  2:555845/1‑150 (MQ=255)
                                                                                 aaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAg                                                                  <  2:1600014/150‑1 (MQ=255)
                                                                                   ccACTGGAAGTAGAAGTGACGCTCGAAACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCACTGTAAGGTGCAACGTAGAGCGATGGTGAGCTAAAAGATGCCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGa                                                                >  2:1172754/1‑150 (MQ=255)
                                                                                     aCTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAATGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAg                                                              <  1:67508/150‑1 (MQ=255)
                                                                                     aCTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAg                                                              <  2:1294144/150‑1 (MQ=255)
                                                                                     aCTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAg                                                              <  1:1232272/150‑1 (MQ=255)
                                                                                     aCTGGAAGTAGAAGTGACGCTCGATACCAGCACATCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAg                                                              >  2:1725129/1‑150 (MQ=255)
                                                                                      cTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGCAGc                                                             >  2:1566849/1‑150 (MQ=255)
                                                                                       tGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGcc                                                            <  2:630741/150‑1 (MQ=255)
                                                                                         gAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCa                                                          <  1:880097/150‑1 (MQ=255)
                                                                                           aGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgc                                                        <  2:1154993/150‑1 (MQ=255)
                                                                                              aGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctgc                                                     >  2:1048326/1‑150 (MQ=255)
                                                                                               gAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctgct                                                    >  2:158372/1‑150 (MQ=255)
                                                                                                aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                                   <  2:811640/150‑1 (MQ=255)
                                                                                                aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                                   <  2:505633/150‑1 (MQ=255)
                                                                                                aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTACTc                                                   >  1:1650668/1‑150 (MQ=255)
                                                                                                    gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                               >  2:1330367/1‑150 (MQ=255)
                                                                                                       gCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctgg                                            <  2:893563/150‑1 (MQ=255)
                                                                                                       gCTCGATACCAGCACCTCTTCGGCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGGGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctgg                                            <  2:117515/150‑1 (MQ=255)
                                                                                                        cTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctggc                                           <  1:616816/150‑1 (MQ=255)
                                                                                                        cTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctggc                                           <  1:636268/150‑1 (MQ=255)
                                                                                                         tCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCAAAAAGATCGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAATCCCAGCTGCTCAggctggct                                          >  1:701463/1‑150 (MQ=255)
                                                                                                           gATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTg                                        <  1:1291383/150‑1 (MQ=255)
                                                                                                            aTACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGACAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGc                                       <  2:392100/150‑1 (MQ=255)
                                                                                                               ccAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaa                                    >  2:1565182/1‑150 (MQ=255)
                                                                                                               ccAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGAGCAACGTTGAGCGATGGGCAGCGAAAAGATGGAGGCAAAGGTTTTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGaaaa                                    >  1:878846/1‑150 (MQ=255)
                                                                                                                 aGCACCTCTTCGTCGGCCAGCGCGGAAATGCTCACGCCAGCGCTGGAAGGGGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaaa                                  <  1:1064936/150‑1 (MQ=255)
                                                                                                                  gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTTCAACGTTGAGCGATGGTCAGCTAAAAGATGGCTGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAg                                 <  2:193584/150‑1 (MQ=255)
                                                                                                                     ccTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACg                              >  1:171282/1‑150 (MQ=255)
                                                                                                                      ctctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGGAACGTTGCGCGATCGTCAGCAAAAAGATGCCGGCAAACGAATTAAAATCGATGAAGTTTACAAAGGAAGCCCACCTGCTCAGCCTGGCTTCGAAAAAGCCGa                             >  2:1608155/1‑150 (MQ=255)
                                                                                                                        ctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGACCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtg                           <  1:1053212/150‑1 (MQ=255)
                                                                                                                         ttCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtgt                          <  1:158372/150‑1 (MQ=255)
                                                                                                                          tcgtcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtgtg                         >  1:1189140/1‑150 (MQ=255)
                                                                                                                          tcgtcgGCCAGCGCTGAAATGATCACGCCAGCGCGGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtgtg                         <  2:1279697/150‑1 (MQ=255)
                                                                                                                              cttccAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCa                     >  2:4555/4‑150 (MQ=255)
                                                                                                                              cttccAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCa                     >  1:46387/4‑150 (MQ=255)
                                                                                                                               ttccAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAt                    >  1:788271/3‑150 (MQ=255)
                                                                                                                                gCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                   <  1:1255581/150‑1 (MQ=255)
                                                                                                                                   aGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGc                <  2:46387/150‑1 (MQ=255)
                                                                                                                                    gcgcTGAAATGATCACGCCAGCGCTGGAAGGGGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCg               <  1:1103970/150‑1 (MQ=255)
                                                                                                                                       cTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa            <  1:194727/150‑1 (MQ=255)
                                                                                                                                         gAAATTCTCACGGAAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCACAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAAc          <  1:81190/150‑1 (MQ=255)
                                                                                                                                         gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAAc          >  1:194997/1‑150 (MQ=255)
                                                                                                                                         gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATAGATGAAGTTGTCAAAGGAAGCCCAGGTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGGCAAc          >  2:1146098/1‑150 (MQ=255)
                                                                                                                                           aaTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcg        >  2:448259/1‑150 (MQ=255)
                                                                                                                                           aaTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcg        >  2:1018350/1‑150 (MQ=255)
                                                                                                                                               aTCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGAt    >  1:810985/1‑150 (MQ=255)
                                                                                                                                                 cACGCCATCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg  <  1:835152/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCG  >  NC_000913/3381710‑3382004

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