Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,088,850 C→A 94.7% R93S (CGT→AGT)  rhaR → transcriptional activator of rhaSR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,088,8500CA94.7% 529.7 / 19.2 181R93S (CGT→AGT) rhaRtranscriptional activator of rhaSR
Reads supporting (aligned to +/- strand):  ref base C (5/5);  new base A (106/65);  total (111/70)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

CTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCATG  >  CP009273/4088705‑4088994
                                                                                                                                                 |                                                                                                                                                
 tGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTc                                                                                                                                               >  1:2256104/1‑148 (MQ=255)
  ggCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCt                                                                                                                                              <  1:230911/148‑1 (MQ=255)
    cgcgGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGa                                                                                                                                            >  1:647980/1‑148 (MQ=255)
      cgGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAg                                                                                                                                          >  2:2372012/1‑148 (MQ=255)
       ggTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGc                                                                                                                                         >  1:275660/1‑148 (MQ=255)
        gTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCt                                                                                                                                        <  2:1692914/148‑1 (MQ=255)
         tAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTg                                                                                                                                       >  1:1183195/1‑148 (MQ=255)
           aTGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGaa                                                                                                                                     >  1:3114181/1‑148 (MQ=255)
             ggCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATc                                                                                                                                   <  2:2295423/148‑1 (MQ=255)
                 tGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGa                                                                                                                               >  1:2451573/1‑148 (MQ=255)
                 tGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGa                                                                                                                               <  1:2381639/148‑1 (MQ=255)
                   cATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACt                                                                                                                             >  1:736227/1‑148 (MQ=255)
                       tACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCa                                                                                                                         <  2:2855541/148‑1 (MQ=255)
                        aCTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAg                                                                                                                        >  1:391305/1‑148 (MQ=255)
                        aCTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAg                                                                                                                        >  1:391244/1‑148 (MQ=255)
                              cGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCg                                                                                                                  >  1:183686/1‑148 (MQ=255)
                              cGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCg                                                                                                                  >  1:2410868/1‑148 (MQ=255)
                              cGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCg                                                                                                                  >  1:2411101/1‑148 (MQ=255)
                                aTCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGAt                                                                                                                >  1:2501165/1‑148 (MQ=255)
                                       ttATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGa                                                                                                         <  1:1381993/148‑1 (MQ=255)
                                       ttATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGa                                                                                                         >  1:849965/1‑148 (MQ=255)
                                       ttATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGa                                                                                                         <  1:2372012/148‑1 (MQ=255)
                                       ttATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGa                                                                                                         >  2:605827/1‑148 (MQ=255)
                                        tATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGAt                                                                                                        >  2:2197885/1‑148 (MQ=255)
                                        tATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTAGTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGAt                                                                                                        <  1:3165444/148‑1 (MQ=255)
                                          tCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGAttt                                                                                                      >  2:844618/1‑148 (MQ=255)
                                              aTTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACg                                                                                                  <  2:335158/148‑1 (MQ=255)
                                               ttACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGc                                                                                                 <  1:1986710/148‑1 (MQ=255)
                                                   ccGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAgc                                                                                             >  1:86988/1‑148 (MQ=255)
                                                     gTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAgcgc                                                                                           <  2:539910/148‑1 (MQ=255)
                                                       ggCGATCTCTTTTACATTCATGCTGACGATTAACACTACTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAg                                                                                         <  1:2701998/148‑1 (MQ=255)
                                                       ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAg                                                                                         <  1:2702033/148‑1 (MQ=255)
                                                       ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAg                                                                                         >  1:1143525/1‑148 (MQ=255)
                                                       ggCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAg                                                                                         >  2:697881/1‑148 (MQ=255)
                                                        gCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAgg                                                                                        <  2:962699/148‑1 (MQ=255)
                                                         cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                       >  2:1201547/1‑148 (MQ=255)
                                                         cGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAggg                                                                                       >  2:2388091/1‑148 (MQ=255)
                                                           aTCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                     >  2:2862005/1‑148 (MQ=255)
                                                           aTCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                     >  1:2470795/1‑148 (MQ=255)
                                                           aTCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCa                                                                                     >  2:3124426/1‑148 (MQ=255)
                                                             ctctTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAc                                                                                   >  1:1845881/1‑148 (MQ=255)
                                                             ctctTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAc                                                                                   <  2:1542704/148‑1 (MQ=255)
                                                             ctctTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAc                                                                                   >  2:2256268/1‑148 (MQ=255)
                                                             ctctTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAc                                                                                   >  2:594107/1‑148 (MQ=255)
                                                              tctTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAAcc                                                                                  >  2:1401192/1‑148 (MQ=255)
                                                                 tttACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACcaca                                                                               >  1:2548692/1‑148 (MQ=255)
                                                                 tttACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACcaca                                                                               >  2:3215117/1‑148 (MQ=255)
                                                                   tACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACt                                                                             <  2:2138828/148‑1 (MQ=255)
                                                                   tACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACt                                                                             <  1:2460739/148‑1 (MQ=255)
                                                                    aCATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTg                                                                            >  2:1928586/1‑148 (MQ=255)
                                                                      aTTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgc                                                                          >  1:1694687/1‑148 (MQ=255)
                                                                       ttCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcg                                                                         <  2:925706/148‑1 (MQ=255)
                                                                        tCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGgcgc                                                                        <  1:231155/148‑1 (MQ=255)
                                                                         cATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCt                                                                       >  1:2291994/1‑148 (MQ=255)
                                                                            gCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAg                                                                    <  2:452275/148‑1 (MQ=255)
                                                                             cTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAgg                                                                   <  1:2849585/148‑1 (MQ=255)
                                                                              tGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGt                                                                  <  2:275660/148‑1 (MQ=255)
                                                                                aCGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAg                                                                >  2:1387793/1‑148 (MQ=255)
                                                                                aCGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAg                                                                <  1:1953476/148‑1 (MQ=255)
                                                                                    tAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCCGGGGGCGATTCCGGGATTTAACGCCCGCGCAGGGCAACCACACTGGCGCTTCGGTAGCATg                                                            >  2:2884870/1‑148 (MQ=255)
                                                                                      aaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATggg                                                          <  2:1028510/148‑1 (MQ=255)
                                                                                      aaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATggg                                                          >  1:64096/1‑148 (MQ=255)
                                                                                      aaCACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATggg                                                          <  1:2332789/148‑1 (MQ=255)
                                                                                       acacTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATgggg                                                         <  2:1972782/148‑1 (MQ=255)
                                                                                         acTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGAt                                                       <  1:269563/148‑1 (MQ=255)
                                                                                         acTCCGACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGAt                                                       <  1:2834748/148‑1 (MQ=255)
                                                                                              tACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGgcgc                                                  >  2:2758103/1‑148 (MQ=255)
                                                                                              tACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGgcgc                                                  >  2:21993/1‑148 (MQ=255)
                                                                                                      cGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCgg                                          <  1:593563/148‑1 (MQ=255)
                                                                                                      cGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCgg                                          >  2:3198229/1‑148 (MQ=255)
                                                                                                      cGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCgg                                          >  2:1579340/1‑148 (MQ=255)
                                                                                                      cGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCgg                                          >  2:1579297/1‑148 (MQ=255)
                                                                                                        ttAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATATTGACTGGCAGCGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGAAGGCGCAGGCGCGGCa                                        >  1:2065246/1‑148 (MQ=255)
                                                                                                          aaCGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAgg                                      >  1:296696/1‑148 (MQ=255)
                                                                                                           aCGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGt                                     <  2:2920026/148‑1 (MQ=255)
                                                                                                              aTCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTAt                                  <  1:1401192/148‑1 (MQ=255)
                                                                                                               tCTGGTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATc                                 >  2:1175034/1‑148 (MQ=255)
                                                                                                                  ggTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGt                              <  2:2561499/148‑1 (MQ=255)
                                                                                                                  ggTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGt                              >  1:2236469/1‑148 (MQ=255)
                                                                                                                   gTTTTGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTc                             <  2:1078629/148‑1 (MQ=255)
                                                                                                                    ttttGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCa                            >  2:1942731/1‑148 (MQ=255)
                                                                                                                      ttGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGc                          <  2:832779/148‑1 (MQ=255)
                                                                                                                       tGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTTGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCt                         <  1:3075613/148‑1 (MQ=255)
                                                                                                                       tGCAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCt                         <  2:849965/148‑1 (MQ=255)
                                                                                                                         cAGAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTg                       >  2:1363134/1‑148 (MQ=255)
                                                                                                                          aGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGa                      <  2:1869410/148‑1 (MQ=255)
                                                                                                                           gAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAg                     >  2:1243931/1‑148 (MQ=255)
                                                                                                                           gAATATTATTTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAg                     <  1:756796/148‑1 (MQ=255)
                                                                                                                            aaTATTATTTATTGCCCGGAGAGTCTGAAGCCGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGc                    >  2:1818772/1‑148 (MQ=255)
                                                                                                                               attattTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATg                 >  2:638126/1‑148 (MQ=255)
                                                                                                                               attattTATTGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATg                 <  2:2548692/148‑1 (MQ=255)
                                                                                                                                 tatttattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGaa               >  1:2992584/1‑148 (MQ=255)
                                                                                                                                   tttattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAag             <  2:3024354/148‑1 (MQ=255)
                                                                                                                                     tattGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagta           >  2:2773669/1‑148 (MQ=255)
                                                                                                                                      attGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtag          >  1:662732/1‑148 (MQ=255)
                                                                                                                                       ttGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtagt         <  2:736227/148‑1 (MQ=255)
                                                                                                                                       ttGCCCGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAagtagt         >  1:274957/1‑148 (MQ=255)
                                                                                                                                          cccGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAg      <  2:647980/148‑1 (MQ=255)
                                                                                                                                            cGGAGAGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCa    >  1:2338549/1‑148 (MQ=255)
                                                                                                                                              gagaGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCAtg  <  2:875226/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGGCGCGGTAATGGCCTGCATGTACTCAACGATCGCCCTTATCGCATTACCCGTGGCGATCTCTTTTACATTCATGCTGACGATAAACACTCCTACGCTTCCGTTAACGATCTGGTTTTGCAGAATATTATTTATTGCCCGGAGCGTCTGAAGCTGAATCTTGACTGGCAGGGGGCGATTCCGGGATTTAACGCCAGCGCAGGGCAACCACACTGGCGCTTAGGTAGCATGGGGATGGCGCAGGCGCGGCAGGTTATCGGTCAGCTTGAGCATGAAAGTAGTCAGCATG  >  CP009273/4088705‑4088994

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: