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,0280CA98.4% 1009.6 / 4.2 316D476Y (GAC→TAC) ravAhexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base A (136/174);  total (137/179)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 316 total aligned reads displayed.

CAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGTA  >  CP009273/3922883‑3923174
                                                                                                                                                 |                                                                                                                                                  
 agcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAg                                                                                                                                               <  2:4074186/150‑1 (MQ=255)
 agcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAAGCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAg                                                                                                                                               >  2:4732664/1‑150 (MQ=255)
  gcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGATGACGGCCAGCATCACATTAAGCGTATCCATCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAgg                                                                                                                                              <  2:2446556/150‑1 (MQ=255)
   cagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAggg                                                                                                                                             <  2:4846430/150‑1 (MQ=255)
   cagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAggg                                                                                                                                             >  2:3081624/1‑150 (MQ=255)
    agTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGa                                                                                                                                            <  1:3137044/150‑1 (MQ=255)
    agTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGa                                                                                                                                            >  2:779094/1‑150 (MQ=255)
    agTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGa                                                                                                                                            <  1:3990232/150‑1 (MQ=255)
    agTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGa                                                                                                                                            <  1:3470943/150‑1 (MQ=255)
     gTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGAt                                                                                                                                           <  2:2750209/150‑1 (MQ=255)
      tGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATa                                                                                                                                          <  2:5005583/150‑1 (MQ=255)
      tGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATa                                                                                                                                          <  2:2816737/150‑1 (MQ=255)
      tGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATa                                                                                                                                          <  1:4494783/150‑1 (MQ=255)
          atgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACa                                                                                                                                      >  1:3627236/1‑150 (MQ=255)
           tgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAg                                                                                                                                     <  1:1973582/150‑1 (MQ=255)
           tgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAg                                                                                                                                     <  1:2332988/150‑1 (MQ=255)
           tgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAg                                                                                                                                     >  1:4181971/1‑150 (MQ=255)
             atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   <  1:4288845/150‑1 (MQ=255)
                  ctctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtttt                                                                                                                              <  1:2201047/150‑1 (MQ=255)
                  ctctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtttt                                                                                                                              >  1:3973943/1‑150 (MQ=255)
                  ctctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtttt                                                                                                                              >  2:1130352/1‑150 (MQ=255)
                  ctctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtttt                                                                                                                              <  2:2640041/150‑1 (MQ=255)
                   tcgtCGATCAATCCCTCTTCGCTGACGGCCAGCATCCCATTTCGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGACTCAGTGCACCGACATCCTGGAGGCTGGCTTCAATGCGTCCTAACCAGGAGCCAGGGATAAACAGCCAttttt                                                                                                                             <  2:3328832/147‑1 (MQ=255)
                     ttCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGc                                                                                                                           >  2:199500/1‑150 (MQ=255)
                       cGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGTCTTCAATGCGACCTAACCAGTAGCCATGGATAAACAGACATTTTTGCTg                                                                                                                         <  2:4582786/150‑1 (MQ=255)
                       cGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                         <  2:5113559/150‑1 (MQ=255)
                        gATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGt                                                                                                                        <  1:3705205/150‑1 (MQ=255)
                         aTCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGtt                                                                                                                       <  2:1376620/150‑1 (MQ=255)
                             aTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCt                                                                                                                   <  1:3580262/150‑1 (MQ=255)
                               cccTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGa                                                                                                                 >  2:173650/1‑150 (MQ=255)
                                ccTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGaa                                                                                                                <  2:4476858/150‑1 (MQ=255)
                                  tctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATa                                                                                                              >  2:2309692/1‑150 (MQ=255)
                                      cGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATACACAGACATTTTTGCTGTTCGCTGAATAAAGc                                                                                                          >  1:2634028/1‑150 (MQ=255)
                                        cTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgt                                                                                                        <  2:1617823/150‑1 (MQ=255)
                                        cTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgt                                                                                                        >  2:3533708/1‑150 (MQ=255)
                                        cTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgt                                                                                                        <  2:3933681/150‑1 (MQ=255)
                                        cTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgt                                                                                                        <  1:674945/150‑1 (MQ=255)
                                         tGACGGCCAGCATCACATTCAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgtg                                                                                                       <  1:150615/150‑1 (MQ=255)
                                          gACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgtgt                                                                                                      <  1:2753701/150‑1 (MQ=255)
                                            cGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttg                                                                                                    >  1:4645895/1‑150 (MQ=255)
                                            cGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttg                                                                                                    <  1:5096482/150‑1 (MQ=255)
                                             ggCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgc                                                                                                   <  2:4944404/150‑1 (MQ=255)
                                              gCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgct                                                                                                  >  1:2449556/1‑150 (MQ=255)
                                              gCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgct                                                                                                  >  2:1116019/1‑150 (MQ=255)
                                               ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 >  2:3851635/1‑150 (MQ=255)
                                               ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 <  1:2954354/150‑1 (MQ=255)
                                               ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 <  1:3434295/150‑1 (MQ=255)
                                                        caTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTc                                                                                        <  1:223852/150‑1 (MQ=255)
                                                         aTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCg                                                                                       >  1:1472524/1‑150 (MQ=255)
                                                         aTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCg                                                                                       <  1:958426/150‑1 (MQ=255)
                                                            aaGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCtt                                                                                    >  2:3807361/1‑150 (MQ=255)
                                                            aaGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCtt                                                                                    <  2:3491360/150‑1 (MQ=255)
                                                             aGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCttt                                                                                   <  1:820700/150‑1 (MQ=255)
                                                                    ccAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGctc                                                                            >  1:1308378/1‑150 (MQ=255)
                                                                      aGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCtt                                                                          >  1:4890548/1‑150 (MQ=255)
                                                                        cGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTAGCTGTTCGCTGAATAAATCGTGTTGATTGCGCCAGTCGCTTTCAAGCTCTTcc                                                                        >  2:4641392/1‑150 (MQ=255)
                                                                          ttAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAg                                                                      >  2:3628949/1‑150 (MQ=255)
                                                                                ttgttgTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttg                                                                <  2:1033011/150‑1 (MQ=255)
                                                                                 tgttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttgt                                                               >  1:1272189/1‑150 (MQ=255)
                                                                                   ttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGttt                                                             >  2:1923638/1‑150 (MQ=255)
                                                                                   ttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGttt                                                             <  2:2920123/150‑1 (MQ=255)
                                                                                   ttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGttt                                                             <  1:2370677/150‑1 (MQ=255)
                                                                                    tgtgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGtttt                                                            >  1:4753817/1‑150 (MQ=255)
                                                                                    tgtgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGtttt                                                            >  1:3373556/1‑150 (MQ=255)
                                                                                     gtgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTa                                                           <  1:753348/150‑1 (MQ=255)
                                                                                        ccTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAttt                                                        <  1:173299/150‑1 (MQ=255)
                                                                                        ccTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTGCAAGCTCTTCCAGTTGTTGTTGTAttt                                                        >  2:4393327/1‑150 (MQ=255)
                                                                                         cTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTc                                                       >  2:4786528/1‑150 (MQ=255)
                                                                                              gAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAg                                                  >  2:907504/1‑150 (MQ=255)
                                                                                               aaTCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAgg                                                 >  2:1468835/1‑150 (MQ=255)
                                                                                                aTCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGc                                                <  1:1647505/150‑1 (MQ=255)
                                                                                                   tgtgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGa                                             >  2:4583995/1‑150 (MQ=255)
                                                                                                     tgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGAcc                                           >  2:2242761/1‑150 (MQ=255)
                                                                                                       cACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAGCAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCtt                                         <  1:3906478/150‑1 (MQ=255)
                                                                                                         ccGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCa                                       <  1:2067626/150‑1 (MQ=255)
                                                                                                         ccGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCa                                       <  2:2517961/150‑1 (MQ=255)
                                                                                                          cGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAg                                      <  1:709905/150‑1 (MQ=255)
                                                                                                             cATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCg                                   <  2:3270346/150‑1 (MQ=255)
                                                                                                              aTCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGa                                  <  1:538032/150‑1 (MQ=255)
                                                                                                               tCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGAt                                 <  1:1541744/150‑1 (MQ=255)
                                                                                                                 cTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGa                               <  2:2862920/150‑1 (MQ=255)
                                                                                                                    tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                            <  2:3151372/150‑1 (MQ=255)
                                                                                                                    tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                            >  1:1053522/1‑150 (MQ=255)
                                                                                                                    tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                            >  1:3725985/1‑150 (MQ=255)
                                                                                                                     aGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCg                           >  1:1595270/1‑150 (MQ=255)
                                                                                                                      ggctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTGCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAACgg                          >  2:5148787/1‑150 (MQ=255)
                                                                                                                       gctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCggg                         >  1:5008977/1‑150 (MQ=255)
                                                                                                                        ctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGa                        >  1:3802787/1‑150 (MQ=255)
                                                                                                                        ctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGa                        <  1:4327887/150‑1 (MQ=255)
                                                                                                                         tggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGaa                       <  2:3981449/150‑1 (MQ=255)
                                                                                                                               cAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCa                 >  1:4167422/1‑150 (MQ=255)
                                                                                                                                aaTGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAg                <  1:468189/150‑1 (MQ=255)
                                                                                                                                aaTGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGGTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAACGGGAAGTGCCAg                >  2:963919/1‑150 (MQ=255)
                                                                                                                                   gCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGt             <  1:1877456/150‑1 (MQ=255)
                                                                                                                                     gACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTAc           >  1:4233907/1‑150 (MQ=255)
                                                                                                                                         tAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGcc       >  2:5071175/1‑150 (MQ=255)
                                                                                                                                         tAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGcc       >  2:2140535/1‑150 (MQ=255)
                                                                                                                                          aaCCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCt      >  1:1122755/1‑150 (MQ=255)
                                                                                                                                            ccAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTg    >  2:3581750/1‑150 (MQ=255)
                                                                                                                                            ccAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTg    >  2:2897168/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGTA  >  CP009273/3922883‑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: