Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,126,576 A→C 32.4% T35P (ACC→CCC)  flgB → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,126,5760AC32.4% 188.3 / 252.4 229T35P (ACC→CCC) flgBflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (70/85);  new base C (41/33);  total (111/118)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 230 total aligned reads displayed.

CATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAG  >  CP009273/1126432‑1126723
                                                                                                                                                |                                                                                                                                                   
cATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATAccc                                                                                                                                                  >  1:3548433/1‑148 (MQ=255)
cATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATAccc                                                                                                                                                  <  1:2875684/148‑1 (MQ=255)
 aTTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATAcccc                                                                                                                                                 <  2:3482952/148‑1 (MQ=255)
  tttAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCt                                                                                                                                                <  1:1617347/148‑1 (MQ=255)
   ttAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTg                                                                                                                                               >  2:144195/1‑148 (MQ=255)
    tAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTgg                                                                                                                                              >  1:2089476/1‑148 (MQ=255)
     aaGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGt                                                                                                                                             <  2:2576137/148‑1 (MQ=255)
      aGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGtt                                                                                                                                            >  1:867396/1‑148 (MQ=255)
      aGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGtt                                                                                                                                            >  1:1293587/1‑148 (MQ=255)
          gATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCa                                                                                                                                        <  2:1664036/148‑1 (MQ=255)
            tGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAgg                                                                                                                                      <  2:2057994/148‑1 (MQ=255)
               cAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGgcgc                                                                                                                                   >  1:2260698/1‑148 (MQ=255)
                  aaTTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGcgcgcg                                                                                                                                >  2:2509409/1‑148 (MQ=255)
                   aTTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGa                                                                                                                               <  2:2464692/148‑1 (MQ=255)
                      ttGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGatat                                                                                                                            <  1:716041/148‑1 (MQ=255)
                        gATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCg                                                                                                                          <  1:3425156/148‑1 (MQ=255)
                         aTACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGa                                                                                                                         <  2:824749/148‑1 (MQ=255)
                         aTACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGa                                                                                                                         >  1:1389323/1‑148 (MQ=255)
                          tACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGAt                                                                                                                        <  1:768994/148‑1 (MQ=255)
                           aCCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGAtt                                                                                                                       <  1:2090663/148‑1 (MQ=255)
                           aCCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGAtt                                                                                                                       >  1:2385596/1‑148 (MQ=255)
                           aCCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGAtt                                                                                                                       <  2:2556380/148‑1 (MQ=255)
                           aCCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGAtt                                                                                                                       >  2:1110226/1‑148 (MQ=255)
                            ccTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGAttt                                                                                                                      >  2:1395086/1‑148 (MQ=255)
                             cTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGAtttt                                                                                                                     >  1:465940/1‑148 (MQ=255)
                             cTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGAtttt                                                                                                                     <  2:3353536/148‑1 (MQ=255)
                             cTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGAtttt                                                                                                                     >  1:2454543/1‑148 (MQ=255)
                                cGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGcc                                                                                                                  >  1:2992437/1‑148 (MQ=255)
                                    ggagATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTg                                                                                                              <  1:3487167/148‑1 (MQ=255)
                                      agaTATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGaa                                                                                                            <  1:2225462/148‑1 (MQ=255)
                                      agaTATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGaa                                                                                                            >  2:2472272/1‑148 (MQ=255)
                                            gCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCACCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaa                                                                                                      >  2:2411018/1‑148 (MQ=255)
                                              tCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaa                                                                                                    <  2:734742/148‑1 (MQ=255)
                                              tCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaa                                                                                                    <  2:734653/148‑1 (MQ=255)
                                              tCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaa                                                                                                    <  2:31403/148‑1 (MQ=255)
                                               cGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaaa                                                                                                   >  2:1041918/1‑148 (MQ=255)
                                               cGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaaa                                                                                                   <  2:3442745/148‑1 (MQ=255)
                                                 aTAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGt                                                                                                 >  1:986724/1‑148 (MQ=255)
                                                 aTAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGt                                                                                                 >  1:1115364/1‑148 (MQ=255)
                                                  tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTc                                                                                                <  1:2118543/148‑1 (MQ=255)
                                                  tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTc                                                                                                >  1:1507932/1‑148 (MQ=255)
                                                  tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTc                                                                                                >  1:64358/1‑148 (MQ=255)
                                                  tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTc                                                                                                >  1:3274746/1‑148 (MQ=255)
                                                    aGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTCAAAAAGTCAt                                                                                              >  2:2723464/1‑148 (MQ=255)
                                                      cTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGGTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGc                                                                                            <  2:1461967/148‑1 (MQ=255)
                                                          aCGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACg                                                                                        <  1:1782288/148‑1 (MQ=255)
                                                          aCGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACg                                                                                        <  1:2045264/148‑1 (MQ=255)
                                                           cgccgcCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGt                                                                                       >  2:3195970/1‑148 (MQ=255)
                                                            gccgccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTg                                                                                      <  1:824151/148‑1 (MQ=255)
                                                            gccgccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTg                                                                                      <  1:2606305/148‑1 (MQ=255)
                                                             ccgccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTgg                                                                                     <  2:1692/148‑1 (MQ=255)
                                                               gccTTACGTTTTCATCAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGAc                                                                                   <  1:500245/148‑1 (MQ=255)
                                                               gccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGAc                                                                                   >  1:2406120/1‑148 (MQ=255)
                                                                 cTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACgg                                                                                 >  1:2189245/1‑148 (MQ=255)
                                                                     cGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATg                                                                             <  1:3181174/148‑1 (MQ=255)
                                                                         ttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAAc                                                                         <  2:663483/148‑1 (MQ=255)
                                                                              caaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAgtg                                                                    >  1:600827/1‑148 (MQ=255)
                                                                              caaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAgtg                                                                    >  1:732357/1‑148 (MQ=255)
                                                                                 gagGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTgg                                                                 >  2:521828/1‑148 (MQ=255)
                                                                                    gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTg                                                              >  1:2654472/1‑148 (MQ=255)
                                                                                       cTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCAc                                                           <  2:2312099/148‑1 (MQ=255)
                                                                                        tCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACt                                                          >  1:538506/1‑148 (MQ=255)
                                                                                         cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTg                                                         >  2:429644/1‑148 (MQ=255)
                                                                                         cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTg                                                         >  1:1969306/1‑148 (MQ=255)
                                                                                         cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAATAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTg                                                         >  2:30777/1‑148 (MQ=255)
                                                                                         cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTg                                                         <  2:2556639/148‑1 (MQ=255)
                                                                                           aTCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                       >  2:72377/1‑148 (MQ=255)
                                                                                            tCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACg                                                      <  2:939865/148‑1 (MQ=255)
                                                                                             cTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGa                                                     <  2:2252244/148‑1 (MQ=255)
                                                                                             cTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGa                                                     <  1:1707962/148‑1 (MQ=255)
                                                                                               gcgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATg                                                   <  2:1201831/148‑1 (MQ=255)
                                                                                               gcgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATg                                                   <  2:1068674/148‑1 (MQ=255)
                                                                                                cgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGa                                                  >  2:3006774/1‑148 (MQ=255)
                                                                                                 gcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAc                                                 <  2:895814/148‑1 (MQ=255)
                                                                                                  cgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAcc                                                <  1:1204532/148‑1 (MQ=255)
                                                                                                  cgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAcc                                                >  1:369326/1‑148 (MQ=255)
                                                                                                     ccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCa                                             <  1:155865/148‑1 (MQ=255)
                                                                                                     ccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCa                                             >  2:1870673/1‑148 (MQ=255)
                                                                                                      cAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCaa                                            >  1:1885062/1‑148 (MQ=255)
                                                                                                        gCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACg                                          <  2:2671990/148‑1 (MQ=255)
                                                                                                             aGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAaca                                     <  2:1425531/148‑1 (MQ=255)
                                                                                                             aGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAaca                                     <  2:2801124/148‑1 (MQ=255)
                                                                                                              ggAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAacag                                    <  1:951062/148‑2 (MQ=255)
                                                                                                              ggAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAacac                                    <  1:2462860/148‑1 (MQ=255)
                                                                                                               gAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAAcaca                                   <  2:3263099/148‑1 (MQ=255)
                                                                                                                aaGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACAt                                  <  1:875535/148‑1 (MQ=255)
                                                                                                                  gTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTc                                >  2:2752601/1‑148 (MQ=255)
                                                                                                                        gcagcaAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgcgc                          >  2:428986/1‑148 (MQ=255)
                                                                                                                        gcagcaAACATCGCCAATGCCGATACCCCTGGTTATCGGGCGCGCGAGATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgcgc                          >  2:1452437/1‑148 (MQ=255)
                                                                                                                         cagcaAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCa                         <  1:2418261/148‑1 (MQ=255)
                                                                                                                            caAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgc                      >  2:492048/1‑148 (MQ=255)
                                                                                                                            caAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgc                      >  2:361065/1‑148 (MQ=255)
                                                                                                                              aaCATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                    <  1:2733867/148‑1 (MQ=255)
                                                                                                                              aaCATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                    <  1:1035126/148‑1 (MQ=255)
                                                                                                                                  tCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGAc                >  2:3157369/1‑148 (MQ=255)
                                                                                                                                   cGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACg               >  2:393032/1‑148 (MQ=255)
                                                                                                                                       aaTGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctc           <  1:3147995/148‑1 (MQ=255)
                                                                                                                                          gCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTa        >  1:2022326/1‑148 (MQ=255)
                                                                                                                                          gCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTa        <  2:127452/148‑1 (MQ=255)
                                                                                                                                              aTACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGc    <  2:1839100/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAG  >  CP009273/1126432‑1126723

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: