Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,126,576 A→C 14.3% 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,5760AC14.3% 471.0 / 85.7 215T35P (ACC→CCC) flgBflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (94/90);  new base C (19/12);  total (113/102)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 218 total aligned reads displayed.

AGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAG  >  CP009273/1126429‑1126723
                                                                                                                                                   |                                                                                                                                                   
 ggCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATAc                                                                                                                                                    >  2:1126845/1‑148 (MQ=255)
    aTTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATAcccc                                                                                                                                                 <  1:1003856/148‑1 (MQ=255)
     tttAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCt                                                                                                                                                <  1:736107/148‑1 (MQ=255)
      ttAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATTCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGACAGGACGTGCTGGCAGCAAACATCTCCAATGCCGATACCCCTg                                                                                                                                               >  1:401015/1‑148 (MQ=255)
      ttAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTg                                                                                                                                               >  1:2600519/1‑148 (MQ=255)
       tAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTgg                                                                                                                                              <  1:2803369/148‑1 (MQ=255)
          gCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTa                                                                                                                                           <  1:186210/148‑1 (MQ=255)
          gCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTa                                                                                                                                           >  1:3051092/1‑148 (MQ=255)
          gCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTa                                                                                                                                           >  1:2069207/1‑148 (MQ=255)
           cTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTAt                                                                                                                                          <  1:934942/148‑1 (MQ=255)
            tGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATc                                                                                                                                         <  1:2212333/148‑1 (MQ=255)
             gATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCa                                                                                                                                        >  1:176391/1‑148 (MQ=255)
                      aTTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGa                                                                                                                               <  1:279275/148‑1 (MQ=255)
                        tttGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGata                                                                                                                             <  2:3583187/148‑1 (MQ=255)
                        tttGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGata                                                                                                                             <  1:1521663/148‑1 (MQ=255)
                         ttGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGatat                                                                                                                            <  2:2343289/148‑1 (MQ=255)
                          tGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATc                                                                                                                           <  1:21449/148‑1 (MQ=255)
                           gATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCg                                                                                                                          <  1:385102/148‑1 (MQ=255)
                            aTACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGa                                                                                                                         >  2:2979189/1‑148 (MQ=255)
                            aTACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGa                                                                                                                         >  1:2655024/1‑148 (MQ=255)
                              aCCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGAtt                                                                                                                       >  1:1881508/1‑148 (MQ=255)
                              aCCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGAtt                                                                                                                       <  2:572792/148‑1 (MQ=255)
                                 tGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTg                                                                                                                    >  2:2273777/1‑148 (MQ=255)
                                  gCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGc                                                                                                                   <  2:1720336/148‑1 (MQ=255)
                                      aggagATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGt                                                                                                               >  1:1617945/1‑148 (MQ=255)
                                       ggagATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTg                                                                                                              >  1:2093812/1‑148 (MQ=255)
                                        gagaTATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGa                                                                                                             >  2:671934/1‑148 (MQ=255)
                                             atGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTa                                                                                                        >  1:3442362/1‑148 (MQ=255)
                                             atGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTa                                                                                                        >  1:3727307/1‑148 (MQ=255)
                                              tGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaa                                                                                                       >  1:3406402/1‑148 (MQ=255)
                                              tGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaa                                                                                                       <  1:2551744/148‑1 (MQ=255)
                                               gCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaa                                                                                                      >  2:408069/1‑148 (MQ=255)
                                                 tCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaa                                                                                                    >  2:897615/1‑148 (MQ=255)
                                                  cGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaaa                                                                                                   >  1:1315520/1‑148 (MQ=255)
                                                  cGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTaaaaaa                                                                                                   >  1:1668452/1‑148 (MQ=255)
                                                   gATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAg                                                                                                  <  1:480848/148‑1 (MQ=255)
                                                   gATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAg                                                                                                  <  1:646758/148‑1 (MQ=255)
                                                   gATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAg                                                                                                  <  1:3139748/148‑1 (MQ=255)
                                                   gATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAg                                                                                                  <  1:2388301/148‑1 (MQ=255)
                                                    aTAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGCTATCGATTTTGCCAGTGAACTTAAAAAAGt                                                                                                 >  1:322019/1‑148 (MQ=255)
                                                     tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTc                                                                                                <  1:1163417/148‑1 (MQ=255)
                                                      aaGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCa                                                                                               <  2:2448124/148‑1 (MQ=255)
                                                        gCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATg                                                                                             >  1:2873166/1‑148 (MQ=255)
                                                           cGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCaa                                                                                          >  1:905278/1‑148 (MQ=255)
                                                              cgccgcCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGt                                                                                       <  2:1518730/148‑1 (MQ=255)
                                                                   ccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACg                                                                                  >  1:3054299/1‑148 (MQ=255)
                                                                   ccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTGAAAAAAGTCATGCAACGTGGACg                                                                                  >  2:2548761/1‑148 (MQ=255)
                                                                    cTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACgg                                                                                 >  1:2219395/1‑148 (MQ=255)
                                                                     ttACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACggg                                                                                >  2:3741202/1‑148 (MQ=255)
                                                                     ttAAGTATTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCACCAAACAACGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTAGCCAGTGAACTTAAAATAGTCATGCAACGTGGACggg                                                                                >  2:3022017/1‑148 (MQ=255)
                                                                      tACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGa                                                                               >  2:1951837/1‑148 (MQ=255)
                                                                      tACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGa                                                                               >  2:3351642/1‑148 (MQ=255)
                                                                          ttttCAACTAGAGGCTCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGTCAGCAAACATCGCCAATGCCGATACCACTGGTTATAAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAAAGTGGACGGGATGCa                                                                           <  2:76005/148‑1 (MQ=255)
                                                                          ttttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCa                                                                           <  1:2667552/148‑1 (MQ=255)
                                                                          ttttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCa                                                                           >  1:2828993/1‑148 (MQ=255)
                                                                          ttttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCa                                                                           >  2:1738066/1‑148 (MQ=255)
                                                                           tttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCaa                                                                          <  2:1564116/148‑1 (MQ=255)
                                                                           tttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCaa                                                                          <  2:2939253/148‑1 (MQ=255)
                                                                           tttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCaa                                                                          >  1:3115767/1‑148 (MQ=255)
                                                                           tttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCaa                                                                          >  1:628190/1‑148 (MQ=255)
                                                                            ttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAAc                                                                         >  2:2249251/1‑148 (MQ=255)
                                                                                acaaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAgt                                                                     <  2:472664/148‑1 (MQ=255)
                                                                                  aaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAgtgt                                                                   <  2:3725973/148‑1 (MQ=255)
                                                                                   agagGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGtgtg                                                                  <  1:1349019/148‑1 (MQ=255)
                                                                                         gcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCa                                                            <  2:260672/148‑1 (MQ=255)
                                                                                         gcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCa                                                            <  2:1661038/148‑1 (MQ=255)
                                                                                           tCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAACAAGTCATGCAACGGGGACGGGATGCAACCAGTGTGGGTGCAcg                                                          >  2:2080130/1‑147 (MQ=255)
                                                                                            cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTg                                                         >  2:3818919/1‑148 (MQ=255)
                                                                                              aTCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                       >  2:3153370/1‑148 (MQ=255)
                                                                                               tCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACg                                                      >  1:1253632/1‑148 (MQ=255)
                                                                                               tCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACg                                                      <  1:1904837/148‑1 (MQ=255)
                                                                                               tCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACg                                                      >  1:2280817/1‑148 (MQ=255)
                                                                                                  gcgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATg                                                   <  2:1253443/148‑1 (MQ=255)
                                                                                                      gcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCt                                               <  2:564127/148‑1 (MQ=255)
                                                                                                         cAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCaa                                            >  1:3278143/1‑148 (MQ=255)
                                                                                                               cAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAac                                      <  1:1795562/148‑1 (MQ=255)
                                                                                                               cAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAac                                      <  1:3687792/148‑1 (MQ=255)
                                                                                                                 ggAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAacac                                    >  1:1065810/1‑148 (MQ=255)
                                                                                                                 ggAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAacac                                    >  1:2947413/1‑148 (MQ=255)
                                                                                                                  taaGTGCTGGCAGCAAACATCGCCAATGCCGATCCCCCTTTTTATCAGTCGCGCGATATCTATTTTGCCATTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGTTTTCACTGACGATGACCTCAACGCAAcaca                                   <  2:2898531/147‑1 (MQ=255)
                                                                                                                  gAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAAcaca                                   <  2:970409/148‑1 (MQ=255)
                                                                                                                   aaGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACAt                                  >  2:932939/1‑148 (MQ=255)
                                                                                                                     gTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTc                                <  1:3311144/148‑1 (MQ=255)
                                                                                                                     gTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTc                                <  2:3582458/148‑1 (MQ=255)
                                                                                                                          ggCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgcg                           >  2:2922508/1‑148 (MQ=255)
                                                                                                                            cagcaAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCa                         >  1:2872110/1‑148 (MQ=255)
                                                                                                                                aaaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcg                     >  1:521671/1‑148 (MQ=255)
                                                                                                                                aaaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcg                     >  1:3089080/1‑148 (MQ=255)
                                                                                                                                 aaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                    >  2:2590093/1‑148 (MQ=255)
                                                                                                                                 aaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                    >  1:203452/1‑148 (MQ=255)
                                                                                                                                   cATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCACCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                  >  2:2818444/1‑148 (MQ=255)
                                                                                                                                   cATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                  >  1:1908479/1‑148 (MQ=255)
                                                                                                                                      cGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACg               <  2:1315520/148‑1 (MQ=255)
                                                                                                                                      cGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACg               >  1:1687234/1‑148 (MQ=255)
                                                                                                                                          aaTGCCGATACCCCTGGTTATGAGGCGCGCGATATCGATTTTGCCAGTGAACGTTAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGTCCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctc           <  2:401015/148‑1 (MQ=255)
                                                                                                                                          aaTGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctc           <  1:3818903/148‑1 (MQ=255)
                                                                                                                                           aTGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctcc          <  1:3775996/148‑1 (MQ=255)
                                                                                                                                            tGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctcct         >  1:3053357/1‑148 (MQ=255)
                                                                                                                                            tGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctcct         >  2:2813166/1‑148 (MQ=255)
                                                                                                                                                   aCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAg  >  2:1302235/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAG  >  CP009273/1126429‑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: