Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,984,987 C→A 100% T38T (ACG→ACT)  yecA ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,984,9870CA97.9% 789.5 / 5.9 250T38T (ACG→ACT) yecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (122/120);  total (125/125)
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: Frequency below/above cutoff threshold.

Only 100 of 255 total aligned reads displayed.

CGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAG  >  CP009273/1984842‑1985131
                                                                                                                                                 |                                                                                                                                                
   tATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCaa                                                                                                                                             >  1:1877946/1‑148 (MQ=255)
    aTCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAAt                                                                                                                                            >  2:904418/1‑148 (MQ=255)
      cGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACTCCGCAGTCAATaa                                                                                                                                          <  1:1356544/148‑1 (MQ=255)
      cGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATaa                                                                                                                                          >  1:333496/1‑148 (MQ=255)
       ggCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATaaa                                                                                                                                         <  2:336290/148‑1 (MQ=255)
           atatGTTGAAAAGCCATATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCATCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGt                                                                                                                                     <  2:3306963/148‑1 (MQ=255)
            tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTc                                                                                                                                    >  1:3562209/1‑148 (MQ=255)
            tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTc                                                                                                                                    <  2:1615259/148‑1 (MQ=255)
            tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTc                                                                                                                                    <  2:1727370/148‑1 (MQ=255)
               gTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAg                                                                                                                                 >  2:3509152/1‑148 (MQ=255)
               gTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTGTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAg                                                                                                                                 <  2:1961639/148‑1 (MQ=255)
                   aaaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTcc                                                                                                                             >  2:2462694/1‑148 (MQ=255)
                    aaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGCCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAAAAAACCGTCCAGCTCCg                                                                                                                            >  1:781638/1‑148 (MQ=255)
                     aaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGc                                                                                                                           >  2:1003482/1‑148 (MQ=255)
                      aGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGcc                                                                                                                          >  1:2200572/1‑148 (MQ=255)
                          aGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACAt                                                                                                                      <  2:2241570/148‑1 (MQ=255)
                           gATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCACTCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATc                                                                                                                     >  1:1698003/1‑148 (MQ=255)
                              tCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCa                                                                                                                    >  2:2160161/1‑146 (MQ=255)
                              tCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCa                                                                                                                    <  1:2160161/146‑1 (MQ=255)
                                aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCTGTCAATAAACCGTCCAGCTCCGCCACATCAAGGa                                                                                                                <  2:2472236/148‑1 (MQ=255)
                                aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGa                                                                                                                >  1:607495/1‑148 (MQ=255)
                                aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGa                                                                                                                <  2:1368408/148‑1 (MQ=255)
                                 tAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                               >  2:631978/1‑148 (MQ=255)
                                  aaaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATg                                                                                                              <  1:1461863/148‑1 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                            >  2:69848/1‑148 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                            >  2:2971145/1‑148 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                            <  1:1364612/148‑1 (MQ=255)
                                    aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                            >  2:924093/1‑148 (MQ=255)
                                     gcgcGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCg                                                                                                           >  2:3070357/1‑148 (MQ=255)
                                     gcgcGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCg                                                                                                           >  1:2088755/1‑148 (MQ=255)
                                          tCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTc                                                                                                      <  1:2325422/148‑1 (MQ=255)
                                           cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                     <  2:3566665/148‑1 (MQ=255)
                                           cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                     >  2:2599371/1‑148 (MQ=255)
                                           cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                     >  2:2599165/1‑148 (MQ=255)
                                           cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                     >  1:702213/1‑148 (MQ=255)
                                           cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                     >  1:3064137/1‑148 (MQ=255)
                                           cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                     >  1:2031756/1‑148 (MQ=255)
                                            aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAg                                                                                                    <  1:3554944/148‑1 (MQ=255)
                                            aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAg                                                                                                    <  1:2231150/148‑1 (MQ=255)
                                             tctcTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                   >  1:3370641/1‑148 (MQ=255)
                                             tctcTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                   <  2:548262/148‑1 (MQ=255)
                                               tctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtgt                                                                                                 >  1:564507/1‑148 (MQ=255)
                                               tctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtgt                                                                                                 <  2:1455099/148‑1 (MQ=255)
                                               tctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGATCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtgt                                                                                                 <  1:2966993/148‑1 (MQ=255)
                                                ctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgtt                                                                                                >  1:3702959/1‑148 (MQ=255)
                                                  ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                              <  2:1573593/148‑1 (MQ=255)
                                                  ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                              >  1:1126517/1‑148 (MQ=255)
                                                  ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                              >  1:2414436/1‑148 (MQ=255)
                                                      ctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAttt                                                                                          >  2:2138299/1‑148 (MQ=255)
                                                         aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGt                                                                                       >  1:2461315/1‑148 (MQ=255)
                                                         aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGt                                                                                       >  1:2553978/1‑148 (MQ=255)
                                                         aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGt                                                                                       >  2:446654/1‑148 (MQ=255)
                                                         aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGt                                                                                       >  2:1420844/1‑148 (MQ=255)
                                                                gcgcGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAata                                                                                >  1:491503/1‑148 (MQ=255)
                                                                gcgcGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAata                                                                                >  1:2144022/1‑148 (MQ=255)
                                                                 cgcgGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAatat                                                                               <  1:2865135/148‑1 (MQ=255)
                                                                    ggCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGt                                                                            <  1:1377510/148‑1 (MQ=255)
                                                                      cacaTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTcc                                                                          <  2:1758084/148‑1 (MQ=255)
                                                                              cAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGACTATCGTCCAGCCATTc                                                                  >  2:1859940/1‑148 (MQ=255)
                                                                              cAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTc                                                                  <  1:3572419/148‑1 (MQ=255)
                                                                              cAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTc                                                                  <  2:874999/148‑1 (MQ=255)
                                                                                gCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCa                                                                <  2:2793765/148‑1 (MQ=255)
                                                                                gCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCa                                                                <  2:879341/148‑1 (MQ=255)
                                                                                gCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCa                                                                <  1:607922/148‑1 (MQ=255)
                                                                                 ccccaccccaCACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCaa                                                               <  2:1561291/148‑1 (MQ=255)
                                                                                       cccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACt                                                         <  2:1129219/148‑1 (MQ=255)
                                                                                         cacaCGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACttt                                                       <  2:3596146/148‑1 (MQ=255)
                                                                                         cacaCGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACttt                                                       >  2:3731718/1‑148 (MQ=255)
                                                                                          acacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTc                                                      >  2:2695213/1‑148 (MQ=255)
                                                                                           cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCg                                                     <  1:631978/148‑1 (MQ=255)
                                                                                             cGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGtt                                                   <  1:3691373/148‑1 (MQ=255)
                                                                                             cGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGtt                                                   <  1:1800877/148‑1 (MQ=255)
                                                                                                cAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTac                                                <  1:69446/148‑2 (MQ=255)
                                                                                                      gCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTcc                                          <  2:2748857/148‑1 (MQ=255)
                                                                                                        cACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCg                                        >  2:2021760/1‑148 (MQ=255)
                                                                                                         aCTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGt                                       <  2:1427980/148‑1 (MQ=255)
                                                                                                          cTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGtt                                      >  1:3060601/1‑148 (MQ=255)
                                                                                                           tGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGttt                                     <  2:1809038/148‑1 (MQ=255)
                                                                                                              tCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCa                                  >  1:1187215/1‑148 (MQ=255)
                                                                                                                cGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATa                                >  2:3515421/1‑148 (MQ=255)
                                                                                                                   ttCAAACTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGc                             >  1:2525149/1‑148 (MQ=255)
                                                                                                                    tCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCt                            <  2:2116951/148‑1 (MQ=255)
                                                                                                                    tCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCt                            >  2:3153124/1‑148 (MQ=255)
                                                                                                                     cAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtt                           <  2:955322/148‑1 (MQ=255)
                                                                                                                       aTCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtttt                         >  2:1689145/1‑148 (MQ=255)
                                                                                                                        tctcTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTc                        >  2:793554/1‑148 (MQ=255)
                                                                                                                        tctcTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTc                        >  2:2156250/1‑148 (MQ=255)
                                                                                                                          tctTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCt                      <  1:1003482/148‑1 (MQ=255)
                                                                                                                          tctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCt                      <  1:1576279/148‑1 (MQ=255)
                                                                                                                          tctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCt                      >  1:1444324/1‑148 (MQ=255)
                                                                                                                           ctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCtt                     >  1:3157474/1‑148 (MQ=255)
                                                                                                                            ttGCGGAGAACTCAACACCGCAGTCTATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtg                    <  1:2775144/148‑1 (MQ=255)
                                                                                                                                  agaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtc              >  1:2730777/1‑148 (MQ=255)
                                                                                                                                   gaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtct             >  2:1926297/1‑148 (MQ=255)
                                                                                                                                     aCTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCa           <  2:3063274/148‑1 (MQ=255)
                                                                                                                                          acacCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCg      >  1:3628741/1‑148 (MQ=255)
                                                                                                                                           cacCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGGTTTCATACGCTTTTCCTTGTGTATCTCAACTCGc     >  1:1788067/1‑148 (MQ=255)
                                                                                                                                            acCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGcc    >  1:450570/1‑148 (MQ=255)
                                                                                                                                              cgcagTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAg  >  1:976172/1‑148 (MQ=255)
                                                                                                                                              cgcagTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAg  <  2:717857/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAG  >  CP009273/1984842‑1985131

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: ‑