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% 1787.1 / NA 456A44A (GCG→GCC) pscLSctL family type III secretion system stator protein PscL
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (221/234);  total (222/234)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 474 total aligned reads displayed.

GCGATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCCGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTGT  >  NC_008463/3769033‑3769326
                                                                                                                                                   |                                                                                                                                                  
 cGATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCggcg                                                                                                                                                 >  1:8707205/1‑150 (MQ=255)
  gATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTcggcgg                                                                                                                                                <  2:2026233/150‑1 (MQ=255)
   ctCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCggcggc                                                                                                                                               <  1:1876750/149‑1 (MQ=255)
    tCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCa                                                                                                                                              <  1:4276634/150‑1 (MQ=255)
    tCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCa                                                                                                                                              >  1:9507288/1‑150 (MQ=255)
      gACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACg                                                                                                                                            >  1:10244155/1‑150 (MQ=255)
      gACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACg                                                                                                                                            >  2:3432280/1‑150 (MQ=255)
        ccTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCt                                                                                                                                          <  1:4746513/150‑1 (MQ=255)
        ccTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCt                                                                                                                                          <  2:10340871/150‑1 (MQ=255)
           gCCGGTAGTCGCGATTGCAGCGCAGTAGCGAGTCCTGGATCAGCCCGGCCTGGCGCAGCCGCGCCTCCTCCAGGCCGGCCTCCCATCCCAGGCGCTTCTGCTCCTGGGAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTcgc                                                                                                                                       <  2:8712636/150‑1 (MQ=255)
           gCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTcgc                                                                                                                                       >  1:7475503/1‑150 (MQ=255)
           gCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTcgc                                                                                                                                       <  2:914379/150‑1 (MQ=255)
             cGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGt                                                                                                                                     >  1:212419/1‑150 (MQ=255)
                   tagCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCt                                                                                                                               <  1:9331192/150‑1 (MQ=255)
                     gCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTc                                                                                                                             >  1:10011547/1‑150 (MQ=255)
                      cGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCAAGCTCGCGTGCCGCTTcc                                                                                                                            >  2:2899108/1‑150 (MQ=255)
                      cGATTGCAGCGCAGTAGCGTTTCCTGGAACAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTcc                                                                                                                            <  2:7992550/150‑1 (MQ=255)
                       gATTGCAGCGCAGTAGCGTTTCCTGGATCATCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTcca                                                                                                                           >  2:944622/1‑150 (MQ=255)
                       gATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTcca                                                                                                                           >  2:8961550/1‑150 (MQ=255)
                             aGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCGTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGAGTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGc                                                                                                                     >  2:8181527/1‑150 (MQ=255)
                                 cAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGtt                                                                                                                 <  2:9893706/150‑1 (MQ=255)
                                 cAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGtt                                                                                                                 <  2:2344743/150‑1 (MQ=255)
                                  aGTAGCGTTTCCTGGATCAGCCCGGCCAGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTg                                                                                                                >  1:10427228/1‑150 (MQ=255)
                                      gCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGgcgc                                                                                                            <  2:419551/150‑1 (MQ=255)
                                       cGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCt                                                                                                           >  2:653066/1‑150 (MQ=255)
                                         tttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTCAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGa                                                                                                         >  2:3790521/1‑150 (MQ=255)
                                         tttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGa                                                                                                         <  1:1376245/150‑1 (MQ=255)
                                         tttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCGGa                                                                                                         >  2:624326/1‑150 (MQ=255)
                                          ttCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAg                                                                                                        <  2:10149966/150‑1 (MQ=255)
                                              tGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAg                                                                                                    <  2:1928469/150‑1 (MQ=255)
                                              tGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAg                                                                                                    >  2:6024648/1‑150 (MQ=255)
                                               ggATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCCCCAGGCGGGTGGGGCGGGGGTAGt                                                                                                   >  1:5431864/1‑150 (MQ=255)
                                               ggATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCCGGCGGGTGGCGCTGAGGTAGt                                                                                                   >  2:8375230/1‑150 (MQ=255)
                                                gATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTc                                                                                                  >  2:8553164/1‑150 (MQ=255)
                                                    aGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTgg                                                                                              >  1:3067277/1‑150 (MQ=255)
                                                       ccGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAg                                                                                           >  2:3207002/1‑150 (MQ=255)
                                                           ccTGGCGCAGACGCGCCTCCTCCCGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTcgc                                                                                       <  1:4482810/150‑1 (MQ=255)
                                                            cTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTcgcg                                                                                      <  2:7617424/150‑1 (MQ=255)
                                                             tGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCgg                                                                                     <  1:3301613/150‑1 (MQ=255)
                                                             tGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCgg                                                                                     >  2:10685797/1‑150 (MQ=255)
                                                             tGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTACGGGTAGTCGCgg                                                                                     >  2:4855274/1‑150 (MQ=255)
                                                                   aGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGGGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGgcgcgc                                                                               >  2:5780095/1‑150 (MQ=255)
                                                                    gACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCa                                                                              >  1:2560498/1‑150 (MQ=255)
                                                                        cgcCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAg                                                                          <  1:8358719/150‑1 (MQ=255)
                                                                          cctcctCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAggg                                                                        >  1:2392224/1‑150 (MQ=255)
                                                                          cctcctCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAggg                                                                        >  2:4267457/1‑150 (MQ=255)
                                                                            tcctccAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGcc                                                                      <  2:1374286/150‑1 (MQ=255)
                                                                            tcctccAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGcc                                                                      >  2:4673320/1‑150 (MQ=255)
                                                                                 cAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTCGTCGTGGGCGTCGCGGGCGCGCAACAGGGCCTGccc                                                                 >  2:10064039/1‑150 (MQ=255)
                                                                                      cGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAACCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCggggg                                                            <  1:171004/150‑1 (MQ=255)
                                                                                          cTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAg                                                        <  2:7309918/150‑1 (MQ=255)
                                                                                          cTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAg                                                        >  2:7786116/1‑150 (MQ=255)
                                                                                               agcccagGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTa                                                   <  1:5184602/150‑1 (MQ=255)
                                                                                                gcccagGCGCTTCTGCTCCTGGTAACCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTAc                                                  <  2:6060930/150‑1 (MQ=255)
                                                                                                 cccagGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTAcg                                                 <  2:8291354/150‑1 (MQ=255)
                                                                                                  cgagGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTGGCGCTCGTTTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTAcgc                                                <  1:554577/148‑1 (MQ=255)
                                                                                                  ccagGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTAcgc                                                <  2:9315213/150‑1 (MQ=255)
                                                                                                   cagGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTAcgcg                                               >  2:5058472/1‑150 (MQ=255)
                                                                                                    agGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGa                                              <  1:7216108/150‑1 (MQ=255)
                                                                                                     ggCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGAc                                             <  1:4019909/150‑1 (MQ=255)
                                                                                                       cgcTTCTGCTCCTGTTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTg                                           <  1:4880861/150‑1 (MQ=255)
                                                                                                         cTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGc                                         >  1:6395106/1‑150 (MQ=255)
                                                                                                         cTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGc                                         >  2:4513972/1‑150 (MQ=255)
                                                                                                          ttCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCg                                        <  2:7424724/150‑1 (MQ=255)
                                                                                                           tCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGt                                       >  1:956617/1‑150 (MQ=255)
                                                                                                           tCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGt                                       >  1:9840783/1‑150 (MQ=255)
                                                                                                            cTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTc                                      >  1:3874821/1‑150 (MQ=255)
                                                                                                              gCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTa                                    >  1:6050396/1‑150 (MQ=255)
                                                                                                              gCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTa                                    >  2:5351293/1‑150 (MQ=255)
                                                                                                               cTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAg                                   >  1:6839981/1‑150 (MQ=255)
                                                                                                                tCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGt                                  >  2:3866868/1‑150 (MQ=255)
                                                                                                                 ccTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGtt                                 <  1:7895647/150‑1 (MQ=255)
                                                                                                                  cTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTc                                <  2:9931918/150‑1 (MQ=255)
                                                                                                                  cTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTc                                <  2:2160779/150‑1 (MQ=255)
                                                                                                                   tGGTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTca                               <  2:1110232/150‑1 (MQ=255)
                                                                                                                    ggTAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTcaa                              >  1:1950715/1‑150 (MQ=255)
                                                                                                                      tAAACCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTcaaca                            >  2:6295384/1‑150 (MQ=255)
                                                                                                                        aaCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACaaa                          <  1:8885180/150‑1 (MQ=255)
                                                                                                                         aCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAAt                         >  2:2764705/1‑150 (MQ=255)
                                                                                                                         aCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAAt                         >  2:9061443/1‑150 (MQ=255)
                                                                                                                         aCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAAt                         <  1:8224216/150‑1 (MQ=255)
                                                                                                                         aCCTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAAt                         <  1:1511585/150‑1 (MQ=255)
                                                                                                                          ccTCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATg                        >  2:795558/1‑150 (MQ=255)
                                                                                                                            tCGTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGa                      >  2:178300/1‑150 (MQ=255)
                                                                                                                              gTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAg                    <  1:9985770/150‑1 (MQ=255)
                                                                                                                              gTGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAg                    >  1:6669611/1‑150 (MQ=255)
                                                                                                                               tGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGc                   >  2:9830867/1‑150 (MQ=255)
                                                                                                                               tGGGCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGc                   >  1:6991012/1‑150 (MQ=255)
                                                                                                                                gggCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCa                  <  1:7346561/150‑1 (MQ=255)
                                                                                                                                  gCCTCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCAtt                >  1:9061769/1‑150 (MQ=255)
                                                                                                                                     tCGCGCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCAt             <  2:9804275/150‑1 (MQ=255)
                                                                                                                                      cgcgCTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCAtc            >  1:629735/1‑150 (MQ=255)
                                                                                                                                         gcTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCAtcgct         <  2:9674426/150‑1 (MQ=255)
                                                                                                                                         gcTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCAtcgct         >  2:2622468/1‑150 (MQ=255)
                                                                                                                                          cTCGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCAtcgctc        >  1:1535504/1‑150 (MQ=255)
                                                                                                                                            cGATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCAtcgctcgc      >  2:10448575/1‑150 (MQ=255)
                                                                                                                                             gATTTCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATcgctcgct     <  1:3513366/150‑1 (MQ=255)
                                                                                                                                               tttCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTg   >  1:10084224/1‑150 (MQ=255)
                                                                                                                                                ttCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTGt  <  2:6996610/150‑1 (MQ=255)
                                                                                                                                                ttCGGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGGAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTGt  <  1:2767165/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGATCGACCTGCCGGTAGTAGCGATTGCAGCGCAGTAGCGTTTCCTGGATCAGCCCGGCCTGGCGCAGACGCGCCTCCTCCAGTCCGGCCTCCCAGCCCAGGCGCTTCTGCTCCTGGTAAACCTCGTGGGCCTCGCGCTCGATTTCCGCGGCACGCTCGCGTGCCGCTTCCACCAGGCGGTTGGCGCTGAGGTAGTCCTGGTAGTCGCGGGCGCGCAACAGGGCCTGCCCGGGGGCCAGCCGTACGCGACTGGCGTCTAGTTCAACAAATGGAAGCATTCATCGCTCGCTTGT  >  NC_008463/3769033‑3769326

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