Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,769,180 C→G A44A (GCG→GCC)  pscL ← SctL family type III secretion system stator protein PscL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,769,1800CG99.8% 2026.8 / NA 515A44A (GCG→GCC) pscLSctL family type III secretion system stator protein PscL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (266/248);  minor base T (1/0);  total (267/248)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 534 total aligned reads displayed.

GCGATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCCGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTGT  >  NC_008463/3769033‑3769326
                                                                                                                                                   |                                                                                                                                                  
 cGATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCggcg                                                                                                                                                 >  2:2799548/1‑150 (MQ=255)
  gATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTACTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTcggcgg                                                                                                                                                >  2:5391185/1‑150 (MQ=255)
  gATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGGAGCGTTTCCTGGATCAGCCCGGCCAGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGAAAACCTCGTGGGCCTCGCGCTCGATTTcggcgg                                                                                                                                                <  1:9676193/150‑1 (MQ=255)
   aTCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCggcggc                                                                                                                                               >  1:10388444/1‑150 (MQ=255)
    tCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCa                                                                                                                                              <  2:120002/150‑1 (MQ=255)
       aCCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGc                                                                                                                                           >  1:7077829/1‑150 (MQ=255)
         cTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTc                                                                                                                                         >  1:10470247/1‑150 (MQ=255)
          tGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTcg                                                                                                                                        <  1:5001130/150‑1 (MQ=255)
            ccGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTcgcg                                                                                                                                      >  2:5008689/1‑150 (MQ=255)
             cGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGt                                                                                                                                     <  2:7646768/150‑1 (MQ=255)
              ggTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTg                                                                                                                                    >  1:1135967/1‑150 (MQ=255)
               gtagtaGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGc                                                                                                                                   >  1:6259713/1‑150 (MQ=255)
                tagtagCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGcc                                                                                                                                  >  1:8715639/1‑150 (MQ=255)
                   tagCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCt                                                                                                                               >  2:8485660/1‑150 (MQ=255)
                       gATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCGTcca                                                                                                                           >  1:1555637/1‑150 (MQ=255)
                        aTTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTccac                                                                                                                          >  1:9138053/1‑150 (MQ=255)
                               cgcAGTAGCGTGTCCTGGATCAGCCCGGCCTGTCGGAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCTGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGATCGATTTCGGCGGCACGCTCGCGTGCCGCTTCAACCAGTCgg                                                                                                                   <  2:7945126/150‑1 (MQ=255)
                                gcAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGt                                                                                                                  <  2:5904264/150‑1 (MQ=255)
                                 cAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCCGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGtt                                                                                                                 <  2:1407925/150‑1 (MQ=255)
                                  aGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTg                                                                                                                <  2:2630823/150‑1 (MQ=255)
                                   gTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTgg                                                                                                               <  2:8919610/150‑1 (MQ=255)
                                   gTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTgg                                                                                                               >  1:7381738/1‑150 (MQ=255)
                                   gTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTgg                                                                                                               <  1:558229/150‑1 (MQ=255)
                                    tAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGgc                                                                                                              >  2:3789923/1‑150 (MQ=255)
                                       cGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGgcgcg                                                                                                           >  1:4797887/1‑149 (MQ=255)
                                        gTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGgcgcgg                                                                                                          >  1:5890998/1‑148 (MQ=255)
                                        gTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTg                                                                                                          <  1:9036300/150‑1 (MQ=255)
                                         tttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGa                                                                                                         >  2:9731668/1‑150 (MQ=255)
                                         tttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGa                                                                                                         <  2:5314556/150‑1 (MQ=255)
                                         tttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGa                                                                                                         >  2:8048006/1‑150 (MQ=255)
                                         tttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGACCGGCCTCCCAGCCCCGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATGTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGa                                                                                                         <  1:2933551/150‑1 (MQ=255)
                                         tgtCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGAGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCTCTCGATTTCGGCGGCACGCTCGCGTGCCGATTCCTCCAGGCGGTTGGCGCTGa                                                                                                         <  2:6650156/148‑1 (MQ=255)
                                          ttCCTGGTTCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGAAAACCGCGAGGGCCTCGCGCGCGAGTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAg                                                                                                        >  2:5753799/1‑150 (MQ=255)
                                          ttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAg                                                                                                        >  1:3846852/1‑150 (MQ=255)
                                           tCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGTCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAgg                                                                                                       <  2:426942/150‑1 (MQ=255)
                                           tCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCAAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCACGATTTCGGCGGCACGCACGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAgt                                                                                                       >  2:3554402/1‑149 (MQ=255)
                                                 acccccccGGCCTGGCGCAGCCCCGCCTCCCCCAGACCGCCCTCCCCGCCCAGGCGCTTCTGCTCCTGGTAACCCTCCTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTcc                                                                                                 <  1:9218798/145‑1 (MQ=255)
                                                   cAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGGGCGGCTCCTCCTCGTAAAGCTCGTGGGCCTCGCGCGCGATTTCTGCCGCTCGCTAGCGTGCCGCTGACTCCAGGCGGTTGGCGCGGAGGTAGACCTg                                                                                               >  1:8499693/1‑150 (MQ=255)
                                                   cAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTg                                                                                               >  1:10147280/1‑150 (MQ=255)
                                                   cAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGGTGGCGCTGAGGTAGTCCTg                                                                                               >  1:2100480/1‑150 (MQ=255)
                                                     gCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGt                                                                                             <  2:3899603/150‑1 (MQ=255)
                                                     gCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGt                                                                                             <  2:9588636/150‑1 (MQ=255)
                                                      cccGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTa                                                                                            >  1:8227551/1‑150 (MQ=255)
                                                        cGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTTGTAGt                                                                                          >  1:8949834/1‑150 (MQ=255)
                                                         ggCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTc                                                                                         >  2:2344057/1‑150 (MQ=255)
                                                          gCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTcg                                                                                        <  1:8061749/150‑1 (MQ=255)
                                                            cTGGCGCAGACGCGCCTCCTCCCGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTcgcg                                                                                      <  2:10617766/150‑1 (MQ=255)
                                                             tGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCgg                                                                                     >  2:10215548/1‑150 (MQ=255)
                                                              ggCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGGTGGCGCTGAGGTAGTCCTGGTAGTCGCggg                                                                                    >  1:8776943/1‑150 (MQ=255)
                                                               gcgcAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGgc                                                                                   >  1:7713506/1‑150 (MQ=255)
                                                               gcgcAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGgc                                                                                   <  2:5550697/150‑1 (MQ=255)
                                                                 gcAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGgcgc                                                                                 <  2:1686959/150‑1 (MQ=255)
                                                                             cctccAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGCTTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCt                                                                     >  1:7434591/1‑150 (MQ=255)
                                                                              ctccAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTg                                                                    >  1:4132611/1‑150 (MQ=255)
                                                                                ccAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGcc                                                                  >  1:9979994/1‑150 (MQ=255)
                                                                                 cAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGccc                                                                 >  1:2505983/1‑150 (MQ=255)
                                                                                 cAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGccc                                                                 <  1:3114255/150‑1 (MQ=255)
                                                                                  aGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCg                                                                >  2:8228083/1‑150 (MQ=255)
                                                                                   gTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCgg                                                               <  1:6014867/150‑1 (MQ=255)
                                                                                   gTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCGGGTAGTCGCGGGCGCGCAACACGGCCGGCCCgg                                                               >  1:1005473/1‑150 (MQ=255)
                                                                                        gCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGcc                                                          <  2:2084838/150‑1 (MQ=255)
                                                                                            cccagcccagGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCAAGCGTGACGCTTCCACCAGGCGGTTGGAGATGTGGAAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGcc                                                      >  2:8365700/1‑150 (MQ=255)
                                                                                             ccagcccagGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCg                                                     >  2:2247997/1‑150 (MQ=255)
                                                                                                     ggCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGAc                                             <  1:2688437/150‑1 (MQ=255)
                                                                                                     ggCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGAc                                             >  1:6939121/1‑150 (MQ=255)
                                                                                                     ggCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGAc                                             <  1:9509238/150‑1 (MQ=255)
                                                                                                     ggCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGAc                                             <  2:10448000/150‑1 (MQ=255)
                                                                                                     ggCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGAc                                             >  2:4841976/1‑150 (MQ=255)
                                                                                                        gcTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTgg                                          <  2:2983523/150‑1 (MQ=255)
                                                                                                         cTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGc                                         >  2:9334338/1‑150 (MQ=255)
                                                                                                          ttCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCg                                        >  1:3648633/1‑150 (MQ=255)
                                                                                                          ttCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCGGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCg                                        >  1:10627237/1‑150 (MQ=255)
                                                                                                           tCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGt                                       >  1:7848075/1‑150 (MQ=255)
                                                                                                             tGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCt                                     <  2:1833387/150‑1 (MQ=255)
                                                                                                               cTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAg                                   >  1:6762407/1‑150 (MQ=255)
                                                                                                                tCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGt                                  >  1:5182486/1‑150 (MQ=255)
                                                                                                                 ccTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGtt                                 <  2:1596127/150‑1 (MQ=255)
                                                                                                                  cTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTc                                >  1:6609644/1‑150 (MQ=255)
                                                                                                                  cTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTc                                <  2:5285785/150‑1 (MQ=255)
                                                                                                                   tGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTca                               >  2:6255834/1‑150 (MQ=255)
                                                                                                                    ggTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTcaa                              <  2:6472704/150‑1 (MQ=255)
                                                                                                                    ggTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTcaa                              <  1:535997/150‑1 (MQ=255)
                                                                                                                       aacccTCGTGGGCCTCGCGCTCGAGTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTcaacaa                           <  1:7757689/147‑1 (MQ=255)
                                                                                                                       aaaCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTcaacaa                           <  1:5089617/150‑1 (MQ=255)
                                                                                                                         aCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAAt                         <  1:3812744/150‑1 (MQ=255)
                                                                                                                         aCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAAt                         <  2:7047323/150‑1 (MQ=255)
                                                                                                                           cTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATgg                       <  2:9128029/150‑1 (MQ=255)
                                                                                                                           cTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATgg                       <  1:6427926/150‑1 (MQ=255)
                                                                                                                             cGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGaa                     <  2:4052773/150‑1 (MQ=255)
                                                                                                                              gTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAg                    >  2:5703719/1‑150 (MQ=255)
                                                                                                                              gTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAg                    >  1:4363822/1‑150 (MQ=255)
                                                                                                                                gggCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCa                  <  1:2209656/150‑1 (MQ=255)
                                                                                                                                gggCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCa                  >  2:6774848/1‑150 (MQ=255)
                                                                                                                                  gCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCAtt                <  2:4407415/150‑1 (MQ=255)
                                                                                                                                   ccTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTc               <  2:8695673/150‑1 (MQ=255)
                                                                                                                                    cTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCa              >  1:10714253/1‑150 (MQ=255)
                                                                                                                                     ccgcgCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGATGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCAt             <  2:3323996/149‑1 (MQ=255)
                                                                                                                                          cTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCCGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCACCAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGACGCATTCAgcgctc        >  2:7701244/1‑150 (MQ=255)
                                                                                                                                              aTTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCtt    <  2:7375829/150‑1 (MQ=255)
                                                                                                                                               tttCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTg   >  1:1907767/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCCGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTGT  >  NC_008463/3769033‑3769326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑