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,9600TA97.7% 689.2 / 3.2 218F32Y (TTC→TAC) cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base A (111/102);  total (112/106)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
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 219 total aligned reads displayed.

TCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCG  >  CP009273/3890815‑3891104
                                                                                                                                                 |                                                                                                                                                
   cTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaa                                                                                                                                             >  2:1021612/1‑148 (MQ=255)
    ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTACAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGTACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCATCTACATATTCTACGCaaa                                                                                                                                            >  1:684408/1‑148 (MQ=255)
    ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaaa                                                                                                                                            <  2:1470726/148‑1 (MQ=255)
    ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaaa                                                                                                                                            <  1:2176931/148‑1 (MQ=255)
    ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaaa                                                                                                                                            >  1:460596/1‑148 (MQ=255)
    ttACCACTCTTAAATTGGAGCTCAAAGGTAATTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTTCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaaa                                                                                                                                            <  2:2977808/148‑1 (MQ=255)
        cccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGaa                                                                                                                                        <  2:1576494/148‑1 (MQ=255)
        cccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGaa                                                                                                                                        <  2:2912181/148‑1 (MQ=255)
         cctcTTAAATTGTCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTTATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAAc                                                                                                                                       <  1:124694/145‑1 (MQ=255)
           tcTTACATGGGCGCTAACAGGTAAGTAAAGGGAGTTTGTTATGTCTGTTTCACGTCGGGGAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACAGATTCTACGCAAAAGAACCg                                                                                                                                     <  2:2539064/147‑1 (MQ=255)
           gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCg                                                                                                                                     <  1:2234650/148‑1 (MQ=255)
                 aTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGtt                                                                                                                               <  2:2438983/148‑1 (MQ=255)
                   tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCt                                                                                                                             <  1:104487/148‑1 (MQ=255)
                   tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCt                                                                                                                             <  1:1903447/148‑1 (MQ=255)
                   tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCt                                                                                                                             <  1:2513341/148‑1 (MQ=255)
                   tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTAAGCAAAAGAACCGCTGGTTCt                                                                                                                             >  1:1579027/1‑148 (MQ=255)
                   tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATAAATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCt                                                                                                                             >  1:1504127/1‑148 (MQ=255)
                    ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCtt                                                                                                                            <  1:510828/148‑1 (MQ=255)
                      cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTg                                                                                                                          >  2:2839338/1‑148 (MQ=255)
                      cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTg                                                                                                                          <  1:913832/148‑1 (MQ=255)
                          cAACGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                      <  2:200050/148‑1 (MQ=255)
                          cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                      >  2:647986/1‑148 (MQ=255)
                           acaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGa                                                                                                                     <  2:2281813/146‑1 (MQ=255)
                            aaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGat                                                                                                                    >  2:1806325/1‑148 (MQ=255)
                             aGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGata                                                                                                                   >  1:1720749/1‑148 (MQ=255)
                             aGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGata                                                                                                                   >  1:351925/1‑148 (MQ=255)
                              ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGataa                                                                                                                  >  2:2664158/1‑148 (MQ=255)
                              ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGataa                                                                                                                  <  2:1186540/148‑1 (MQ=255)
                                 aagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATAc                                                                                                               <  1:1076122/148‑1 (MQ=255)
                                  agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACa                                                                                                              <  2:2203974/148‑1 (MQ=255)
                                    taaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                            >  2:2835694/1‑148 (MQ=255)
                                    taaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                            >  2:2438568/1‑148 (MQ=255)
                                        gggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAt                                                                                                        <  2:1720749/148‑1 (MQ=255)
                                         ggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAtt                                                                                                       >  1:1725541/1‑148 (MQ=255)
                                          gAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAttt                                                                                                      <  1:3215566/148‑1 (MQ=255)
                                          gAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAttt                                                                                                      >  1:2138630/1‑148 (MQ=255)
                                            gTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAttttt                                                                                                    >  1:1535030/1‑148 (MQ=255)
                                              ttGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTtct                                                                                                  >  2:1538628/1‑148 (MQ=255)
                                                gATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTctct                                                                                                <  2:2975161/148‑1 (MQ=255)
                                                 atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCtt                                                                                               <  2:1246511/148‑1 (MQ=255)
                                                  tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                              <  1:195716/148‑1 (MQ=255)
                                                  tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                              <  2:879060/148‑1 (MQ=255)
                                                  tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                              >  2:2766275/1‑148 (MQ=255)
                                                   atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAt                                                                                             <  1:2518061/148‑1 (MQ=255)
                                                   atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAt                                                                                             <  1:3539984/148‑1 (MQ=255)
                                                      tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGt                                                                                          <  2:2249488/148‑1 (MQ=255)
                                                           ttCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCtata                                                                                     >  2:2081492/1‑148 (MQ=255)
                                                           ttCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATa                                                                                     <  1:2073873/148‑1 (MQ=255)
                                                             cACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAAc                                                                                   >  2:372226/1‑148 (MQ=255)
                                                                 tCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                               >  2:2263872/1‑148 (MQ=255)
                                                                     gTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                           <  1:3124737/148‑1 (MQ=255)
                                                                     gTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                           >  1:2578599/1‑148 (MQ=255)
                                                                     gTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                           >  1:798442/1‑148 (MQ=255)
                                                                     gTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                           >  1:2579928/1‑148 (MQ=255)
                                                                      tAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGc                                                                          <  2:3436107/148‑1 (MQ=255)
                                                                        aTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTa                                                                        >  2:235995/1‑148 (MQ=255)
                                                                        aTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTa                                                                        >  2:324174/1‑148 (MQ=255)
                                                                        aTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTa                                                                        <  1:1327444/148‑1 (MQ=255)
                                                                         tACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTAt                                                                       <  1:2259627/148‑1 (MQ=255)
                                                                           cATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTc                                                                     >  2:127497/1‑148 (MQ=255)
                                                                                   aCTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttg                                                             <  1:296803/148‑1 (MQ=255)
                                                                                    cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgt                                                            >  2:347532/1‑148 (MQ=255)
                                                                                    cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgt                                                            >  1:2171421/1‑148 (MQ=255)
                                                                                     tttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                           >  1:407646/1‑148 (MQ=255)
                                                                                      ttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                          <  1:974245/148‑1 (MQ=255)
                                                                                      ttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                          >  2:1890080/1‑148 (MQ=255)
                                                                                      ttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                          >  2:593859/1‑148 (MQ=255)
                                                                                       tATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTaa                                                         >  2:1548070/1‑148 (MQ=255)
                                                                                          tttGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCg                                                      <  1:372226/148‑1 (MQ=255)
                                                                                          tttGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGGTGTTAACCg                                                      >  1:1963957/1‑148 (MQ=255)
                                                                                            tGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgt                                                    <  2:3056253/148‑1 (MQ=255)
                                                                                               aGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTa                                                 <  2:2744858/148‑1 (MQ=255)
                                                                                                      ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTg                                          <  2:346125/148‑1 (MQ=255)
                                                                                                       gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                         <  1:2716769/148‑1 (MQ=255)
                                                                                                        aCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                        >  2:3534677/1‑148 (MQ=255)
                                                                                                         ccGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                       <  1:2796629/148‑1 (MQ=255)
                                                                                                            ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                    >  2:216233/1‑148 (MQ=255)
                                                                                                            ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                    <  2:414767/148‑1 (MQ=255)
                                                                                                            ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                    >  2:1464119/1‑148 (MQ=255)
                                                                                                             tAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCt                                   >  1:1420882/1‑148 (MQ=255)
                                                                                                              aaTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTCCGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTg                                  <  1:1290757/148‑1 (MQ=255)
                                                                                                               aTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGc                                 >  2:2388027/1‑148 (MQ=255)
                                                                                                               aTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGc                                 >  1:3018134/1‑148 (MQ=255)
                                                                                                                 tGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCg                               >  2:607742/1‑148 (MQ=255)
                                                                                                                   gtgtTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGa                             >  1:2380477/1‑148 (MQ=255)
                                                                                                                     gtTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaa                           <  1:560615/148‑1 (MQ=255)
                                                                                                                      ttCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaaa                          >  1:546223/1‑148 (MQ=255)
                                                                                                                       tCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAg                         <  1:3317661/148‑1 (MQ=255)
                                                                                                                           tttCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCg                     >  2:481257/1‑148 (MQ=255)
                                                                                                                           tttCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCg                     >  2:243006/1‑148 (MQ=255)
                                                                                                                           tttCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCg                     <  2:2945941/148‑1 (MQ=255)
                                                                                                                                 gTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTcac               <  1:2664158/148‑1 (MQ=255)
                                                                                                                                  tCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTcaca              >  2:762874/1‑148 (MQ=255)
                                                                                                                                    ctctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGa            >  1:3052338/1‑148 (MQ=255)
                                                                                                                                     tctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa           <  1:764927/148‑1 (MQ=255)
                                                                                                                                      ctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa          >  2:1378392/1‑148 (MQ=255)
                                                                                                                                           tatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa     <  2:351925/148‑1 (MQ=255)
                                                                                                                                           tatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa     <  2:788179/148‑1 (MQ=255)
                                                                                                                                           tatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa     >  2:918586/1‑148 (MQ=255)
                                                                                                                                            atTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTAt    >  1:593560/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCG  >  CP009273/3890815‑3891104

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: