Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,923,028 C→A 100% 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,0280CA99.4% 1123.2 / ‑2.1 347D476Y (GAC→TAC) ravAhexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (177/166);  total (179/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 347 total aligned reads displayed.

GCCAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGTA  >  CP009273/3922881‑3923174
                                                                                                                                                   |                                                                                                                                                  
gCCAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGc                                                                                                                                                  >  1:1466083/1‑150 (MQ=255)
  cagcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCa                                                                                                                                                <  1:5000305/150‑1 (MQ=255)
  cagcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCa                                                                                                                                                >  2:3701021/1‑150 (MQ=255)
   agcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAg                                                                                                                                               >  2:3441844/1‑150 (MQ=255)
   agcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAg                                                                                                                                               >  1:3211987/1‑150 (MQ=255)
       gTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGAt                                                                                                                                           <  1:998978/150‑1 (MQ=255)
        tGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATTCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATa                                                                                                                                          >  1:775295/1‑150 (MQ=255)
          cGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATaaa                                                                                                                                        <  2:3784976/150‑1 (MQ=255)
          cGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATaaa                                                                                                                                        >  1:3226048/1‑150 (MQ=255)
           gatgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                       >  2:3188954/1‑150 (MQ=255)
           gatgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                       >  2:1213642/1‑150 (MQ=255)
           gatgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                       >  2:5164004/1‑150 (MQ=255)
           gatgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                       >  1:4020659/1‑150 (MQ=255)
            atgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACa                                                                                                                                      <  1:4786805/150‑1 (MQ=255)
             tgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAg                                                                                                                                     <  2:2436701/150‑1 (MQ=255)
               atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTCACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   >  1:276346/1‑150 (MQ=255)
               atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   <  1:3894390/150‑1 (MQ=255)
               atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   >  1:1973804/1‑150 (MQ=255)
               atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   <  1:4517604/150‑1 (MQ=255)
                 catcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAt                                                                                                                                 >  2:2501179/1‑150 (MQ=255)
                 catcTCTTAGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATACTGTAGGCTGGCTTCAATGCGACCTAACCAGTTGCCAGGGATAAACAGACAt                                                                                                                                 <  2:2069055/150‑1 (MQ=255)
                  atcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtt                                                                                                                                >  1:2930024/1‑150 (MQ=255)
                   tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAttt                                                                                                                               <  1:957930/150‑1 (MQ=255)
                        tCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCt                                                                                                                          >  2:2342938/1‑150 (MQ=255)
                         cGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                         >  1:2357322/1‑150 (MQ=255)
                         cGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                         <  1:4126547/150‑1 (MQ=255)
                          gATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGt                                                                                                                        >  2:1604050/1‑150 (MQ=255)
                           aTCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGtt                                                                                                                       <  1:1374101/150‑1 (MQ=255)
                              aaTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGc                                                                                                                    >  2:2053614/1‑150 (MQ=255)
                               aTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCt                                                                                                                   >  1:2044297/1‑150 (MQ=255)
                               aTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCt                                                                                                                   <  1:2265809/150‑1 (MQ=255)
                                    tctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATa                                                                                                              <  2:5300954/150‑1 (MQ=255)
                                      ttCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATaaa                                                                                                            <  1:4315641/150‑1 (MQ=255)
                                           tGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgtg                                                                                                       >  1:3746287/1‑150 (MQ=255)
                                               ggCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgc                                                                                                   >  1:5226340/1‑150 (MQ=255)
                                                 ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 <  2:1677676/150‑1 (MQ=255)
                                                  cAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgcttg                                                                                                <  1:2572432/150‑1 (MQ=255)
                                                  cAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgcttg                                                                                                >  1:1735581/1‑150 (MQ=255)
                                                    gCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTgcg                                                                                              <  1:3015293/150‑1 (MQ=255)
                                                     cATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTgcgc                                                                                             >  1:3553320/1‑150 (MQ=255)
                                                      aTCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTCCGCTGAATAAAGCGTGTTGCTTGCGcc                                                                                            >  1:3299163/1‑150 (MQ=255)
                                                       tCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCa                                                                                           >  1:867407/1‑150 (MQ=255)
                                                         acaTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGt                                                                                         >  2:4692021/1‑150 (MQ=255)
                                                           aTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCg                                                                                       >  1:4931645/1‑150 (MQ=255)
                                                                 cGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCa                                                                                 <  1:3665711/150‑1 (MQ=255)
                                                                 cGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCa                                                                                 >  2:3635439/1‑150 (MQ=255)
                                                                  gTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCaa                                                                                >  2:4574713/1‑150 (MQ=255)
                                                                   tATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAg                                                                               <  1:1547304/150‑1 (MQ=255)
                                                                     tCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGct                                                                             <  2:4223751/150‑1 (MQ=255)
                                                                      ccAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGTACCGACATCCTGTAGGCTGGCTTCAATGCGACATAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGctc                                                                            <  1:246710/150‑1 (MQ=255)
                                                                         gCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTc                                                                         <  2:3392805/150‑1 (MQ=255)
                                                                          cGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTcc                                                                        <  1:1124923/150‑1 (MQ=255)
                                                                          cGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTcc                                                                        <  1:2201985/150‑1 (MQ=255)
                                                                          cGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTcc                                                                        <  1:4706073/150‑1 (MQ=255)
                                                                           gTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTcct                                                                       >  2:4324640/1‑149 (MQ=255)
                                                                           gTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCa                                                                       >  1:538789/1‑150 (MQ=255)
                                                                             tAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgt                                                                     <  2:3097860/150‑1 (MQ=255)
                                                                               gCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttg                                                                   <  1:739471/150‑1 (MQ=255)
                                                                                 aTTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgtt                                                                 <  1:1873702/150‑1 (MQ=255)
                                                                                 aTTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgtt                                                                 >  1:5144522/1‑150 (MQ=255)
                                                                                  ttgttgTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttg                                                                >  2:2910419/1‑150 (MQ=255)
                                                                                   tgttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttgt                                                               >  2:5079376/1‑150 (MQ=255)
                                                                                       gtgCCTGGCGAATCTGTGCACCGACATTCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTa                                                           <  2:775295/150‑1 (MQ=255)
                                                                                        tgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAt                                                          <  1:2882461/150‑1 (MQ=255)
                                                                                         gCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAtt                                                         <  2:3031020/150‑1 (MQ=255)
                                                                                          ccTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAttt                                                        >  2:4988160/1‑150 (MQ=255)
                                                                                           cTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTc                                                       >  2:4810780/1‑150 (MQ=255)
                                                                                                gAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAg                                                  <  2:4055509/150‑1 (MQ=255)
                                                                                                gAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAg                                                  <  2:3325612/150‑1 (MQ=255)
                                                                                                 aaTCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAgg                                                 <  2:2732867/150‑1 (MQ=255)
                                                                                                    cTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAg                                              >  2:1396533/1‑150 (MQ=255)
                                                                                                    cTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAg                                              <  1:506546/150‑1 (MQ=255)
                                                                                                      gtgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGAc                                            >  1:630511/1‑150 (MQ=255)
                                                                                                        gCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCt                                          >  1:1563601/1‑150 (MQ=255)
                                                                                                         cACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCtt                                         <  1:3639835/150‑1 (MQ=255)
                                                                                                         cACCGACATCCGGTAGGCTGGCTTCTTTGCGACCTATCCCGTAGCCTGGGATATACAGACATTTTTGCTGTTCGCTGTATATAGCGTGTTGCTTGCGCCAGTCGCTTTCATGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCtt                                         <  2:2293502/150‑1 (MQ=255)
                                                                                                          aCCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTc                                        >  2:691083/1‑150 (MQ=255)
                                                                                                           ccGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCa                                       <  1:2097879/150‑1 (MQ=255)
                                                                                                                aTCCTGTAGGATGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAAAGACATTTTTGCTGTTCGCGGAATAAAGCGAGTTGCTTGAGACAGTCGCTATCAAGCTCTTCCAGTTGTTGTTTTAATTCAACAGGCAGACCTTCAGCAGa                                  >  2:3128912/1‑150 (MQ=255)
                                                                                                                    tGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGaa                              >  1:4595042/1‑150 (MQ=255)
                                                                                                                      tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                            <  2:3462843/150‑1 (MQ=255)
                                                                                                                        ggctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCgg                          >  1:1299078/1‑150 (MQ=255)
                                                                                                                         gctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCggg                         >  1:1707907/1‑150 (MQ=255)
                                                                                                                              ctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTg                    >  1:2473880/1‑150 (MQ=255)
                                                                                                                                 cAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCa                 <  2:3038277/150‑1 (MQ=255)
                                                                                                                                 cAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCa                 <  2:2164515/150‑1 (MQ=255)
                                                                                                                                  aaTGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAg                >  2:49305/1‑150 (MQ=255)
                                                                                                                                    tGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGACGTGCCAGCg              >  2:458915/1‑150 (MQ=255)
                                                                                                                                    tGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACAATGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCg              >  2:3928275/1‑150 (MQ=255)
                                                                                                                                     gCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGt             <  1:269913/150‑1 (MQ=255)
                                                                                                                                       gACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTAc           <  2:3193558/150‑1 (MQ=255)
                                                                                                                                       gACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTAc           >  1:4336488/1‑150 (MQ=255)
                                                                                                                                       gACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTAGCGCCAGTCGCTTTCAAGCACGTCCAGATGTTGTTTTTTTTCACCAGGCAGACCTTCAGACCAATAACCCCGAAATTCCAGCGAAc           >  1:4075436/1‑150 (MQ=255)
                                                                                                                                         ccTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCTGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTg         <  1:1156139/150‑1 (MQ=255)
                                                                                                                                          cTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGc        >  2:577494/1‑150 (MQ=255)
                                                                                                                                            aaCCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCt      >  1:3319054/1‑150 (MQ=255)
                                                                                                                                             aCCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCtt     <  1:3441845/150‑1 (MQ=255)
                                                                                                                                              ccAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTg    >  2:1363101/1‑150 (MQ=255)
                                                                                                                                               cAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGt   <  2:2803484/150‑1 (MQ=255)
                                                                                                                                                aGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGTa  >  2:1788150/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGTA  >  CP009273/3922881‑3923174

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: