Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,007,090 G→T 20.0% V127V (GTG→GTT fliF → flagellar basal‑body MS‑ring and collar protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,007,0900GT20.0% 976.2 / 163.9 354V127V (GTG→GTTfliFflagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base G (140/143);  new base T (32/39);  total (172/182)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 355 total aligned reads displayed.

AAGCCAGCGGCGCTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAAT  >  CP009273/2006942‑2007237
                                                                                                                                                    |                                                                                                                                                   
  gCCAGCGGCGCTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGaa                                                                                                                                                   <  2:3875982/149‑1 (MQ=255)
   ccAGCGGCGCTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAAt                                                                                                                                                  <  2:1518264/149‑1 (MQ=255)
   ccAGCGGCGCTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAAt                                                                                                                                                  <  1:4710757/149‑1 (MQ=255)
   ccAGCGGCGCTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAAt                                                                                                                                                  <  2:4332924/149‑1 (MQ=255)
     aGCGGCGCTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTa                                                                                                                                                <  1:2375117/149‑1 (MQ=255)
         gcgcTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAg                                                                                                                                            >  1:4471235/1‑149 (MQ=255)
          cgcTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGc                                                                                                                                           <  1:760734/149‑1 (MQ=255)
           gcTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCg                                                                                                                                          >  2:785181/1‑149 (MQ=255)
            cTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCgg                                                                                                                                         >  2:2736734/1‑149 (MQ=255)
            cTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCgg                                                                                                                                         >  1:645185/1‑149 (MQ=255)
            cTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCgg                                                                                                                                         >  2:1056219/1‑149 (MQ=255)
            cTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCgg                                                                                                                                         >  1:1078447/1‑149 (MQ=255)
              aTTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGACAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGgc                                                                                                                                       >  2:4442307/1‑149 (MQ=255)
              aTTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGgc                                                                                                                                       >  1:3572339/1‑149 (MQ=255)
               ttGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGgcg                                                                                                                                      >  2:175984/1‑149 (MQ=255)
                  aaGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTg                                                                                                                                   <  1:1170751/149‑1 (MQ=255)
                  aaGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTg                                                                                                                                   <  2:4709605/149‑1 (MQ=255)
                   aGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTgg                                                                                                                                  >  2:2530882/1‑149 (MQ=255)
                       ccGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAgg                                                                                                                              <  1:2641226/149‑1 (MQ=255)
                          gCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGa                                                                                                                           <  1:1766078/149‑1 (MQ=255)
                           cAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAg                                                                                                                          <  1:1868458/149‑1 (MQ=255)
                            aGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGc                                                                                                                         >  2:1864102/1‑149 (MQ=255)
                              aTAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCtt                                                                                                                       >  1:1019338/1‑149 (MQ=255)
                                aaaGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTtc                                                                                                                     >  1:3594304/1‑149 (MQ=255)
                                 aaGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTtct                                                                                                                    <  2:4100284/149‑1 (MQ=255)
                                 aaGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTtct                                                                                                                    >  2:712668/1‑149 (MQ=255)
                                 aaGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTtct                                                                                                                    <  2:3100955/149‑1 (MQ=255)
                                      cACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTAc                                                                                                               >  1:4605388/1‑149 (MQ=255)
                                      cACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTAc                                                                                                               >  1:408691/1‑149 (MQ=255)
                                       aCGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTAcc                                                                                                              <  2:3354017/149‑1 (MQ=255)
                                          aaCTGCGTCTGCGCCTGGCACAACAAGGTATGCCAAAAGGAGGCGCGGTCGTTTTCGAACTGCTTTATCATGAAAAGTTTGGTATCATCCATTTAAGAGAAAAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTATCTTACCATc                                                                                                           >  1:1244301/1‑149 (MQ=255)
                                            cTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGa                                                                                                         <  1:1157613/149‑1 (MQ=255)
                                               cGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAAc                                                                                                      <  1:2687905/149‑1 (MQ=255)
                                                gTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACt                                                                                                     >  1:1732281/1‑149 (MQ=255)
                                                 tCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTa                                                                                                    >  1:502380/1‑149 (MQ=255)
                                                  cTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTAt                                                                                                   >  1:4214722/1‑149 (MQ=255)
                                                    gcgcCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCg                                                                                                 <  2:3489338/149‑1 (MQ=255)
                                                       ccTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCATGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAAATATCGGcc                                                                                              >  2:2386745/1‑149 (MQ=255)
                                                       ccTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGcc                                                                                              <  2:4122482/149‑1 (MQ=255)
                                                        cTGGCACAACAAGGTTTTCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGccc                                                                                             <  1:250838/149‑1 (MQ=255)
                                                        cTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGccc                                                                                             <  1:3700044/149‑1 (MQ=255)
                                                          ggCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCgg                                                                                           >  1:695332/1‑149 (MQ=255)
                                                           gCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGt                                                                                          >  2:2979005/1‑149 (MQ=255)
                                                           gCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGt                                                                                          <  2:1294956/149‑1 (MQ=255)
                                                           gCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTGTGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGt                                                                                          <  2:2008194/149‑1 (MQ=255)
                                                               aacaaGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGACACTATCGGCCCGGTaaaa                                                                                      >  1:829912/1‑149 (MQ=255)
                                                               aacaaGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTaaaa                                                                                      >  1:3055928/1‑149 (MQ=255)
                                                               aacaaGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTaaaa                                                                                      <  2:3315412/149‑1 (MQ=255)
                                                                acaaGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAg                                                                                     >  1:3409428/1‑149 (MQ=255)
                                                                   aGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAgggg                                                                                  >  2:4654303/1‑149 (MQ=255)
                                                                   aGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAgggg                                                                                  >  2:3625016/1‑149 (MQ=255)
                                                                    ggTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGgc                                                                                 >  1:3106877/1‑149 (MQ=255)
                                                                        tGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGgcgcgc                                                                             <  2:374006/149‑1 (MQ=255)
                                                                        tGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGgcgcgc                                                                             >  1:972804/1‑149 (MQ=255)
                                                                           cAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTa                                                                          >  2:1398278/1‑149 (MQ=255)
                                                                             aaaGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACa                                                                        >  2:305534/1‑149 (MQ=255)
                                                                              aaGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGGGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACAt                                                                       >  2:1565075/1‑149 (MQ=255)
                                                                               aGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATc                                                                      <  2:1287410/149‑1 (MQ=255)
                                                                                ggcggcGCGGTCTGTTTCGAACTGGTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCTGCGAACAGGTTATTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCt                                                                     <  2:4240585/149‑1 (MQ=255)
                                                                                ggcggcGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCt                                                                     <  2:1575939/149‑1 (MQ=255)
                                                                                 gcggcgCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTg                                                                    <  2:2243977/149‑1 (MQ=255)
                                                                                  cggcgCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTgg                                                                   <  2:1141701/149‑1 (MQ=255)
                                                                                      gcgGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAAt                                                               <  2:864846/149‑1 (MQ=255)
                                                                                      gcgGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAAt                                                               >  2:3290723/1‑149 (MQ=255)
                                                                                      gcgGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAAt                                                               <  2:609215/149‑1 (MQ=255)
                                                                                         gtcggtTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGcc                                                            >  1:4164272/1‑149 (MQ=255)
                                                                                          tcggtTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCg                                                           >  1:272230/1‑149 (MQ=255)
                                                                                           cggtTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGa                                                          <  1:2270383/149‑1 (MQ=255)
                                                                                             gtTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGTGGCGCGCGTACATCTGGCAATGCCGaaa                                                        >  1:2872264/1‑149 (MQ=255)
                                                                                             gtTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGaaa                                                        <  2:2242087/149‑1 (MQ=255)
                                                                                                tCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCg                                                     >  2:2616296/1‑149 (MQ=255)
                                                                                                  gAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTc                                                   <  1:3257629/149‑1 (MQ=255)
                                                                                                   aaCTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCt                                                  <  2:4032533/149‑1 (MQ=255)
                                                                                                      tGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTa                                               <  1:86873/149‑1 (MQ=255)
                                                                                                         ttGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTc                                            >  1:3838960/1‑149 (MQ=255)
                                                                                                         ttGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTc                                            >  2:4399793/1‑149 (MQ=255)
                                                                                                         ttGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTc                                            >  2:4544889/1‑149 (MQ=255)
                                                                                                          tGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCg                                           <  2:1339362/149‑1 (MQ=255)
                                                                                                                 gAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTACTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGACTATAATCGTCCGtca                                    >  2:183479/1‑147 (MQ=255)
                                                                                                                          ggTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATc                           <  1:1667336/149‑1 (MQ=255)
                                                                                                                           ggaTCAGCCAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATcc                          <  2:1628890/147‑1 (MQ=255)
                                                                                                                            tATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATccc                         <  2:3418460/149‑1 (MQ=255)
                                                                                                                                aGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTc                     <  1:1831312/149‑1 (MQ=255)
                                                                                                                                 gCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCt                    >  1:3159529/1‑149 (MQ=255)
                                                                                                                                 gCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCt                    <  2:3846816/149‑1 (MQ=255)
                                                                                                                                   cAGTTCAGCGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGc                  <  1:1358574/149‑1 (MQ=255)
                                                                                                                                   cAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGc                  >  2:3948981/1‑149 (MQ=255)
                                                                                                                                    aGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCa                 >  1:3146655/1‑149 (MQ=255)
                                                                                                                                     gTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCAt                >  1:4725356/1‑149 (MQ=255)
                                                                                                                                         aGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGt            >  2:4783176/1‑149 (MQ=255)
                                                                                                                                          gCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAATTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTg           >  1:555643/1‑149 (MQ=255)
                                                                                                                                          gCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTg           >  1:3381780/1‑149 (MQ=255)
                                                                                                                                          gCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTg           <  1:785181/149‑1 (MQ=255)
                                                                                                                                           cGAACAGGTTAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGa          <  2:3106877/149‑1 (MQ=255)
                                                                                                                                            gAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGAc         <  2:2302069/149‑1 (MQ=255)
                                                                                                                                             aaCAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACg        >  1:4523784/1‑149 (MQ=255)
                                                                                                                                               cAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGt      >  1:43504/1‑149 (MQ=255)
                                                                                                                                                  gTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTaaa   <  2:2052711/149‑1 (MQ=255)
                                                                                                                                                   tGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAAt  <  1:300939/149‑1 (MQ=255)
                                                                                                                                                   tGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAAt  >  1:3517354/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAGCCAGCGGCGCTATTGAAGTTCCGGCAGATAAAGTTCACGAACTGCGTCTGCGCCTGGCACAACAAGGTTTGCCAAAAGGCGGCGCGGTCGGTTTCGAACTGCTTGATCAGGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAAT  >  CP009273/2006942‑2007237

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: