Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,358,943 G→T 5.6% D155Y (GAT→TAT)  fxsA → suppressor of F exclusion of phage T7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,358,9430GT 5.6% 940.7 / 19.3 246D155Y (GAT→TAT) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base G (83/149);  new base T (6/8);  total (89/157)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 251 total aligned reads displayed.

AAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGAACGCATA  >  CP009273/4358795‑4359091
                                                                                                                                                    |                                                                                                                                                    
aaaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                                                      >  1:2439009/1‑149 (MQ=255)
 aaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGa                                                                                                                                                     >  2:3298783/1‑149 (MQ=255)
  aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGAt                                                                                                                                                    <  1:2910522/149‑1 (MQ=255)
   gCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATc                                                                                                                                                   >  1:2231454/1‑149 (MQ=255)
    cATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCg                                                                                                                                                  <  2:312523/149‑1 (MQ=255)
    cATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCAGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCg                                                                                                                                                  <  2:4557773/149‑1 (MQ=255)
      tCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGcc                                                                                                                                                <  2:1028779/149‑1 (MQ=255)
      tCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGcc                                                                                                                                                <  1:2860660/149‑1 (MQ=255)
       cTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCa                                                                                                                                               <  1:1308208/149‑1 (MQ=255)
        tGACAGTGATGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGATGATCGCCAg                                                                                                                                              <  1:4418265/149‑1 (MQ=255)
        tGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAg                                                                                                                                              >  1:38408/1‑149 (MQ=255)
         gACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAgg                                                                                                                                             <  1:996571/149‑1 (MQ=255)
          aCAGTGAAGTTGCTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGa                                                                                                                                            <  2:4582869/149‑1 (MQ=255)
          aCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCTACCGCCTTGATCATAAAGACGATCGCCAGGa                                                                                                                                            >  2:1968290/1‑149 (MQ=255)
          aCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGa                                                                                                                                            >  2:3566051/1‑149 (MQ=255)
            aGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACTATCGCCAGGAtt                                                                                                                                          <  2:2260617/149‑1 (MQ=255)
               gAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAAt                                                                                                                                       >  1:696812/1‑149 (MQ=255)
                 aGTTGCTGCCGCTTTTGCGCTTTTCTCGCTTGCGTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAATGACGATCGCCAGGATTAATGt                                                                                                                                     <  1:9517/149‑1 (MQ=255)
                      aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaaa                                                                                                                                >  1:3866232/1‑149 (MQ=255)
                      aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaaa                                                                                                                                >  1:1476159/1‑149 (MQ=255)
                      aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaaa                                                                                                                                <  2:1382780/149‑1 (MQ=255)
                      aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGATAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaaa                                                                                                                                <  2:3671613/149‑1 (MQ=255)
                         ccGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGc                                                                                                                             >  2:797839/1‑149 (MQ=255)
                          cGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGcc                                                                                                                            <  2:564121/149‑1 (MQ=255)
                          cGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGcc                                                                                                                            >  2:754748/1‑149 (MQ=255)
                          cGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGcc                                                                                                                            <  2:3240751/149‑1 (MQ=255)
                          cGCATTTGCGCTTTTCTCGCAGGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGcc                                                                                                                            <  2:2956397/149‑1 (MQ=255)
                             aTTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGa                                                                                                                         <  1:3064978/149‑1 (MQ=255)
                             aTTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATAGCCAGGATTAATGTCGAAACGCCGGa                                                                                                                         >  1:2805733/1‑149 (MQ=255)
                              tttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAt                                                                                                                        <  2:4767151/149‑1 (MQ=255)
                              tttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAt                                                                                                                        <  2:3400913/149‑1 (MQ=255)
                               ttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAtt                                                                                                                       <  1:784660/149‑1 (MQ=255)
                               tggcgcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATAATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAtt                                                                                                                       <  1:864473/147‑1 (MQ=255)
                                tGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTa                                                                                                                      <  2:652448/149‑1 (MQ=255)
                                 gcgcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAt                                                                                                                     <  2:2709645/149‑1 (MQ=255)
                                    cTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtgtg                                                                                                                  >  1:4758455/1‑149 (MQ=255)
                                     ttttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTgg                                                                                                                 >  2:3000694/1‑149 (MQ=255)
                                      tttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGCATTATGTGGt                                                                                                                >  1:3562559/1‑149 (MQ=255)
                                        tctcGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTa                                                                                                              >  1:4746197/1‑149 (MQ=255)
                                            gCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGcc                                                                                                          <  2:3942799/149‑1 (MQ=255)
                                              aTGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACTATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAt                                                                                                        <  2:2003767/149‑1 (MQ=255)
                                              aTGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAt                                                                                                        <  2:1550927/149‑1 (MQ=255)
                                               tGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAtt                                                                                                       <  2:1784849/149‑1 (MQ=255)
                                                   tGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTcc                                                                                                   <  1:2006978/149‑1 (MQ=255)
                                                        gcggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                                              <  2:2971499/149‑1 (MQ=255)
                                                        gcggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                                              <  1:4656790/149‑1 (MQ=255)
                                                         cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCTACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGtt                                                                                             <  1:1968463/149‑1 (MQ=255)
                                                          ggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttt                                                                                            <  2:4050645/149‑1 (MQ=255)
                                                           gTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGtttt                                                                                           <  2:2648295/149‑1 (MQ=255)
                                                           gTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGtttt                                                                                           >  1:3142498/1‑149 (MQ=255)
                                                            ttttCGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttttt                                                                                          <  1:4080367/149‑1 (MQ=255)
                                                            ttttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttttt                                                                                          <  2:1804408/149‑1 (MQ=255)
                                                             tttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTc                                                                                         <  1:397864/149‑1 (MQ=255)
                                                                aGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAACGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGtt                                                                                      <  1:562322/149‑1 (MQ=255)
                                                                 gcgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGttt                                                                                     <  2:4348779/149‑1 (MQ=255)
                                                                 gcgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGttt                                                                                     <  1:4253691/149‑1 (MQ=255)
                                                                     cGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgc                                                                                 >  2:2081779/1‑149 (MQ=255)
                                                                      gggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgca                                                                                >  1:2502250/1‑149 (MQ=255)
                                                                       ggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgcag                                                                               >  2:241353/1‑149 (MQ=255)
                                                                          ccGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCaa                                                                            >  2:4501496/1‑149 (MQ=255)
                                                                              tGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGc                                                                        <  2:3864613/149‑1 (MQ=255)
                                                                              tGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGc                                                                        >  1:1727109/1‑149 (MQ=255)
                                                                                  ggTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATa                                                                    >  2:1277556/1‑149 (MQ=255)
                                                                                   gTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATaa                                                                   >  2:2525069/1‑149 (MQ=255)
                                                                                               aTGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCAt                                                       <  1:2739931/149‑1 (MQ=255)
                                                                                                 ggTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGa                                                     >  2:372032/1‑149 (MQ=255)
                                                                                                   tGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACTATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATc                                                   <  2:1480848/149‑1 (MQ=255)
                                                                                                          cAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCAt                                            <  1:1256488/149‑1 (MQ=255)
                                                                                                          cAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGAGCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCAt                                            <  1:3378940/149‑1 (MQ=255)
                                                                                                           aGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATa                                           >  1:3184706/1‑149 (MQ=255)
                                                                                                            gCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACTATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATag                                          >  1:2874245/1‑149 (MQ=255)
                                                                                                            gCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATag                                          >  1:1085108/1‑149 (MQ=255)
                                                                                                             cGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga                                         <  1:2699956/149‑1 (MQ=255)
                                                                                                             cGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga                                         <  1:2630436/149‑1 (MQ=255)
                                                                                                               aaaGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAgaga                                       >  1:4302705/1‑149 (MQ=255)
                                                                                                                aaGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGAc                                      <  1:1204076/149‑1 (MQ=255)
                                                                                                                aaGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGAc                                      <  1:4587685/149‑1 (MQ=255)
                                                                                                                   gatgatGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACttt                                   >  1:3677716/1‑149 (MQ=255)
                                                                                                                      gatgaGCTCGACCGCCTTGATCATAAGGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTcc                                <  2:4670267/149‑1 (MQ=255)
                                                                                                                         gaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGATACTTTTCCAGc                             <  1:2927968/149‑1 (MQ=255)
                                                                                                                         gaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGc                             <  1:2600869/149‑1 (MQ=255)
                                                                                                                         gaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGc                             <  1:2345987/149‑1 (MQ=255)
                                                                                                                            ctcgACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGa                          <  2:2805787/147‑1 (MQ=255)
                                                                                                                              cgACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGAcc                        <  1:4558763/149‑1 (MQ=255)
                                                                                                                              cgACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGAcc                        <  1:3005570/149‑1 (MQ=255)
                                                                                                                              cgACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGAcc                        <  2:1691650/149‑1 (MQ=255)
                                                                                                                               gACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCa                       >  2:1934863/1‑149 (MQ=255)
                                                                                                                                 ccGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCAtt                     >  2:3158553/1‑149 (MQ=255)
                                                                                                                                  cGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTTGCAGCAACAGCCATAAACCCGCCAGCATGATCAGCGCATAGAGACTTTTACAGCCGACCAATg                    >  1:3628295/1‑149 (MQ=255)
                                                                                                                                   gCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGc                   >  1:2465867/1‑149 (MQ=255)
                                                                                                                                   gCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGc                   >  2:522437/1‑149 (MQ=255)
                                                                                                                                         aTCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTaa             >  1:281682/1‑149 (MQ=255)
                                                                                                                                          tCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAAc            <  1:1160555/149‑1 (MQ=255)
                                                                                                                                           cATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACa           >  2:1685585/1‑149 (MQ=255)
                                                                                                                                           cATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACa           <  1:342744/149‑1 (MQ=255)
                                                                                                                                             tAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGa         >  2:474257/1‑149 (MQ=255)
                                                                                                                                             tAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTACGGCGTTTTTCGTTTTGGCAGCAAAAGCCATAACCACGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGa         >  1:2591470/1‑149 (MQ=255)
                                                                                                                                               aaGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGAAc       <  2:2178817/149‑1 (MQ=255)
                                                                                                                                                  aCGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGAACGCa    >  1:4387793/1‑149 (MQ=255)
                                                                                                                                                    gATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGAACGCATa  >  1:1334190/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGAACGCATA  >  CP009273/4358795‑4359091

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: