Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,890,960 T→A 100% F32Y (TTC→TAC)  cbrB → PRK09823 family inner membrane protein, creBC regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,890,9600TA98.7% 830.8 / ‑0.5 255F32Y (TTC→TAC) cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base A (134/117);  total (135/120)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 256 total aligned reads displayed.

TTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGT  >  CP009273/3890814‑3891105
                                                                                                                                                  |                                                                                                                                                 
ttCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACg                                                                                                                                                 <  1:3319958/149‑1 (MQ=255)
 tCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGc                                                                                                                                                <  2:1966146/149‑1 (MQ=255)
   aCTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaa                                                                                                                                              >  1:1513346/1‑149 (MQ=255)
     ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaaaa                                                                                                                                            <  1:3150790/149‑1 (MQ=255)
     ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaaaa                                                                                                                                            <  1:1191395/149‑1 (MQ=255)
     ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaaaa                                                                                                                                            >  1:1658746/1‑149 (MQ=255)
       aCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGa                                                                                                                                          <  2:2557809/149‑1 (MQ=255)
       aCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGa                                                                                                                                          >  2:1240295/1‑149 (MQ=255)
        ccccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGaa                                                                                                                                         <  1:3728665/149‑1 (MQ=255)
          ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAAcc                                                                                                                                       <  2:2387150/149‑1 (MQ=255)
                 aaTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGtt                                                                                                                                <  1:2783843/149‑1 (MQ=255)
                  aTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTc                                                                                                                               <  2:3268456/149‑1 (MQ=255)
                     ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCttt                                                                                                                            >  1:2056888/1‑149 (MQ=255)
                     ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCttt                                                                                                                            >  1:1386287/1‑149 (MQ=255)
                       cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTgg                                                                                                                          <  2:1416572/149‑1 (MQ=255)
                        gcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTggg                                                                                                                         >  2:6208/1‑149 (MQ=255)
                           cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGa                                                                                                                      <  1:2705079/149‑1 (MQ=255)
                           cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGa                                                                                                                      <  1:453738/149‑1 (MQ=255)
                            aaaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGat                                                                                                                     <  2:1626408/149‑1 (MQ=255)
                              aGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGataa                                                                                                                   >  1:2668277/1‑149 (MQ=255)
                               ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGataat                                                                                                                  >  1:1176313/1‑149 (MQ=255)
                                 taagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATAc                                                                                                                >  1:1852800/1‑149 (MQ=255)
                                  aagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACa                                                                                                               >  2:823601/1‑149 (MQ=255)
                                    gtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                             >  1:2382824/1‑149 (MQ=255)
                                    gtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                             >  1:494775/1‑149 (MQ=255)
                                     taaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATcc                                                                                                            <  2:1613248/149‑1 (MQ=255)
                                      aaaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAAtacataca                                                                                                           >  2:1855427/1‑149 (MQ=255)
                                      aaaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCa                                                                                                           <  2:946137/149‑1 (MQ=255)
                                       aaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGGGATAATACATCCaa                                                                                                          >  2:2251206/1‑149 (MQ=255)
                                       aaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTGAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCaa                                                                                                          <  2:366619/149‑1 (MQ=255)
                                         gggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAATAACCGCTGGTTCTTTGGGTGATAATACATCCAAtt                                                                                                        >  2:1786325/1‑149 (MQ=255)
                                         gggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAtt                                                                                                        >  1:2330459/1‑149 (MQ=255)
                                          ggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAttt                                                                                                       >  1:1563778/1‑149 (MQ=255)
                                           gAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAtttt                                                                                                      >  1:1993260/1‑149 (MQ=255)
                                           gAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAtttt                                                                                                      <  2:3299877/149‑1 (MQ=255)
                                           gAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAtttt                                                                                                      >  1:3983057/1‑149 (MQ=255)
                                             gattGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTGATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTtc                                                                                                    <  2:3717663/147‑1 (MQ=255)
                                                 gATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCtt                                                                                                <  2:1529756/149‑1 (MQ=255)
                                                  atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                               <  1:2923706/149‑1 (MQ=255)
                                                  atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTATCTTa                                                                                               >  2:2569232/1‑149 (MQ=255)
                                                   tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAt                                                                                              >  1:808782/1‑149 (MQ=255)
                                                    atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAtt                                                                                             >  2:2691645/1‑149 (MQ=255)
                                                    atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAtt                                                                                             <  2:883354/149‑1 (MQ=255)
                                                      gTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGt                                                                                           <  1:368839/149‑1 (MQ=255)
                                                          gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGAt                                                                                       >  1:987935/1‑149 (MQ=255)
                                                          gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGAt                                                                                       <  2:2576419/149‑1 (MQ=255)
                                                          gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGAt                                                                                       >  1:3223236/1‑149 (MQ=255)
                                                           tttCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATa                                                                                      >  1:563027/1‑149 (MQ=255)
                                                             tCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAAc                                                                                    <  2:2066224/149‑1 (MQ=255)
                                                              cACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACt                                                                                   >  2:4080555/1‑149 (MQ=255)
                                                                cGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTAc                                                                                 <  1:2534100/149‑1 (MQ=255)
                                                                  tCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTTTCGATAACTACgg                                                                               >  2:3801642/1‑149 (MQ=255)
                                                                  tCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                               >  2:1930319/1‑149 (MQ=255)
                                                                    gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGa                                                                             >  1:3332697/1‑149 (MQ=255)
                                                                    gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGa                                                                             <  2:135909/149‑1 (MQ=255)
                                                                    gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACAGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGTGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGa                                                                             <  2:3564107/149‑1 (MQ=255)
                                                                         aTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTATTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCTATAACTACGGGAGCTAt                                                                        <  1:3455251/149‑1 (MQ=255)
                                                                         aTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTAt                                                                        >  1:172249/1‑149 (MQ=255)
                                                                         aTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTAt                                                                        >  2:3087466/1‑149 (MQ=255)
                                                                          tACATCACGGACTTTATTTTGCAGCTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTAtt                                                                       >  2:1618457/1‑149 (MQ=255)
                                                                              tCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTg                                                                   >  2:2890463/1‑149 (MQ=255)
                                                                                aCGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCg                                                                 >  2:3104495/1‑149 (MQ=255)
                                                                                 cGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                                <  2:1774784/149‑1 (MQ=255)
                                                                                 cGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                                <  2:2425173/149‑1 (MQ=255)
                                                                                    aCTTTATTTTGCAGTTTTAGGACCGTTAATTGTTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgt                                                             >  2:3046501/1‑149 (MQ=255)
                                                                                      tttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                           >  1:3881247/1‑149 (MQ=255)
                                                                                        tATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAc                                                         <  1:647989/149‑1 (MQ=255)
                                                                                           tttGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCgg                                                      >  2:3883488/1‑149 (MQ=255)
                                                                                            ttGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgt                                                     <  1:1211715/149‑1 (MQ=255)
                                                                                              gCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtgt                                                   >  1:3073234/1‑149 (MQ=255)
                                                                                              gCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtgt                                                   >  2:4058965/1‑149 (MQ=255)
                                                                                                aGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTaa                                                 <  1:879070/149‑1 (MQ=255)
                                                                                                 gTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAAt                                                <  2:1176313/149‑1 (MQ=255)
                                                                                                    ttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGt                                             <  2:1658746/149‑1 (MQ=255)
                                                                                                    ttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGt                                             <  2:797740/149‑1 (MQ=255)
                                                                                                     tAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGtt                                            <  2:1261840/149‑1 (MQ=255)
                                                                                                       ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                          >  1:1217980/1‑149 (MQ=255)
                                                                                                        gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                         <  1:301013/149‑1 (MQ=255)
                                                                                                            gTTAATTGGTGTTCTGGTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                     <  2:2946460/149‑1 (MQ=255)
                                                                                                               aaTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGc                                  <  2:2217169/149‑1 (MQ=255)
                                                                                                                aTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                 >  1:589612/1‑149 (MQ=255)
                                                                                                                 ttGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCg                                >  1:1156071/1‑149 (MQ=255)
                                                                                                                        tCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGa                         >  2:3647509/1‑149 (MQ=255)
                                                                                                                         cTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGAt                        <  1:474810/149‑1 (MQ=255)
                                                                                                                           gTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGTTTGCCTGTCTGCCGGAAAAGATCg                      >  1:2233510/1‑149 (MQ=255)
                                                                                                                              tCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGt                   >  2:3579722/1‑149 (MQ=255)
                                                                                                                                   tCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAg              >  2:2037398/1‑149 (MQ=255)
                                                                                                                                    ccTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGa             >  1:1913506/1‑149 (MQ=255)
                                                                                                                                    ccTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGa             >  1:3561292/1‑149 (MQ=255)
                                                                                                                                    actctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATAAATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGa             >  1:2931967/2‑149 (MQ=255)
                                                                                                                                     ctctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa            <  2:1100911/149‑1 (MQ=255)
                                                                                                                                      tctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa           >  1:2640360/1‑149 (MQ=255)
                                                                                                                                      tctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa           <  1:3630444/149‑1 (MQ=255)
                                                                                                                                          cATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGtt       >  1:596150/1‑149 (MQ=255)
                                                                                                                                          cATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGtt       <  1:3830174/149‑1 (MQ=255)
                                                                                                                                          cATATTCTACGCAAAAGAACCGCTGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa      <  1:222123/149‑1 (MQ=255)
                                                                                                                                           atatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa      >  1:621140/1‑149 (MQ=255)
                                                                                                                                            tatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTAt     <  1:2842468/149‑1 (MQ=255)
                                                                                                                                            tatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTAt     >  1:2015839/1‑149 (MQ=255)
                                                                                                                                               tctacGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGt  >  2:2481365/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGT  >  CP009273/3890814‑3891105

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: