Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,381,849 A→C 5.5% E369D (GAA→GAC) ‡ degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,381,8490AC 5.5% 333.6 / 8.1 109E369D (GAA→GAC) ‡degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (57/46);  new base C (4/2);  total (61/48)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 111 total aligned reads displayed.

CTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCG  >  NC_000913/3381706‑3381998
                                                                                                                                               |                                                                                                                                                     
cTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATc                                                                                                                                                 <  2:1694975/150‑1 (MQ=255)
 tGAATAGCTTTGCTGAGTTGCGATCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCa                                                                                                                                                <  2:1196803/150‑1 (MQ=255)
 ggAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCa                                                                                                                                                <  2:332178/149‑1 (MQ=255)
    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)
                              aTCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGc                                                                                                                   <  2:1298847/150‑1 (MQ=255)
                                 gcgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCCGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGAt                                                                                                                >  1:496778/1‑150 (MQ=255)
                                 gcgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGAt                                                                                                                >  2:1594872/1‑150 (MQ=255)
                                  cgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACTGCAAACCACTGGAAGTAGAAGTGACGCACGATACCAGCACATCTTCGTCGGCCAGCGCTGAAATGATAACCCCTGCGCTGGAAGGTGCAACGTTGAGccaag                                                                                                               >  2:1286415/1‑146 (MQ=255)
                                  cgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATg                                                                                                               >  2:1366872/1‑150 (MQ=255)
                                    accaccGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                                                                             >  1:893563/1‑150 (MQ=255)
                                     ccaccGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                                                                            >  2:636268/1‑150 (MQ=255)
                                      caccGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCa                                                                                                           >  2:880097/1‑150 (MQ=255)
                                       accGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGATGTAGAAGTGACGCTCGATATCAGCCCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                                                                          >  1:918569/1‑150 (MQ=255)
                                       accGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                                                                          >  1:1013554/1‑150 (MQ=255)
                                       accGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                                                                          <  2:895217/150‑1 (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)
                                                     gAAAGTGATGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATggcgg                                                                                            <  1:849581/150‑1 (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:67508/1‑150 (MQ=255)
                                                                                cGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGt                                                                 >  2:1232272/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)
                                                                                       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)
                                                                                                  aGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctgc                                               >  2:1048326/1‑150 (MQ=255)
                                                                                                   gAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctgct                                              >  2:158372/1‑150 (MQ=255)
                                                                                                    aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                             <  2:505633/150‑1 (MQ=255)
                                                                                                    aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                             <  2:811640/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:636268/150‑1 (MQ=255)
                                                                                                            cTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctggc                                     <  1:616816/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)
                                                                                                                     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)
                                                                                                                                 tcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATc                >  2:1118200/1‑150 (MQ=255)
                                                                                                                                  cttccAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCa               >  1:46387/4‑150 (MQ=255)
                                                                                                                                  cttccAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCa               >  2:4555/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)
                                                                                                                                               |                                                                                                                                                     
CTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCG  >  NC_000913/3381706‑3381998

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: