Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 997,965 G→T 100% D235Y (GAT→TAT)  elfG → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273997,9650GT98.6% 917.5 / 1.3 283D235Y (GAT→TAT) elfGputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (143/136);  total (145/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 283 total aligned reads displayed.

ATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCA  >  CP009273/997820‑998109
                                                                                                                                                 |                                                                                                                                                
aTTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTAt                                                                                                                                                >  2:1292967/1‑148 (MQ=255)
aTTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACATGGCTCATGTTACGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTAt                                                                                                                                                >  1:3146143/1‑148 (MQ=255)
  tataTGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGt                                                                                                                                              <  2:953814/148‑1 (MQ=255)
  tataTGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGt                                                                                                                                              >  2:650065/1‑148 (MQ=255)
   atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTc                                                                                                                                             <  1:466326/148‑1 (MQ=255)
   atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTc                                                                                                                                             <  2:1633521/148‑1 (MQ=255)
     atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                           >  2:727166/1‑148 (MQ=255)
     atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                           <  1:636586/148‑1 (MQ=255)
     atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                           <  2:223224/148‑1 (MQ=255)
     atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                           <  2:2494121/148‑1 (MQ=255)
      tGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTc                                                                                                                                          >  2:2079075/1‑148 (MQ=255)
       gAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGTTTTAGGTTATGTCTTCg                                                                                                                                         >  2:1590250/1‑148 (MQ=255)
             tACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCa                                                                                                                                   >  2:396241/1‑148 (MQ=255)
                     aCCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgt                                                                                                                           >  1:1454190/1‑148 (MQ=255)
                            ggTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAg                                                                                                                    >  1:2050085/1‑148 (MQ=255)
                               gATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGc                                                                                                                 >  2:432759/1‑148 (MQ=255)
                                aTTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCa                                                                                                                >  1:2350436/1‑148 (MQ=255)
                                 ttCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCAt                                                                                                               >  2:1624318/1‑148 (MQ=255)
                                 ttCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCAt                                                                                                               <  1:382291/148‑1 (MQ=255)
                                  tCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATa                                                                                                              >  1:2226904/1‑148 (MQ=255)
                                   ccTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATaa                                                                                                             <  2:556385/148‑1 (MQ=255)
                                     tGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAAc                                                                                                           >  2:3751157/1‑148 (MQ=255)
                                     tGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAAc                                                                                                           >  2:2798680/1‑148 (MQ=255)
                                     tGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAAc                                                                                                           >  2:3751197/1‑148 (MQ=255)
                                       cACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACct                                                                                                         <  1:378509/148‑1 (MQ=255)
                                        aCGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACctc                                                                                                        <  1:1146043/148‑1 (MQ=255)
                                                 ggTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCa                                                                                               >  1:1070880/1‑148 (MQ=255)
                                                 ggTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCa                                                                                               >  2:8346/1‑148 (MQ=255)
                                                  gTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAg                                                                                              >  1:2278784/1‑148 (MQ=255)
                                                  gTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAg                                                                                              >  1:2571206/1‑148 (MQ=255)
                                                   tACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGa                                                                                             <  2:54770/148‑1 (MQ=255)
                                                      acCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACg                                                                                          >  2:1499424/1‑148 (MQ=255)
                                                       ccGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                         >  1:2635410/1‑148 (MQ=255)
                                                       ccGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                         >  1:546901/1‑148 (MQ=255)
                                                       ccGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                         >  1:1515048/1‑148 (MQ=255)
                                                       ccGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                         >  1:717532/1‑148 (MQ=255)
                                                          gTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAg                                                                                      <  2:2821505/148‑1 (MQ=255)
                                                           ttttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGa                                                                                     >  2:1309806/1‑148 (MQ=255)
                                                           ttttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGa                                                                                     >  1:1278119/1‑148 (MQ=255)
                                                           ttttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGa                                                                                     >  1:1278220/1‑148 (MQ=255)
                                                            tttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGaa                                                                                    <  1:657770/148‑1 (MQ=255)
                                                            tttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGaa                                                                                    >  2:491283/1‑148 (MQ=255)
                                                             ttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAAc                                                                                   <  2:3086216/148‑1 (MQ=255)
                                                             ttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAAc                                                                                   >  1:213042/1‑148 (MQ=255)
                                                              tAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACt                                                                                  <  1:1406099/148‑1 (MQ=255)
                                                              tAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACt                                                                                  <  1:684724/148‑1 (MQ=255)
                                                                   tATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCaa                                                                             <  1:749072/148‑1 (MQ=255)
                                                                      tACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTc                                                                          <  1:3070795/148‑1 (MQ=255)
                                                                      tACCTGTCTGGTACGATCAAGTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTc                                                                          <  2:2287263/148‑1 (MQ=255)
                                                                        ccTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCa                                                                        >  1:586030/1‑148 (MQ=255)
                                                                        ccTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCa                                                                        >  1:349817/1‑148 (MQ=255)
                                                                            tCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCa                                                                    >  1:192293/1‑148 (MQ=255)
                                                                             cTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAg                                                                   <  1:3273931/148‑1 (MQ=255)
                                                                             cTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAg                                                                   <  1:3272242/148‑1 (MQ=255)
                                                                               ggTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgt                                                                 >  2:3703820/1‑148 (MQ=255)
                                                                                gTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgtg                                                                >  1:994082/1‑148 (MQ=255)
                                                                                 tACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgtgt                                                               <  2:2842342/148‑1 (MQ=255)
                                                                                  aCGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTa                                                              <  1:3728507/148‑1 (MQ=255)
                                                                                    gATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAac                                                            <  1:1737925/148‑1 (MQ=255)
                                                                                    gATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAac                                                            <  1:3630593/148‑1 (MQ=255)
                                                                                     aTCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAaca                                                           <  1:1013643/148‑1 (MQ=255)
                                                                                      tCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCa                                                                              <  2:353026/128‑1 (MQ=255)
                                                                                      tCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAacac                                                          >  1:3531043/1‑148 (MQ=255)
                                                                                      tCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAacac                                                          >  1:425056/1‑148 (MQ=255)
                                                                                      tCAATTCACTTGGCTAATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCAAATTTAAGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAacac                                                          >  1:3020476/1‑148 (MQ=255)
                                                                                          ttCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGa                                                      >  1:1632559/1‑148 (MQ=255)
                                                                                          ttCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGAAACACGGGa                                                      >  1:3214468/1‑148 (MQ=255)
                                                                                           tCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGaa                                                     <  1:820543/148‑1 (MQ=255)
                                                                                           tCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGaa                                                     <  2:675876/148‑1 (MQ=255)
                                                                                            cACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGaaa                                                    >  2:2132739/1‑148 (MQ=255)
                                                                                                tGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAcgcg                                                <  2:2146291/148‑1 (MQ=255)
                                                                                                 ggCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCgg                                               <  1:1378418/148‑1 (MQ=255)
                                                                                                  gCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCggg                                              <  1:1621094/148‑1 (MQ=255)
                                                                                                   cTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGa                                             <  2:3574069/148‑1 (MQ=255)
                                                                                                    tCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGAt                                            <  1:396241/148‑1 (MQ=255)
                                                                                                     cATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGAtt                                           >  2:1935329/1‑148 (MQ=255)
                                                                                                      aTGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTa                                          <  1:3625541/148‑1 (MQ=255)
                                                                                                        gTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                        <  1:1700642/148‑1 (MQ=255)
                                                                                                        gTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                        <  2:2935059/148‑1 (MQ=255)
                                                                                                          tCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTa                                      <  1:495937/148‑1 (MQ=255)
                                                                                                           ccGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTaa                                     <  2:1019389/148‑1 (MQ=255)
                                                                                                            cGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAAt                                    <  2:1066681/148‑1 (MQ=255)
                                                                                                              tCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGt                                  <  1:1292967/148‑1 (MQ=255)
                                                                                                              tCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGt                                  <  2:3249384/148‑1 (MQ=255)
                                                                                                               cAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTc                                 >  2:3472944/1‑148 (MQ=255)
                                                                                                                aaTGCCGGAGAGACAATTGAAGTTGATTTAGGTTTTGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATATACCTCTTGGGGGCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGTTTAGTTAATGTCa                                <  1:922812/148‑1 (MQ=255)
                                                                                                                 aTGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCaa                               >  1:2777135/1‑148 (MQ=255)
                                                                                                                 aTGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCaa                               <  2:1448898/148‑1 (MQ=255)
                                                                                                                   gCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAAcc                             <  2:2619054/148‑1 (MQ=255)
                                                                                                                   gCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAAcc                             <  2:2350234/148‑1 (MQ=255)
                                                                                                                       gagagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                         >  2:1941466/1‑148 (MQ=255)
                                                                                                                          agaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCt                      >  2:3627956/1‑148 (MQ=255)
                                                                                                                          agaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCt                      >  2:1737374/1‑148 (MQ=255)
                                                                                                                             cAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGAc                   >  1:1836799/1‑148 (MQ=255)
                                                                                                                                      ttGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACcac          <  1:2817662/148‑1 (MQ=255)
                                                                                                                                      ttGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACcac          >  2:2565088/1‑148 (MQ=255)
                                                                                                                                         aTTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGa       >  1:160390/1‑148 (MQ=255)
                                                                                                                                            tAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAccc    >  2:289065/1‑148 (MQ=255)
                                                                                                                                              ggTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCa  >  1:2650904/1‑148 (MQ=255)
                                                                                                                                              ggTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCa  >  2:3312850/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCA  >  CP009273/997820‑998109

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: