Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,129,972 G→T 12.4% D239Y (GAT→TAT)  flgF → 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,129,9720GT12.4% 754.1 / 51.2 230D239Y (GAT→TAT) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base G (109/92);  new base T (21/8);  total (130/100)
Fisher's exact test for biased strand distribution p-value = 7.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

GGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGG  >  CP009273/1129825‑1130120
                                                                                                                                                   |                                                                                                                                                    
gggCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCgatga                                                                                                                                                     <  2:539380/149‑1 (MQ=255)
gggCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCgatga                                                                                                                                                     >  2:2825238/1‑149 (MQ=255)
gggCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGATTTGAAATGCAGATGAAGGTGATCAGAAGCGTCgatga                                                                                                                                                     >  2:468877/1‑149 (MQ=255)
 ggCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCgatgat                                                                                                                                                    <  2:3065443/149‑1 (MQ=255)
       tACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAATTGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCa                                                                                                                                              >  2:931913/1‑149 (MQ=255)
       tACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCGGATGAAGGTGATCAGCAGCGTCGAGGATAACGCa                                                                                                                                              >  2:1251718/1‑149 (MQ=255)
       tACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCa                                                                                                                                              >  2:1251589/1‑149 (MQ=255)
        aCTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAg                                                                                                                                             >  1:2592180/1‑149 (MQ=255)
          tGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGc                                                                                                                                           >  2:3589045/1‑149 (MQ=255)
          tGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGc                                                                                                                                           <  1:1332278/149‑1 (MQ=255)
          tGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGc                                                                                                                                           <  1:387974/149‑1 (MQ=255)
           gcaggcagATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGcc                                                                                                                                          <  1:116392/149‑1 (MQ=255)
             aggcagATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGt                                                                                                                                        >  2:1279925/1‑149 (MQ=255)
             aggcagATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGt                                                                                                                                        >  1:2948644/1‑149 (MQ=255)
              ggcagATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAATGTGATCAGCAGCGTCGATGATAACGCAGGCCGTg                                                                                                                                       >  2:2269960/1‑149 (MQ=255)
               gcagATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGc                                                                                                                                      >  2:2353888/1‑149 (MQ=255)
                       aaCCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCcaaccaac                                                                                                                              >  1:1405239/1‑149 (MQ=255)
                        aCCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAAct                                                                                                                             >  2:225151/1‑149 (MQ=255)
                        aCCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAAct                                                                                                                             >  2:1582071/1‑149 (MQ=255)
                          cTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAActgc                                                                                                                           >  2:2077208/1‑149 (MQ=255)
                           ttGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAGTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGGGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAActgct                                                                                                                          <  1:3691291/149‑1 (MQ=255)
                            tGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACTAActgctg                                                                                                                         <  1:612249/149‑1 (MQ=255)
                            tGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAActgctg                                                                                                                         <  1:3450681/149‑1 (MQ=255)
                             gCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGt                                                                                                                        <  2:570740/149‑1 (MQ=255)
                                tgtgATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGa                                                                                                                     >  2:731409/1‑149 (MQ=255)
                                tgtgATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGa                                                                                                                     >  2:731315/1‑149 (MQ=255)
                                 gtgATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGAt                                                                                                                    <  1:900930/149‑1 (MQ=255)
                                      gTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGtt                                                                                                               <  2:2932337/149‑1 (MQ=255)
                                       tCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTa                                                                                                              >  2:3516453/1‑149 (MQ=255)
                                        cGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTaa                                                                                                             >  2:2675593/1‑149 (MQ=255)
                                         gggggTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAAt                                                                                                            >  2:1719670/1‑149 (MQ=255)
                                            ggTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGa                                                                                                         >  1:3012736/1‑149 (MQ=255)
                                               tCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAg                                                                                                      <  1:2358433/149‑1 (MQ=255)
                                               tCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAg                                                                                                      <  2:695753/149‑1 (MQ=255)
                                                   gAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATa                                                                                                  <  1:3855782/149‑1 (MQ=255)
                                                   gAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATa                                                                                                  >  1:1890295/1‑149 (MQ=255)
                                                     aGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACa                                                                                                >  2:3381236/1‑149 (MQ=255)
                                                      ggCTGTAACGTCTATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACAt                                                                                               <  1:1907762/149‑1 (MQ=255)
                                                       gCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATg                                                                                              >  2:2533120/1‑149 (MQ=255)
                                                       gCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATg                                                                                              >  2:1937688/1‑149 (MQ=255)
                                                          gTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACa                                                                                           >  1:208071/1‑149 (MQ=255)
                                                              cGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGta                                                                                       >  2:2203878/1‑149 (MQ=255)
                                                               gTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGtat                                                                                      >  1:381222/1‑149 (MQ=255)
                                                                tCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAAAGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAATGATACATGAACAGtata                                                                                     >  1:1110515/1‑149 (MQ=255)
                                                                 cAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATaa                                                                                    <  1:3655347/149‑1 (MQ=255)
                                                                   aTGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTTCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGt                                                                                  >  1:1813297/1‑149 (MQ=255)
                                                                      ccGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACAGCTGTCGATGAGTTAATTGAAAGGATAAAAGAAAAATATAAGtttc                                                                               >  1:1775869/1‑147 (MQ=255)
                                                                      ccGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGc                                                                               <  1:802116/149‑1 (MQ=255)
                                                                       cGTTGCGGCAATGAGCGACATGAGTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGcc                                                                              <  2:119781/149‑1 (MQ=255)
                                                                        gTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGccc                                                                             >  1:3979123/1‑149 (MQ=255)
                                                                          tGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGa                                                                           <  2:1056233/149‑1 (MQ=255)
                                                                            cGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCTGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATg                                                                         <  2:3089280/149‑1 (MQ=255)
                                                                             ggCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATgc                                                                        >  2:895882/1‑149 (MQ=255)
                                                                             ggCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATgc                                                                        <  2:2650383/149‑1 (MQ=255)
                                                                             ggCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATgc                                                                        >  2:3985576/1‑149 (MQ=255)
                                                                              gCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATgcg                                                                       >  2:504889/1‑149 (MQ=255)
                                                                               cAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCCACCAACTGCTGTCGATGAGTTAATTGGAAGGCTACATGAAAAGGATAAGTTGGCCGATgcgc                                                                      >  1:2130508/1‑149 (MQ=255)
                                                                                aaTGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCa                                                                     >  2:1576164/1‑149 (MQ=255)
                                                                                       gACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCACAAGTTTAt                                                              <  2:694700/149‑1 (MQ=255)
                                                                                        aCATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATc                                                             >  2:2163066/1‑149 (MQ=255)
                                                                                        aCATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATc                                                             <  1:2827717/149‑1 (MQ=255)
                                                                                         cATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCg                                                            <  2:1993227/149‑1 (MQ=255)
                                                                                           tGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCggg                                                          >  2:1493740/1‑149 (MQ=255)
                                                                                            gATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGt                                                         >  2:1287169/1‑149 (MQ=255)
                                                                                              ttGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGCTACATGTCAAGTCGGCGTTGCCCGATGCGCCAGTTTATCGGGTcg                                                       >  2:2808649/1‑148 (MQ=255)
                                                                                               tGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTa                                                      >  1:837637/1‑149 (MQ=255)
                                                                                                gCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAt                                                     <  1:3428881/149‑1 (MQ=255)
                                                                                                   aGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGAAAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATggg                                                  >  1:1082708/1‑149 (MQ=255)
                                                                                                      gcgcgGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGAGAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGc                                               <  2:2065882/149‑1 (MQ=255)
                                                                                                         cggcgTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAAt                                            <  2:3191305/149‑1 (MQ=255)
                                                                                                           gcgTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCg                                          >  1:627984/1‑149 (MQ=255)
                                                                                                             gTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCa                                        >  2:1441804/1‑149 (MQ=255)
                                                                                                               tttGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAAt                                      >  1:2497679/1‑149 (MQ=255)
                                                                                                               tttGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAAt                                      >  1:2512228/1‑149 (MQ=255)
                                                                                                                ttGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAAGGGGGCAATCGCAAta                                     >  2:3473848/1‑148 (MQ=255)
                                                                                                                  gAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGTCATTTACAACCCAATTATCTGGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTa                                   >  1:4034709/1‑149 (MQ=255)
                                                                                                                   aaaTGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTAt                                  >  1:1881596/1‑149 (MQ=255)
                                                                                                                     aTGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCg                                >  2:949928/1‑149 (MQ=255)
                                                                                                                      tGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGa                               >  1:3138552/1‑149 (MQ=255)
                                                                                                                      tGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGa                               >  1:3137488/1‑149 (MQ=255)
                                                                                                                          gATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGAtttt                           <  1:3612761/149‑1 (MQ=255)
                                                                                                                           aTGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTg                          >  1:409193/1‑149 (MQ=255)
                                                                                                                                  tGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCAc                   <  1:931913/149‑1 (MQ=255)
                                                                                                                                  tGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCAc                   <  1:1251671/149‑1 (MQ=255)
                                                                                                                                  tGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCAc                   <  2:1151699/149‑1 (MQ=255)
                                                                                                                                  tGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCAc                   <  2:812000/149‑1 (MQ=255)
                                                                                                                                  tGATCAGCAGCGTCGATGATAACGCAGGCCGGGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCAc                   <  2:3864914/149‑1 (MQ=255)
                                                                                                                                   gATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACt                  <  1:2825238/149‑1 (MQ=255)
                                                                                                                                          agcGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACCAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGc           >  2:1003931/1‑149 (MQ=255)
                                                                                                                                           gcGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGcc          <  2:3166725/149‑1 (MQ=255)
                                                                                                                                            cGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCg         <  2:293362/149‑1 (MQ=255)
                                                                                                                                            cGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCg         <  1:1279925/149‑1 (MQ=255)
                                                                                                                                             gTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCgg        <  2:3733014/149‑1 (MQ=255)
                                                                                                                                              tCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGa       <  1:3372774/149‑1 (MQ=255)
                                                                                                                                               cGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGAt      <  2:3818673/149‑1 (MQ=255)
                                                                                                                                                gatgatAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATa     >  2:1670488/1‑149 (MQ=255)
                                                                                                                                                 atgatAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATaa    >  2:2380610/1‑149 (MQ=255)
                                                                                                                                                   gatAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAgg  >  2:2530379/1‑149 (MQ=255)
                                                                                                                                                   gatAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAgg  <  1:2153210/149‑1 (MQ=255)
                                                                                                                                                   gatAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAgg  >  2:152590/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGG  >  CP009273/1129825‑1130120

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