Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,585,646 C→A 9.8% R8L (CGC→CTC)  yjjM ← yjjMN operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,585,6460CA 9.8% 678.2 / 31.1 216R8L (CGC→CTC) yjjMyjjMN operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (85/108);  new base A (6/16);  total (91/125)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 217 total aligned reads displayed.

CGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGCAACAAAAAGT  >  CP009273/4585503‑4585793
                                                                                                                                               |                                                                                                                                                   
 gcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                                                                                                                                                >  2:1835912/1‑148 (MQ=255)
 gcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                                                                                                                                                >  1:1776276/1‑148 (MQ=255)
   aTTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAtt                                                                                                                                              <  2:1874269/148‑1 (MQ=255)
    ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAttt                                                                                                                                             >  1:2793711/1‑148 (MQ=255)
    ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAttt                                                                                                                                             >  1:1601723/1‑148 (MQ=255)
     tCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAtttt                                                                                                                                            <  2:2160955/148‑1 (MQ=255)
      cgcgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTtg                                                                                                                                           <  2:1344953/148‑1 (MQ=255)
      cgcgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTtg                                                                                                                                           <  2:2285842/148‑1 (MQ=255)
          tAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGaa                                                                                                                                       >  1:389196/1‑148 (MQ=255)
           aaGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAAc                                                                                                                                      >  1:1871980/1‑148 (MQ=255)
             gTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGa                                                                                                                                    >  2:2894699/1‑148 (MQ=255)
             gTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGa                                                                                                                                    >  2:1007948/1‑148 (MQ=255)
              tGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGAc                                                                                                                                   >  1:1656105/1‑148 (MQ=255)
                  tGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCat                                                                                                                               >  2:1503521/1‑148 (MQ=255)
                   gagaATATGACGCACTGTGGTGCGGGTTTTGTTGTACATTGCCGCCAGCGCACTTTGCGTTGGAAGTGGGGATGGTATTTGACCCCGCGCCATATCATCAATGACCTGGTTAATCAATTTGGGGCTTAAATTTTGTGCACGACTCata                                                                                                                              <  1:309963/148‑1 (MQ=255)
                    agaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGTTGGATGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCatat                                                                                                                             <  1:2038123/148‑1 (MQ=255)
                    agaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCatat                                                                                                                             >  1:919026/1‑148 (MQ=255)
                      aaTATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATAttt                                                                                                                           >  2:49995/1‑148 (MQ=255)
                          tGACGCAGTGTGGTGGGGCTGTTGTTGTACTTTTCCGCCAGCGCACTTTGCGTTGGATGTGGGTAAGGTATATGACCCCGCGCCATATCATCAATGACCTGTTTAATCACTTTGTGGCTTAATTTTTGTGAACGACTCATATTTtcct                                                                                                                       <  1:624221/148‑1 (MQ=255)
                            aCGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAGTGTGGCGTAAATTTTGTGAACGACTCATATTTtcctcc                                                                                                                     >  1:2323577/1‑148 (MQ=255)
                               cACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTtg                                                                                                                  <  2:670311/148‑1 (MQ=255)
                                aCTGTGGTGCGGCTGATGTTTTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTtgt                                                                                                                 <  1:252510/148‑1 (MQ=255)
                                aCTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTtgt                                                                                                                 <  1:1921237/148‑1 (MQ=255)
                                 cTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTtgtg                                                                                                                <  2:1367363/148‑1 (MQ=255)
                                    tggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGt                                                                                                             >  2:1341486/1‑148 (MQ=255)
                                       tgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCa                                                                                                          <  2:2644875/148‑1 (MQ=255)
                                          ggCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTa                                                                                                       <  1:1357929/148‑1 (MQ=255)
                                             tGATTTTGTACATTTCCGCCAGAGCACTTTGCGATGGAAGTTGTGAAGGGATATGACGCCGCGCCATATCATCAATGACCTGTTTAATCACTTTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTaaaa                                                                                                    <  2:1804078/148‑1 (MQ=255)
                                                 gttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCg                                                                                                <  2:1545806/148‑1 (MQ=255)
                                                  ttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGa                                                                                               >  1:1979494/1‑148 (MQ=255)
                                                   tgtACATTTCCGCCAGCGCACTTTGCTTTGGATGTGGGTAAGGGATATGACCCCGCGCCATATCATAAATTACCTGGTTAATCACATTGTGGCTTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGAt                                                                                              <  1:1226929/148‑1 (MQ=255)
                                                   tgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGAt                                                                                              <  1:245150/148‑1 (MQ=255)
                                                      aCATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTa                                                                                           >  1:2171749/1‑148 (MQ=255)
                                                         tttCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTaaaa                                                                                        <  2:3042459/148‑1 (MQ=255)
                                                            ccGCCAGCGCACTTTGCTATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACAt                                                                                     <  2:1135176/148‑1 (MQ=255)
                                                             cGCCAGCGCAGTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCAATATCATCAATGACCTGTTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACAtt                                                                                    <  1:1190793/148‑1 (MQ=255)
                                                             cGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACAtt                                                                                    <  1:1190883/148‑1 (MQ=255)
                                                             cGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACAtt                                                                                    <  2:1562439/148‑1 (MQ=255)
                                                              gCCAGCGCACTTTGCGATGGAAGTGGGGAAGGTATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCTTAAATTTTGTGCACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACAttt                                                                                   <  1:1426905/148‑1 (MQ=255)
                                                              gCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATAATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACAttt                                                                                   <  1:3237653/148‑1 (MQ=255)
                                                               ccAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACAtttt                                                                                  >  1:10077/1‑148 (MQ=255)
                                                                 aGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTa                                                                                <  1:3235895/148‑1 (MQ=255)
                                                                  gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTAt                                                                               <  2:2416859/148‑1 (MQ=255)
                                                                       gtttGCGATGGAAGTGGGGATGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACTTTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCgtg                                                                          <  2:1385385/147‑1 (MQ=255)
                                                                       cTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCgtg                                                                          >  2:2481417/1‑148 (MQ=255)
                                                                       cTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCgtg                                                                          >  2:821231/1‑148 (MQ=255)
                                                                          tGCGATGGATGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATGGTGGCGTAAATTTGGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAACACCCGATTTAAAACATTTTTATGCgtgtgt                                                                       <  1:161359/148‑1 (MQ=255)
                                                                           gCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGtt                                                                      <  2:3249263/148‑1 (MQ=255)
                                                                            cGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGttt                                                                     <  1:2757592/148‑1 (MQ=255)
                                                                              aTGGAAGTGGGGAAGGGGTATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGAATTAAAACATTTTTATGCGTGTGTTTTc                                                                   >  1:2895148/1‑148 (MQ=255)
                                                                                     tGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACaaaaa                                                            <  2:261778/148‑1 (MQ=255)
                                                                                        ggtaGGGGTATGACCCCGCGCCTTATAATCAATGACGTGGTTAATCACATTGTGGAGTAAATTTGGTCAACGACTAATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGa                                                         <  2:2895148/145‑1 (MQ=255)
                                                                                        ggAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGa                                                         <  1:1547734/148‑1 (MQ=255)
                                                                                          aaGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGAtt                                                       <  2:1383489/148‑1 (MQ=255)
                                                                                           aGGGATATGGCCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGAttt                                                      <  1:2434432/148‑1 (MQ=255)
                                                                                            ggcatatGACCCCGCGCCAAATCATCAATGACCTGGTTAATAACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGAtttt                                                     <  2:2286336/145‑1 (MQ=255)
                                                                                             ggATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttc                                                    <  2:2678682/148‑1 (MQ=255)
                                                                                             ggATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttc                                                    <  2:900485/148‑1 (MQ=255)
                                                                                               atatGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAATCATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccg                                                  >  2:1366639/1‑148 (MQ=255)
                                                                                               atatGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccg                                                  >  2:47198/1‑148 (MQ=255)
                                                                                               atatGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccg                                                  <  1:1363247/148‑1 (MQ=255)
                                                                                               atatGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccg                                                  <  2:1991781/148‑1 (MQ=255)
                                                                                               atatGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccg                                                  <  2:3098019/148‑1 (MQ=255)
                                                                                                  tGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccgttc                                               >  2:1414385/1‑148 (MQ=255)
                                                                                                    aCCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccgttccg                                             >  2:1712875/1‑148 (MQ=255)
                                                                                                    aCCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccgttccg                                             <  2:2139354/148‑1 (MQ=255)
                                                                                                    aCCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTttccgttccg                                             <  1:3111973/148‑1 (MQ=255)
                                                                                                     ccccGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGc                                            >  2:2032089/1‑148 (MQ=255)
                                                                                                          cgcCATATCATCAATGCCCTGGTTAATCACATGGTGGCGTAAATTTGGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCa                                       <  2:2315702/148‑1 (MQ=255)
                                                                                                          cgcCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCa                                       >  2:3113290/1‑148 (MQ=255)
                                                                                                           gcaatatCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCAt                                      <  1:2640327/145‑1 (MQ=255)
                                                                                                            ccATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCAtt                                     <  1:497177/148‑1 (MQ=255)
                                                                                                              atatCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCAtttt                                   <  1:43582/148‑1 (MQ=255)
                                                                                                                atcatcAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttc                                 >  1:1042071/1‑148 (MQ=255)
                                                                                                                  catcaATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttcat                               <  2:3152644/148‑1 (MQ=255)
                                                                                                                  catcaATGACCTGGTTAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttcat                               >  2:3161496/1‑148 (MQ=255)
                                                                                                                   atcaATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttcatt                              >  2:1294017/1‑148 (MQ=255)
                                                                                                                    tcaATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttcattt                             >  2:2599074/1‑148 (MQ=255)
                                                                                                                      aaTGACCTGGTTAATCACTTTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttcattttc                           <  1:1503926/148‑1 (MQ=255)
                                                                                                                       aTGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttcattttca                          >  1:1260014/1‑148 (MQ=255)
                                                                                                                       aTGACCTGGTTAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATttttcattttca                          <  2:2203864/148‑1 (MQ=255)
                                                                                                                        tGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTtttcattttcat                         <  2:486800/148‑1 (MQ=255)
                                                                                                                           ccTGGTTAATCACATTGTGGCTTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTTATTGa                      <  1:1449996/148‑1 (MQ=255)
                                                                                                                           ccTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGa                      <  1:1102549/148‑1 (MQ=255)
                                                                                                                            cTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGat                     <  1:1662916/148‑1 (MQ=255)
                                                                                                                                ttAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCt                 >  2:882322/1‑148 (MQ=255)
                                                                                                                                ttAATCACATTGTGGAGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCt                 >  1:863839/1‑148 (MQ=255)
                                                                                                                                 tAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTa                >  2:25277/1‑148 (MQ=255)
                                                                                                                                 tAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTa                >  2:24811/1‑148 (MQ=255)
                                                                                                                                  aaTCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTaa               <  1:1058741/148‑1 (MQ=255)
                                                                                                                                   aTCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAg              >  2:812481/1‑148 (MQ=255)
                                                                                                                                    tCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAgg             <  1:2029022/148‑1 (MQ=255)
                                                                                                                                     cacaTTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGc            >  2:367983/1‑148 (MQ=255)
                                                                                                                                      acaTTGTGGAGTAAATTGTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGca           <  2:2916547/148‑1 (MQ=255)
                                                                                                                                          tgtgGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGcaacaa       >  2:413227/1‑148 (MQ=255)
                                                                                                                                          tgtgGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGcaacaa       >  2:2826720/1‑148 (MQ=255)
                                                                                                                                          tgtgGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGcaacaa       >  2:1798456/1‑148 (MQ=255)
                                                                                                                                           gtgGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGCAACaaa      >  2:192569/1‑148 (MQ=255)
                                                                                                                                            tgGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGCAACaaaa     <  2:1777956/148‑1 (MQ=255)
                                                                                                                                               cGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGCAACAAAAAGt  <  1:283375/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATATTTTCCTCCTTGTGAGTCCATTAAAACCCGATTTAAAACATTTTTATGCGTGTGTTTTCACAAAAAAGATTTTCCGTTCCGCTGCCATTTTTCATTTTCATTGATATCTAAGGCAACAAAAAGT  >  CP009273/4585503‑4585793

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