Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,381,859 C→T 5.5% P373S (CCA→TCA)  degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,381,8590CT 5.5% 350.8 / 5.5 111P373S (CCA→TCA) degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (59/45);  new base T (4/2);  total (63/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

Only 100 of 113 total aligned reads displayed.

ATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTG  >  NC_000913/3381710‑3382008
                                                                                                                                                     |                                                                                                                                                     
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)
                             gcgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGAt                                                                                                                          >  2:1594872/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                                                                                                                    >  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)
                                          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)
                                                                        gTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAg                                                                               <  1:1341136/150‑1 (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)
                                                                                 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)
                                                                                              aGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctgc                                                         >  2:1048326/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)
                                                                                                 aGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGAAACGATGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAACGCGCTTCAGTTGTCCAAGGAACCCCAGCTGCTCa                                                      >  1:272809/1‑150 (MQ=255)
                                                                                                    gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                   >  2:1330367/1‑150 (MQ=255)
                                                                                                    gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                   >  1:673661/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)
                                                                                                               ccAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaa                                        >  2:1565182/1‑150 (MQ=255)
                                                                                                               ccAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGAGCAACGTTGAGCGATGGGCAGCGAAAAGATGGAGGCAAAGGTTTTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGaaaa                                        >  1:878846/1‑150 (MQ=255)
                                                                                                                 aGCACCTCTTCGTCGGCCAGCGCGGAAATGCTCACGCCAGCGCTGGAAGGGGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaaa                                      <  1:1064936/150‑1 (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)
                                                                                                                                tccAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                       >  2:533814/2‑150 (MQ=255)
                                                                                                                                   aGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGc                    <  2:46387/150‑1 (MQ=255)
                                                                                                                                    gcgcTGAAATGATCACGCCAGCGCTGGAAGGGGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCg                   <  1:1103970/150‑1 (MQ=255)
                                                                                                                                       cTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                <  1:194727/150‑1 (MQ=255)
                                                                                                                                       cTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                >  2:651776/1‑150 (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:1018350/1‑150 (MQ=255)
                                                                                                                                           aaTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcg            >  2:448259/1‑150 (MQ=255)
                                                                                                                                               aTCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGAt        >  1:810985/1‑150 (MQ=255)
                                                                                                                                                tCACGCCAGCGCTGGAAGGGGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATc       <  2:223461/150‑1 (MQ=255)
                                                                                                                                                 cACGCCATCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg      <  1:835152/150‑1 (MQ=255)
                                                                                                                                                    gCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGt   <  2:661500/150‑1 (MQ=255)
                                                                                                                                                     ccAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTg  >  2:329919/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTG  >  NC_000913/3381710‑3382008

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