Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,358,920 G→T 100% R147L (CGC→CTC)  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,9200GT98.4% 583.5 / 0.6 185R147L (CGC→CTC) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (77/105);  total (80/105)
Fisher's exact test for biased strand distribution p-value = 7.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 186 total aligned reads displayed.

CTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGC  >  CP009273/4358775‑4359064
                                                                                                                                                 |                                                                                                                                                
cTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGctcg                                                                                                                                                <  2:1082879/148‑1 (MQ=255)
      ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                          >  1:863709/1‑148 (MQ=255)
       tGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCt                                                                                                                                         >  1:78797/1‑148 (MQ=255)
       tGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCt                                                                                                                                         >  2:2097457/1‑148 (MQ=255)
        gccgccGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCtt                                                                                                                                        <  2:1448195/148‑1 (MQ=255)
         ccgccgGTGCACAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATTATGAGCTCGACCGCCTTg                                                                                                                                       <  2:1805454/148‑1 (MQ=255)
          cgccgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGa                                                                                                                                      <  1:1206413/148‑1 (MQ=255)
          cgccgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGa                                                                                                                                      <  1:348772/148‑1 (MQ=255)
                    aaaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGAc                                                                                                                            >  1:1231815/1‑148 (MQ=255)
                     aaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                           <  2:3537559/148‑1 (MQ=255)
                     aaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                           >  2:2990593/1‑148 (MQ=255)
                     aaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                           >  2:2990541/1‑148 (MQ=255)
                      aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGa                                                                                                                          <  2:2003016/148‑1 (MQ=255)
                       gCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTATCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGATTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGAt                                                                                                                         >  1:817879/1‑148 (MQ=255)
                       gCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGAt                                                                                                                         <  1:36210/148‑1 (MQ=255)
                         aTCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCg                                                                                                                       <  1:3118251/148‑1 (MQ=255)
                          tCTGACAGTGAATTTGATGCCTCATTTGCTCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGTGACGGGTGTCGGTAATACTTTTGATGTTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGc                                                                                                                      <  1:3011818/148‑1 (MQ=255)
                            tGACTGTGAAGTTGATGCCGCATTTGCGGTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCa                                                                                                                    <  2:2127823/148‑1 (MQ=255)
                              aCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAgg                                                                                                                  >  1:1839327/1‑148 (MQ=255)
                               cAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTTTGTACCAGCGAATGGATGATGATCTCGACCGCCTTGATCATAAAGACGATCGCCAGGa                                                                                                                 <  2:1796921/148‑1 (MQ=255)
                               cAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGa                                                                                                                 <  2:1797452/148‑1 (MQ=255)
                                  tGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTa                                                                                                              <  1:340293/148‑1 (MQ=255)
                                  tGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTa                                                                                                              <  2:1220834/148‑1 (MQ=255)
                                  tGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTa                                                                                                              >  2:412429/1‑148 (MQ=255)
                                   gAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTaa                                                                                                             <  2:316621/148‑1 (MQ=255)
                                    aaGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAAt                                                                                                            >  1:2919297/1‑148 (MQ=255)
                                          aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                                                                      <  1:3215815/148‑1 (MQ=255)
                                          aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                                                                      >  2:1186145/1‑148 (MQ=255)
                                          aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                                                                      <  2:2030718/148‑1 (MQ=255)
                                          aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                                                                      <  1:2820030/148‑1 (MQ=255)
                                            gCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAAc                                                                                                    >  2:1382522/1‑148 (MQ=255)
                                                 aTTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCgg                                                                                               <  2:3269878/148‑1 (MQ=255)
                                                    tGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAtt                                                                                            >  1:3361803/1‑148 (MQ=255)
                                                       gcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtg                                                                                         <  1:190162/148‑1 (MQ=255)
                                                       gcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtg                                                                                         <  1:189551/148‑1 (MQ=255)
                                                        cTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtgt                                                                                        >  1:1500714/1‑148 (MQ=255)
                                                         ttttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGTTGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtgtg                                                                                       >  1:1540076/1‑148 (MQ=255)
                                                         ttttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtgtg                                                                                       >  1:3127349/1‑148 (MQ=255)
                                                         ttttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtgtg                                                                                       <  2:2416990/148‑1 (MQ=255)
                                                           ttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGTATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGt                                                                                     <  1:256637/148‑1 (MQ=255)
                                                           ttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGt                                                                                     <  1:2056476/148‑1 (MQ=255)
                                                           ttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGt                                                                                     >  2:708230/1‑148 (MQ=255)
                                                            tctcGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGtt                                                                                    <  2:3498632/148‑1 (MQ=255)
                                                                gCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGc                                                                                <  2:1286980/148‑1 (MQ=255)
                                                                     cgtGGGGGCGGTTTTAGCGGCGGGACGGGTGGGGGTTATTCTTTTGTTGGTGTGTACCAGGGTTTGGATGTTGAGCTCGTCGGGGTTGTTCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAttt                                                                           <  2:926005/146‑1 (MQ=255)
                                                                     ccTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGTTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAttt                                                                           <  2:2709090/148‑1 (MQ=255)
                                                                        ggcggcGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTcc                                                                        <  1:879425/148‑1 (MQ=255)
                                                                             cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                   <  2:2169501/148‑1 (MQ=255)
                                                                             cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                   <  2:2169562/148‑1 (MQ=255)
                                                                               gTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttt                                                                 <  1:1695893/148‑1 (MQ=255)
                                                                                 tttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttttt                                                               >  2:2976437/1‑148 (MQ=255)
                                                                                  ttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTc                                                              <  1:1010738/148‑1 (MQ=255)
                                                                                   tAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCg                                                             >  1:1014337/1‑148 (MQ=255)
                                                                                    aGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGt                                                            <  2:1877480/148‑1 (MQ=255)
                                                                                      cgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGttt                                                          <  1:1186145/148‑1 (MQ=255)
                                                                                      cgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGttt                                                          <  2:2919297/148‑1 (MQ=255)
                                                                                      cgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGttt                                                          >  1:1669568/1‑148 (MQ=255)
                                                                                           ggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgca                                                     >  1:1412929/1‑148 (MQ=255)
                                                                                           ggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGgca                                                     >  1:1107970/1‑148 (MQ=255)
                                                                                           ggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTAGTTTTGgca                                                     >  1:315856/1‑148 (MQ=255)
                                                                                               cGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAACCGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCaa                                                 <  1:819977/148‑1 (MQ=255)
                                                                                               cGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCaa                                                 >  1:3450118/1‑148 (MQ=255)
                                                                                                   ggCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGc                                             <  2:3455975/148‑1 (MQ=255)
                                                                                                    gggttAATAGTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAAGGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGcc                                            <  1:2690919/144‑1 (MQ=255)
                                                                                                    gCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGcc                                            >  1:1281382/1‑148 (MQ=255)
                                                                                                    gCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGcc                                            >  2:1874399/1‑148 (MQ=255)
                                                                                                     cGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCa                                           <  1:1335508/148‑1 (MQ=255)
                                                                                                     cGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCa                                           >  2:1473733/1‑148 (MQ=255)
                                                                                                        tAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATaa                                        <  2:2491506/148‑1 (MQ=255)
                                                                                                             cTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCg                                   >  2:3196082/1‑148 (MQ=255)
                                                                                                              ttttGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGc                                  <  1:464224/148‑1 (MQ=255)
                                                                                                                   aTGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCa                             >  1:2417273/1‑148 (MQ=255)
                                                                                                                       tGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGAt                         <  1:2456871/148‑1 (MQ=255)
                                                                                                                         aGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCa                       >  2:556147/1‑148 (MQ=255)
                                                                                                                         aGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCa                       >  1:486074/1‑148 (MQ=255)
                                                                                                                         aGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCa                       >  2:1294749/1‑148 (MQ=255)
                                                                                                                          gTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAg                      >  2:1616012/1‑148 (MQ=255)
                                                                                                                          gTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAg                      >  2:1337092/1‑148 (MQ=255)
                                                                                                                            aCCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAgcg                    <  1:2424191/148‑1 (MQ=255)
                                                                                                                             ccAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAgcgc                   <  1:2172363/148‑1 (MQ=255)
                                                                                                                             ccAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAgcgc                   >  2:872135/1‑148 (MQ=255)
                                                                                                                             ccAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAgcgc                   <  1:1345632/148‑1 (MQ=255)
                                                                                                                             ccAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGGTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCCgcgc                   >  1:1648546/1‑148 (MQ=255)
                                                                                                                              cAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCa                  <  1:914248/148‑1 (MQ=255)
                                                                                                                                gCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATa                >  2:876602/1‑148 (MQ=255)
                                                                                                                                 cGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCATGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATag               <  2:267337/148‑1 (MQ=255)
                                                                                                                                 cGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATAAGCGCATag               >  1:1775893/1‑148 (MQ=255)
                                                                                                                                  gAAAGGATGATGCGCTCCACCGCATTGAACACAAAGACGCTCGCCAGTATTAATGTCGAACCGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga              <  1:1596285/148‑1 (MQ=255)
                                                                                                                                  gAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga              <  1:1597092/148‑1 (MQ=255)
                                                                                                                                  gAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga              <  2:3170291/148‑1 (MQ=255)
                                                                                                                                  gAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga              <  1:1597406/148‑1 (MQ=255)
                                                                                                                                  gAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga              <  1:637898/148‑1 (MQ=255)
                                                                                                                                   aaaGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATagag             >  1:238203/1‑148 (MQ=255)
                                                                                                                                     aGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGAc           >  1:156621/1‑148 (MQ=255)
                                                                                                                                     aGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGAc           >  2:1146447/1‑148 (MQ=255)
                                                                                                                                       gatgatGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGGCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtt         <  1:1983854/148‑1 (MQ=255)
                                                                                                                                         tgatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtttt       <  2:3098194/148‑1 (MQ=255)
                                                                                                                                         tgatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtttt       <  2:2474917/148‑1 (MQ=255)
                                                                                                                                          gatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTc      <  2:3416635/148‑1 (MQ=255)
                                                                                                                                              aGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCCTAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGc  >  2:838244/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGC  >  CP009273/4358775‑4359064

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: