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,9200GT99.5% 669.7 / ‑4.6 204R147L (CGC→CTC) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (82/121);  minor base A (1/0);  total (83/121)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-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 204 total aligned reads displayed.

TTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGC  >  CP009273/4358776‑4359064
                                                                                                                                                |                                                                                                                                                
ttttATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGa                                                                                                                                               <  2:705773/148‑1 (MQ=255)
  ttATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGAcc                                                                                                                                             >  1:976644/1‑148 (MQ=255)
    aTTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGc                                                                                                                                           <  2:1485969/148‑1 (MQ=255)
     ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                          <  1:1612803/148‑1 (MQ=255)
     ttGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGcc                                                                                                                                          <  2:2798827/148‑1 (MQ=255)
        ccgccgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTg                                                                                                                                       >  2:638070/1‑148 (MQ=255)
         cgccgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGa                                                                                                                                      <  2:1915490/148‑1 (MQ=255)
          gccgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGAt                                                                                                                                     <  1:676161/148‑1 (MQ=255)
           ccgGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATc                                                                                                                                    <  1:1829879/148‑1 (MQ=255)
              gTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATa                                                                                                                                 <  2:141447/148‑1 (MQ=255)
               tGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATaa                                                                                                                                >  1:254275/1‑148 (MQ=255)
               tGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATaa                                                                                                                                <  1:1645961/148‑1 (MQ=255)
                  aaaaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGa                                                                                                                             >  2:1129826/1‑148 (MQ=255)
                   aaaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGAc                                                                                                                            >  2:3208866/1‑148 (MQ=255)
                    aaGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                           <  2:2573307/148‑1 (MQ=255)
                    aaGCATCTGACAGTGAAGTTGATGCCGCATGTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACg                                                                                                                           <  2:564373/148‑1 (MQ=255)
                     aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGa                                                                                                                          >  1:1937821/1‑148 (MQ=255)
                     aGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGa                                                                                                                          >  2:2081154/1‑148 (MQ=255)
                      gCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGAt                                                                                                                         >  2:1730511/1‑148 (MQ=255)
                      gCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGAt                                                                                                                         <  1:3089499/148‑1 (MQ=255)
                         tCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGc                                                                                                                      <  1:691477/148‑1 (MQ=255)
                          cTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGcc                                                                                                                     >  2:2660350/1‑148 (MQ=255)
                             aCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAgg                                                                                                                  <  1:1154235/148‑1 (MQ=255)
                             aCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAgg                                                                                                                  <  1:3208866/148‑1 (MQ=255)
                             aCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAgg                                                                                                                  >  1:2107070/1‑148 (MQ=255)
                              cAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGa                                                                                                                 <  2:1182888/148‑1 (MQ=255)
                                 tGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTa                                                                                                              <  2:1070904/148‑1 (MQ=255)
                                  gAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCGTGGCGGCGGTTTTTGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGATAGGATGATGAGCTCGACCGCCTTGATCTTAAAGACGATCGCCAGGATTaa                                                                                                             <  2:3444508/148‑1 (MQ=255)
                                   aaGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAAt                                                                                                            <  1:3224396/148‑1 (MQ=255)
                                   aaGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAAt                                                                                                            <  1:298070/148‑1 (MQ=255)
                                    aGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATg                                                                                                           <  2:2391852/148‑1 (MQ=255)
                                    aGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATg                                                                                                           <  1:2835362/148‑1 (MQ=255)
                                     gTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGt                                                                                                          <  2:362930/148‑1 (MQ=255)
                                         aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                                                                      <  2:906427/148‑1 (MQ=255)
                                         aTGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                                                                      <  2:354107/148‑1 (MQ=255)
                                            ccGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                                                                                   >  2:1292134/1‑148 (MQ=255)
                                            ccGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                                                                                   <  1:1088302/148‑1 (MQ=255)
                                             cgatttggcgcTTTTCTCGCATGCGTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGTATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGc                                                                                                  <  1:3266992/141‑1 (MQ=255)
                                             cGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGc                                                                                                  <  2:3301999/148‑1 (MQ=255)
                                                aTTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGTTAATAATTTTGTTGGTGAGTACCAGCGAAAGGATTATGAGCTCGACCGCCTTGATAATAAAGACGATCGCCAGGATTAATGTCGAAACGCCgg                                                                                               <  1:213154/148‑1 (MQ=255)
                                                 tttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGa                                                                                              <  2:1891896/148‑1 (MQ=255)
                                                  ttGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGAt                                                                                             >  2:3352627/1‑148 (MQ=255)
                                                     cgcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAt                                                                                          <  2:91616/148‑1 (MQ=255)
                                                      gcTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtg                                                                                         <  1:1159262/148‑1 (MQ=255)
                                                        ttttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGTTAATACTTTTGATGGTGAGTACCAGCGAAAGTATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtgtg                                                                                       <  1:641193/148‑1 (MQ=255)
                                                          ttCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGt                                                                                     <  2:1449025/148‑1 (MQ=255)
                                                             tcGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTAt                                                                                  <  1:660233/148‑1 (MQ=255)
                                                             tcGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGGCGAAACGCCGGATTATGTGGTTAt                                                                                  <  2:519589/148‑1 (MQ=255)
                                                              cGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATg                                                                                 <  1:258017/148‑1 (MQ=255)
                                                                 aTGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCa                                                                              <  1:847291/148‑1 (MQ=255)
                                                                  tGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAt                                                                             >  2:3421413/1‑148 (MQ=255)
                                                                  tGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAt                                                                             <  2:3220709/148‑1 (MQ=255)
                                                                   gCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATAATTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTTATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAtt                                                                            <  2:2948409/148‑1 (MQ=255)
                                                                     cTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCGGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCAtttt                                                                          <  1:462622/148‑1 (MQ=255)
                                                                      tGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTc                                                                         <  2:1590998/148‑1 (MQ=255)
                                                                       ggcggcGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCACGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGGCGAAACGCCGGATTATGTGGGTATGCCATTTTcc                                                                        >  2:2073997/1‑148 (MQ=255)
                                                                        gcggcgGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCg                                                                       <  2:1363437/148‑1 (MQ=255)
                                                                            cggTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGt                                                                   >  2:34609/1‑148 (MQ=255)
                                                                                tttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTCCCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttttt                                                               <  1:3423846/148‑1 (MQ=255)
                                                                                tttAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGttttt                                                               >  1:2367143/1‑148 (MQ=255)
                                                                                  tAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCg                                                             <  1:2339273/148‑1 (MQ=255)
                                                                                    gcgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGtt                                                           >  2:2316095/1‑148 (MQ=255)
                                                                                    gcgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGtt                                                           <  1:782130/148‑1 (MQ=255)
                                                                                     cgcCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGttt                                                          <  1:638070/148‑1 (MQ=255)
                                                                                          ggACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTGTTTCGTTTTGgca                                                     >  2:2289804/1‑148 (MQ=255)
                                                                                              cGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCaa                                                 <  1:1292134/148‑1 (MQ=255)
                                                                                                   gCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTCTGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCACCAGcc                                            >  1:1874294/1‑148 (MQ=255)
                                                                                                   gCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGcc                                            >  2:1738251/1‑148 (MQ=255)
                                                                                                      gTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATa                                         >  2:1047982/1‑148 (MQ=255)
                                                                                                        aaTACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAAc                                       >  1:1836453/1‑148 (MQ=255)
                                                                                                        aaTACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAAc                                       >  2:474048/1‑148 (MQ=255)
                                                                                                        aaTACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAAc                                       >  1:1233721/1‑148 (MQ=255)
                                                                                                          tACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAAccc                                     >  2:1454933/1‑148 (MQ=255)
                                                                                                           aCTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAAcccc                                    >  1:561834/1‑148 (MQ=255)
                                                                                                             ttttGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGc                                  >  2:2090401/1‑148 (MQ=255)
                                                                                                             ttttGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGc                                  <  2:2864898/148‑1 (MQ=255)
                                                                                                             ttttGATGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGc                                  >  1:3404239/1‑148 (MQ=255)
                                                                                                                  aTGGTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCa                             <  2:2223530/148‑1 (MQ=255)
                                                                                                                    ggTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATg                           <  1:40318/148‑1 (MQ=255)
                                                                                                                     gTGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACACCCATAACACCGCCAGCCTGa                          >  1:1522198/1‑148 (MQ=255)
                                                                                                                      tGAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGAt                         <  1:2435573/148‑1 (MQ=255)
                                                                                                                       gAGTACCAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATc                        <  1:3020956/148‑1 (MQ=255)
                                                                                                                            ccAGCGAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAgcgc                   >  2:2977150/1‑148 (MQ=255)
                                                                                                                                 gAAAGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATaga              <  1:1738737/148‑1 (MQ=255)
                                                                                                                                  aaaGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATagag             <  2:1776017/148‑1 (MQ=255)
                                                                                                                                   aaGGATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAgaga            <  1:2660350/148‑1 (MQ=255)
                                                                                                                                    aggaacatgaGCTCGACCGCCTTGATCATAAAGACGATCGCTAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGAc           >  1:50677/7‑148 (MQ=255)
                                                                                                                                     ggATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCTAAACGCCGGATTATGTGGTTATGCCATTATCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATATAGACt          >  1:1708050/1‑148 (MQ=255)
                                                                                                                                     ggATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACt          >  2:2373898/1‑148 (MQ=255)
                                                                                                                                     ggATGATGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACt          <  2:1807634/148‑1 (MQ=255)
                                                                                                                                      gatgatGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtt         >  2:3150223/1‑148 (MQ=255)
                                                                                                                                      gatgatGAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtt         >  2:1404912/1‑148 (MQ=255)
                                                                                                                                       atgatgAGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACttt        <  2:2342201/148‑1 (MQ=255)
                                                                                                                                        tgatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtttt       >  2:2102079/1‑148 (MQ=255)
                                                                                                                                        tgatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACtttt       >  2:1933055/1‑148 (MQ=255)
                                                                                                                                         gatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTc      <  2:3309406/148‑1 (MQ=255)
                                                                                                                                         gatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTc      <  2:830647/148‑1 (MQ=255)
                                                                                                                                         gatgaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTc      >  2:2040846/1‑148 (MQ=255)
                                                                                                                                            gaGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAg   >  1:1683600/1‑148 (MQ=255)
                                                                                                                                             aGCTCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGc  <  1:1730511/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGC  >  CP009273/4358776‑4359064

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: