Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,923,028 C→A 94.8% D476Y (GAC→TAC)  ravA ← hexameric AAA+ MoxR family ATPase, putative molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,923,0280CA94.8% 561.1 / 21.2 192D476Y (GAC→TAC) ravAhexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (5/5);  new base A (85/97);  total (90/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 192 total aligned reads displayed.

AGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTG  >  CP009273/3922884‑3923172
                                                                                                                                                |                                                                                                                                                
agcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGcc                                                                                                                                               <  1:816724/148‑1 (MQ=255)
agcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGcc                                                                                                                                               <  2:362615/148‑1 (MQ=255)
 gcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCa                                                                                                                                              <  1:1104622/148‑1 (MQ=255)
 gcagTGCGATGATCATCTCTTCGATCAATCCCTCGTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCa                                                                                                                                              <  2:2489011/148‑1 (MQ=255)
  cagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAg                                                                                                                                             >  1:2130932/1‑148 (MQ=255)
     tGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGa                                                                                                                                          >  2:303339/1‑148 (MQ=255)
      gCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGAt                                                                                                                                         >  2:2343176/1‑148 (MQ=255)
       cGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATa                                                                                                                                        <  1:486843/148‑1 (MQ=255)
       cGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATa                                                                                                                                        <  1:1916513/148‑1 (MQ=255)
            atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAg                                                                                                                                   >  2:2804628/1‑148 (MQ=255)
            atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAg                                                                                                                                   <  1:641267/148‑1 (MQ=255)
             tcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGa                                                                                                                                  <  1:193540/148‑1 (MQ=255)
             tcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTAGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGa                                                                                                                                  <  1:2326712/148‑1 (MQ=255)
                tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAt                                                                                                                               >  1:2657418/1‑148 (MQ=255)
                tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAt                                                                                                                               <  2:483502/148‑1 (MQ=255)
                 ctctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtt                                                                                                                              >  2:3111415/1‑148 (MQ=255)
                 ctctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtt                                                                                                                              >  2:150499/1‑148 (MQ=255)
                   ctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtttt                                                                                                                            >  2:1468457/1‑148 (MQ=255)
                        aTCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                       <  2:2227922/148‑1 (MQ=255)
                        aTCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                       <  2:2961825/148‑1 (MQ=255)
                         tCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGt                                                                                                                      >  2:2948561/1‑148 (MQ=255)
                            aTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCg                                                                                                                   <  1:1567022/148‑1 (MQ=255)
                            aTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCg                                                                                                                   <  2:2877632/148‑1 (MQ=255)
                              cccTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCt                                                                                                                 >  1:2565739/1‑148 (MQ=255)
                                 tctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGaa                                                                                                              <  2:312794/148‑1 (MQ=255)
                                 tctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGaa                                                                                                              >  1:1435028/1‑148 (MQ=255)
                                 tcgtCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGaa                                                                                                              <  2:1980662/145‑1 (MQ=255)
                                  ctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAAt                                                                                                             <  1:1973772/148‑1 (MQ=255)
                                       cTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGc                                                                                                        >  2:1075284/1‑148 (MQ=255)
                                          aCGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgtg                                                                                                     >  2:1882849/1‑148 (MQ=255)
                                             gCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttg                                                                                                  >  2:3097183/1‑148 (MQ=255)
                                               cAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgct                                                                                                >  2:1133900/1‑148 (MQ=255)
                                                aGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                               >  2:3137981/1‑148 (MQ=255)
                                                 gCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgcttg                                                                                              <  1:2179302/148‑1 (MQ=255)
                                                         ttAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTc                                                                                      >  2:2680818/1‑148 (MQ=255)
                                                         ttAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTc                                                                                      >  1:513905/1‑148 (MQ=255)
                                                         ttAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTc                                                                                      >  2:68479/1‑148 (MQ=255)
                                                          tAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCg                                                                                     >  2:1879786/1‑148 (MQ=255)
                                                          tAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCg                                                                                     >  2:645592/1‑148 (MQ=255)
                                                            aGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCt                                                                                   >  2:1169372/1‑148 (MQ=255)
                                                             gCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCtt                                                                                  >  2:2327832/1‑148 (MQ=255)
                                                                  tCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAg                                                                             >  1:1212329/1‑148 (MQ=255)
                                                                   ccAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGc                                                                            >  1:72164/1‑148 (MQ=255)
                                                                       cGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCtt                                                                        >  2:1399158/1‑148 (MQ=255)
                                                                       cGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCtt                                                                        >  2:1088322/1‑148 (MQ=255)
                                                                         ttAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTcc                                                                      >  2:2046898/1‑148 (MQ=255)
                                                                          tAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCa                                                                     >  2:2286927/1‑148 (MQ=255)
                                                                          tAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCa                                                                     >  2:625153/1‑148 (MQ=255)
                                                                          tAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCa                                                                     >  1:2428081/1‑148 (MQ=255)
                                                                                tgttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgtt                                                               >  2:1704240/1‑148 (MQ=255)
                                                                                tgttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgtt                                                               >  1:1903410/1‑148 (MQ=255)
                                                                                tgttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgtt                                                               <  1:1312972/148‑1 (MQ=255)
                                                                                   tgtgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttgtt                                                            >  1:622729/1‑148 (MQ=255)
                                                                                    gtgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGttt                                                           <  2:2162266/148‑1 (MQ=255)
                                                                                     tgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGtttt                                                          <  2:986411/148‑1 (MQ=255)
                                                                                      gCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTa                                                         <  1:856729/148‑1 (MQ=255)
                                                                                       ccTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAt                                                        <  2:3028550/148‑1 (MQ=255)
                                                                                       ccTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAt                                                        <  1:1802070/148‑1 (MQ=255)
                                                                                          ggCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTc                                                     >  1:2482668/1‑148 (MQ=255)
                                                                                           gCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCa                                                    <  2:1666218/148‑1 (MQ=255)
                                                                                           gCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCa                                                    >  2:2119731/1‑148 (MQ=255)
                                                                                            cGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCAc                                                   <  2:411569/148‑1 (MQ=255)
                                                                                               aTCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAg                                                >  1:1940727/1‑148 (MQ=255)
                                                                                                tCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAgg                                               >  1:1708597/1‑148 (MQ=255)
                                                                                                   gtgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAg                                            <  1:2330638/148‑1 (MQ=255)
                                                                                                   gtgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAg                                            <  1:1379474/148‑1 (MQ=255)
                                                                                                    tgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGa                                           <  2:134843/148‑1 (MQ=255)
                                                                                                     gCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGAc                                          <  2:2251461/148‑1 (MQ=255)
                                                                                                     gCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGAc                                          <  2:883845/148‑1 (MQ=255)
                                                                                                      cACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGAcc                                         <  1:221870/148‑1 (MQ=255)
                                                                                                        ccGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTAAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCtt                                       >  1:1611667/1‑148 (MQ=255)
                                                                                                          gACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACAAGGCAGACCTTCa                                     >  2:2415568/1‑148 (MQ=255)
                                                                                                              tCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCg                                 <  1:1992180/148‑1 (MQ=255)
                                                                                                              tCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCg                                 >  2:1132424/1‑148 (MQ=255)
                                                                                                              tCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCg                                 <  1:2809062/148‑1 (MQ=255)
                                                                                                                cTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGAt                               <  2:436910/148‑1 (MQ=255)
                                                                                                                   tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGaa                            >  2:1912641/1‑148 (MQ=255)
                                                                                                                    aGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAc                           <  1:800654/148‑1 (MQ=255)
                                                                                                                     ggctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                          >  2:1188822/1‑148 (MQ=255)
                                                                                                                           ctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAg                    >  1:300799/1‑148 (MQ=255)
                                                                                                                           ctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAg                    >  2:1167996/1‑148 (MQ=255)
                                                                                                                             tCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTg                  <  1:1520698/148‑1 (MQ=255)
                                                                                                                             tCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTg                  >  2:2912895/1‑148 (MQ=255)
                                                                                                                             tCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTg                  >  2:101560/1‑148 (MQ=255)
                                                                                                                              cAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGc                 <  1:1024992/148‑1 (MQ=255)
                                                                                                                              cAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGc                 <  2:1323013/148‑1 (MQ=255)
                                                                                                                                 tGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAg              >  1:2327593/1‑148 (MQ=255)
                                                                                                                                 tGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAg              >  1:1036412/1‑148 (MQ=255)
                                                                                                                                    gACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGGTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGt           >  1:185726/1‑148 (MQ=255)
                                                                                                                                     aCCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTa          >  1:1434268/1‑148 (MQ=255)
                                                                                                                                     aCCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTa          <  2:2130932/148‑1 (MQ=255)
                                                                                                                                     aCCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTa          <  1:1087650/148‑1 (MQ=255)
                                                                                                                                      ccTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTAc         <  2:2492265/148‑1 (MQ=255)
                                                                                                                                        tAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTAAAGCTCTTACAGTTGTTTTTTTTTTTAACCAAGCCAACCTTAAACCGATGTAACGGGAATTGCCcgcgaaaat       >  1:850652/1‑143 (MQ=255)
                                                                                                                                           ccAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCt    <  1:3155038/148‑1 (MQ=255)
                                                                                                                                           ccAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCt    <  1:3267073/148‑1 (MQ=255)
                                                                                                                                           ccAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCt    >  2:773377/1‑148 (MQ=255)
                                                                                                                                           ccAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCt    >  2:1855340/1‑148 (MQ=255)
                                                                                                                                             aGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTg  >  2:259482/1‑148 (MQ=255)
                                                                                                                                             aGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTg  >  2:1630452/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTG  >  CP009273/3922884‑3923172

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: