Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,305,116 C→A 100% intergenic (‑103/+9) apbE ← / ← ompC putative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,305,1160CA95.6% 895.1 / 28.2 272intergenic (‑103/+9)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base A (132/128);  total (136/136)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 275 total aligned reads displayed.

AGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGC  >  CP009273/2304968‑2305262
                                                                                                                                                    |                                                                                                                                                  
agagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCa                                                                                                                                                   <  2:2611052/150‑1 (MQ=255)
   gcgcGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATc                                                                                                                                                <  1:3598331/150‑1 (MQ=255)
    cgcgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCg                                                                                                                                               <  1:3969444/150‑1 (MQ=255)
        gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGAt                                                                                                                                           >  2:3728721/1‑150 (MQ=255)
         cagccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGAtt                                                                                                                                          >  2:823676/1‑150 (MQ=255)
           gccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAg                                                                                                                                        >  1:1142020/1‑150 (MQ=255)
           gccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAg                                                                                                                                        <  2:293974/150‑1 (MQ=255)
            ccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGa                                                                                                                                       <  2:2443934/150‑1 (MQ=255)
              agctgTGCCTCGCGGTTAAAGATTATTTCCATACTTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAAc                                                                                                                                     <  1:992216/150‑1 (MQ=255)
               gcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACt                                                                                                                                    <  1:223562/150‑1 (MQ=255)
                cagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTg                                                                                                                                   <  2:1228418/150‑1 (MQ=255)
                 agTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTgg                                                                                                                                  >  2:3789129/1‑150 (MQ=255)
                 agTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTgg                                                                                                                                  <  2:2489681/150‑1 (MQ=255)
                  gTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTggc                                                                                                                                 >  2:1333534/1‑149 (MQ=255)
                    gCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCaa                                                                                                                               >  1:3206611/1‑150 (MQ=255)
                    gCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCaa                                                                                                                               <  1:2968366/150‑1 (MQ=255)
                      cACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAAc                                                                                                                             >  2:2736792/1‑150 (MQ=255)
                      cACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAAc                                                                                                                             >  2:687143/1‑150 (MQ=255)
                         gcgGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAg                                                                                                                          <  2:638740/150‑1 (MQ=255)
                           tggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGTAAACAATGAAAAAAGGGCCCGCGGGCGGTGTGTTCGATATAAATGGAGATTAGATCTGGCAAGCCAGAc                                                                                                                        >  2:3562389/2‑150 (MQ=255)
                           gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAc                                                                                                                        <  1:796271/150‑1 (MQ=255)
                            ggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGAcc                                                                                                                       <  1:1571367/150‑1 (MQ=255)
                              tAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCa                                                                                                                     >  2:748638/1‑150 (MQ=255)
                              tAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCa                                                                                                                     >  1:696685/1‑150 (MQ=255)
                               aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCag                                                                                                                    <  2:2700440/150‑1 (MQ=255)
                                aaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCaga                                                                                                                   >  2:2921348/1‑150 (MQ=255)
                                aaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCaga                                                                                                                   >  1:2167447/1‑150 (MQ=255)
                                 aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCagag                                                                                                                  <  1:1281037/150‑1 (MQ=255)
                                 aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCagag                                                                                                                  <  2:4578271/150‑1 (MQ=255)
                                 aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCagag                                                                                                                  <  1:5014865/150‑1 (MQ=255)
                                 aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCagag                                                                                                                  <  2:1390806/150‑1 (MQ=255)
                                 aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCagag                                                                                                                  <  2:3767717/150‑1 (MQ=255)
                                 aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCagag                                                                                                                  >  2:3994994/1‑150 (MQ=255)
                                     tATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTAc                                                                                                              <  2:878356/150‑1 (MQ=255)
                                         tCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATg                                                                                                          >  1:1466265/1‑150 (MQ=255)
                                         tCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATg                                                                                                          <  1:2894302/150‑1 (MQ=255)
                                          ccATCCGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACTCTGAAAACAATGAAAATAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGt                                                                                                         <  2:238032/150‑1 (MQ=255)
                                          ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGt                                                                                                         >  2:3946226/1‑150 (MQ=255)
                                          ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTCAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGAttt                                                                                                         <  2:1494633/150‑3 (MQ=255)
                                            aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTa                                                                                                       >  2:3296692/1‑150 (MQ=255)
                                             tACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTAt                                                                                                      <  1:323890/150‑1 (MQ=255)
                                              aCGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATc                                                                                                     >  2:4476178/1‑150 (MQ=255)
                                                gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAg                                                                                                   >  2:855595/1‑150 (MQ=255)
                                                gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGAGGTTATCAg                                                                                                   >  1:4767505/1‑150 (MQ=255)
                                                 tGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                                                  <  2:4998166/150‑1 (MQ=255)
                                                  gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtg                                                                                                 <  2:1743180/150‑1 (MQ=255)
                                                      aTCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGa                                                                                             >  1:2723898/1‑150 (MQ=255)
                                                         cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                          >  2:1294284/1‑150 (MQ=255)
                                                         cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                          <  2:4358527/150‑1 (MQ=255)
                                                          tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGcc                                                                                         >  2:309404/1‑150 (MQ=255)
                                                          tCAGGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGcc                                                                                         <  1:2548186/150‑1 (MQ=255)
                                                           cATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCa                                                                                        >  1:4094057/1‑150 (MQ=255)
                                                            aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAg                                                                                       >  2:1998876/1‑150 (MQ=255)
                                                             tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGc                                                                                      >  1:575822/1‑150 (MQ=255)
                                                               cGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                                                    >  1:2678677/1‑150 (MQ=255)
                                                                gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTc                                                                                   >  2:4112691/1‑150 (MQ=255)
                                                                            gagTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTg                                                                       <  1:3630361/150‑1 (MQ=255)
                                                                             agTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                                      >  1:3720177/1‑150 (MQ=255)
                                                                               tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                                                    >  2:3389296/1‑150 (MQ=255)
                                                                               tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                                                    <  1:5113218/150‑1 (MQ=255)
                                                                                 cacCAAAAAACTGAGTTTGTACGCGGACAACAATGACAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                                                  <  2:1050790/150‑1 (MQ=255)
                                                                                  acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtc                                                                 >  2:2072120/1‑150 (MQ=255)
                                                                                   ccAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcg                                                                >  2:5022144/1‑150 (MQ=255)
                                                                                        aaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAACCCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcag                                                           <  1:2384454/150‑1 (MQ=255)
                                                                                            tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAgg                                                       >  1:942233/1‑150 (MQ=255)
                                                                                            tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAgg                                                       <  2:2261617/150‑1 (MQ=255)
                                                                                            tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAgg                                                       <  1:3131238/150‑1 (MQ=255)
                                                                                             gAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACGGGTTGTCGTCCAGCAGGt                                                      >  1:366878/1‑150 (MQ=255)
                                                                                                  tGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAtt                                                 <  2:1243670/150‑1 (MQ=255)
                                                                                                     aCGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTg                                              <  1:2002159/150‑1 (MQ=255)
                                                                                                       gCTGAAAACAATTAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                                            <  1:1508460/150‑1 (MQ=255)
                                                                                                       gCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                                            <  1:3308306/150‑1 (MQ=255)
                                                                                                       gCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                                            >  1:3821078/1‑150 (MQ=255)
                                                                                                       gCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                                            >  1:3493808/1‑150 (MQ=255)
                                                                                                           aaaaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCa                                        >  2:3367230/1‑150 (MQ=255)
                                                                                                            aaaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCaa                                       >  1:1941264/1‑150 (MQ=255)
                                                                                                               cAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGt                                    >  1:384896/1‑150 (MQ=255)
                                                                                                                aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTa                                   <  2:1479191/150‑1 (MQ=255)
                                                                                                                aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTGGATTTTGTAGTCAACGtc                                   >  2:1904930/1‑149 (MQ=255)
                                                                                                                 aTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAg                                  >  1:219825/1‑150 (MQ=255)
                                                                                                                  tGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAgg                                 >  1:2084573/1‑150 (MQ=255)
                                                                                                                    aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg                               <  1:1469464/150‑1 (MQ=255)
                                                                                                                       aaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGAc                            <  1:447331/150‑1 (MQ=255)
                                                                                                                        aaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACa                           >  1:3361894/1‑150 (MQ=255)
                                                                                                                         aGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACAt                          >  2:4895738/1‑150 (MQ=255)
                                                                                                                          gggCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATg                         >  1:494137/1‑150 (MQ=255)
                                                                                                                           ggCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGt                        >  2:2814861/1‑150 (MQ=255)
                                                                                                                            gCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGtt                       <  1:704820/150‑1 (MQ=255)
                                                                                                                               cGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttttt                    >  2:4295557/1‑150 (MQ=255)
                                                                                                                               cGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttttt                    >  1:3419222/1‑150 (MQ=255)
                                                                                                                                 cAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgt                  >  1:4672548/1‑150 (MQ=255)
                                                                                                                                  aGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgtt                 >  2:2825818/1‑150 (MQ=255)
                                                                                                                                   ggCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                >  1:1033577/1‑150 (MQ=255)
                                                                                                                                      ccTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCCGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTTGGTGGACATGTTTTTGTTGacg             >  2:4787210/1‑148 (MQ=255)
                                                                                                                                       cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt            <  2:3747943/150‑1 (MQ=255)
                                                                                                                                       cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt            >  2:3207404/1‑150 (MQ=255)
                                                                                                                                         ttGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag          >  1:892260/1‑150 (MQ=255)
                                                                                                                                             tCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg      >  1:1354794/1‑150 (MQ=255)
                                                                                                                                             tCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg      >  1:391605/1‑150 (MQ=255)
                                                                                                                                               gATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCATGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggta    <  1:339248/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGC  >  CP009273/2304968‑2305262

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: