Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,923,028 C→A 92.0% 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,0280CA92.0% 901.9 / 70.9 314D476Y (GAC→TAC) ravAhexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (6/19);  new base A (126/163);  total (132/182)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 314 total aligned reads displayed.

GCCAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGTA  >  CP009273/3922881‑3923174
                                                                                                                                                   |                                                                                                                                                  
  cagcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCa                                                                                                                                                <  1:2743660/150‑1 (MQ=255)
  cagcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCa                                                                                                                                                <  1:877599/150‑1 (MQ=255)
  cagcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCa                                                                                                                                                <  1:1472735/150‑1 (MQ=255)
  cagcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCa                                                                                                                                                <  1:684788/150‑1 (MQ=255)
   agcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAg                                                                                                                                               >  1:4248278/1‑150 (MQ=255)
    gcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAgg                                                                                                                                              <  2:3042434/150‑1 (MQ=255)
      agTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGAGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGa                                                                                                                                            >  2:3341490/1‑150 (MQ=255)
           gatgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAAc                                                                                                                                       <  1:1198775/150‑1 (MQ=255)
           gatgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                       >  2:4300727/1‑150 (MQ=255)
           gatgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                       <  2:695354/150‑1 (MQ=255)
              gatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGa                                                                                                                                    <  2:4150766/150‑1 (MQ=255)
              gatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGa                                                                                                                                    <  2:4112221/150‑1 (MQ=255)
               atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   >  2:1574164/1‑150 (MQ=255)
               atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   <  1:1599270/150‑1 (MQ=255)
               atcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGAc                                                                                                                                   <  2:292417/150‑1 (MQ=255)
                  atcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtt                                                                                                                                <  1:2427278/150‑1 (MQ=255)
                   tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAttt                                                                                                                               <  1:3144866/150‑1 (MQ=255)
                   tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAttt                                                                                                                               <  2:1901344/150‑1 (MQ=255)
                   tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAttt                                                                                                                               >  1:560823/1‑150 (MQ=255)
                    ctctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAtttt                                                                                                                              <  1:3243665/150‑1 (MQ=255)
                     tctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACAttttt                                                                                                                             <  2:360611/150‑1 (MQ=255)
                       ttCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGc                                                                                                                           >  1:4153715/1‑150 (MQ=255)
                       ttCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGc                                                                                                                           <  1:191898/150‑1 (MQ=255)
                        tCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCt                                                                                                                          >  2:90741/1‑150 (MQ=255)
                        tCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCt                                                                                                                          <  1:2781109/150‑1 (MQ=255)
                         cGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                         >  2:3380263/1‑150 (MQ=255)
                          gATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGt                                                                                                                        >  1:22619/1‑150 (MQ=255)
                           aTCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGtt                                                                                                                       <  2:2663865/150‑1 (MQ=255)
                            tCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTc                                                                                                                      <  1:4548548/150‑1 (MQ=255)
                            tCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTc                                                                                                                      >  2:2667342/1‑150 (MQ=255)
                              aaTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGc                                                                                                                    >  1:509844/1‑150 (MQ=255)
                                tCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTg                                                                                                                  >  2:679977/1‑150 (MQ=255)
                                tCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTg                                                                                                                  <  2:4004878/150‑1 (MQ=255)
                                 cccTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGa                                                                                                                 <  2:1001425/150‑1 (MQ=255)
                                   ctctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAAt                                                                                                               <  2:4155597/150‑1 (MQ=255)
                                    tctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATa                                                                                                              <  1:2295656/150‑1 (MQ=255)
                                      ttCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATaaa                                                                                                            >  2:4453322/1‑150 (MQ=255)
                                        cGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGc                                                                                                          <  2:4091484/150‑1 (MQ=255)
                                                 ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 <  1:1999356/150‑1 (MQ=255)
                                                 ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 >  1:1002378/1‑150 (MQ=255)
                                                 ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 <  1:3333788/150‑1 (MQ=255)
                                                 ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgctt                                                                                                 <  1:2599133/150‑1 (MQ=255)
                                                   aGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgcttgc                                                                                               >  2:725081/1‑150 (MQ=255)
                                                   aGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATATTTGCTGTTCGCTGAATAAAGCGTGttgcttgc                                                                                               >  2:2492325/1‑150 (MQ=255)
                                                    gCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTgcg                                                                                              <  1:1935281/150‑1 (MQ=255)
                                                         acaTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGt                                                                                         >  2:3134192/1‑150 (MQ=255)
                                                            ttAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGc                                                                                      <  1:2876261/150‑1 (MQ=255)
                                                             tAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCt                                                                                     <  2:1202026/150‑1 (MQ=255)
                                                               aGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCttt                                                                                   >  2:1799860/1‑150 (MQ=255)
                                                                  gTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCaa                                                                                <  1:4406918/150‑1 (MQ=255)
                                                                  gTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCaa                                                                                <  2:841795/150‑1 (MQ=255)
                                                                  gTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCaa                                                                                <  2:1302068/150‑1 (MQ=255)
                                                                    aTCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGc                                                                              <  1:1539899/150‑1 (MQ=255)
                                                                        aGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCtt                                                                          <  1:359034/150‑1 (MQ=255)
                                                                        aGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCtt                                                                          <  2:1943439/150‑1 (MQ=255)
                                                                          cGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTcc                                                                        <  1:2984413/150‑1 (MQ=255)
                                                                           gTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCa                                                                       >  2:2177077/1‑150 (MQ=255)
                                                                            ttATCATTGTTGTGCCTGGCGAATCTGTGAACCGACTTCCTGTAGTCTGGCTTCAATGCGAGATAACCCTTAGCCAGGGATAAACAGAAATTTTTGCTGTTCGCTGAATAACGCTTGTTTCTTGCGCCAGTCGCTTTCAATCTCTTCCAg                                                                      <  2:2536225/150‑1 (MQ=255)
                                                                             tAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgt                                                                     <  1:1870742/150‑1 (MQ=255)
                                                                              aGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgtt                                                                    <  2:1507703/150‑1 (MQ=255)
                                                                                cATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgt                                                                  <  1:4480437/150‑1 (MQ=255)
                                                                                 aTTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgtt                                                                 <  1:4079270/150‑1 (MQ=255)
                                                                                 aTTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgtggtt                                                                 >  2:3410418/1‑150 (MQ=255)
                                                                                    gttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttgtt                                                              >  2:4123638/1‑150 (MQ=255)
                                                                                    gttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttgtt                                                              >  2:3558006/1‑150 (MQ=255)
                                                                                      tgtgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGtttt                                                            >  1:2026043/1‑150 (MQ=255)
                                                                                        tgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAt                                                          >  1:2414335/1‑150 (MQ=255)
                                                                                          ccTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAttt                                                        <  1:159139/150‑1 (MQ=255)
                                                                                           cTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTc                                                       <  2:512938/150‑1 (MQ=255)
                                                                                           cTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTc                                                       <  2:3002847/150‑1 (MQ=255)
                                                                                           cTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGAGTGTTGCTTGCTCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTc                                                       >  2:2729094/1‑150 (MQ=255)
                                                                                              gCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCAcc                                                    <  2:1099488/150‑1 (MQ=255)
                                                                                               cGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCa                                                   >  2:3074856/1‑150 (MQ=255)
                                                                                                gAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAg                                                  >  2:4372310/1‑150 (MQ=255)
                                                                                                 aaTCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAgg                                                 >  1:3536015/1‑150 (MQ=255)
                                                                                                    cTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAg                                              >  1:4093004/1‑150 (MQ=255)
                                                                                                    cTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAg                                              <  1:3241461/150‑1 (MQ=255)
                                                                                                    cTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAg                                              <  1:3006560/150‑1 (MQ=255)
                                                                                                     tgtgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGa                                             <  2:2021698/150‑1 (MQ=255)
                                                                                                     tgtgCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGa                                             <  2:1278200/150‑1 (MQ=255)
                                                                                                        gCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCt                                          <  2:1362407/150‑1 (MQ=255)
                                                                                                              aCATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGcc                                    >  2:734093/1‑150 (MQ=255)
                                                                                                                  ccTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATg                                >  1:3476868/1‑150 (MQ=255)
                                                                                                                    tGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGaa                              >  1:2176173/1‑150 (MQ=255)
                                                                                                                     gTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAc                             <  1:2667342/150‑1 (MQ=255)
                                                                                                                      tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                            >  1:1065491/1‑150 (MQ=255)
                                                                                                                      tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                            <  1:2564512/150‑1 (MQ=255)
                                                                                                                      tAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                            <  1:4333972/150‑1 (MQ=255)
                                                                                                                       aGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCg                           >  1:3073226/1‑150 (MQ=255)
                                                                                                                          ctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGa                        >  1:2768315/1‑150 (MQ=255)
                                                                                                                            ggctTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCGGGCAGACCTTCAGCCGATGAACCGGGAAg                      <  2:591172/150‑1 (MQ=255)
                                                                                                                            ggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAg                      <  2:382951/150‑1 (MQ=255)
                                                                                                                             gctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGt                     >  2:4508434/1‑150 (MQ=255)
                                                                                                                                 cAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCa                 <  1:1767303/150‑1 (MQ=255)
                                                                                                                                   aTGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGc               <  1:854770/150‑1 (MQ=255)
                                                                                                                                        aCCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACt          <  1:1664105/150‑1 (MQ=255)
                                                                                                                                         ccTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTg         >  2:1442469/1‑150 (MQ=255)
                                                                                                                                             aCCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCtt     >  2:3243126/1‑150 (MQ=255)
                                                                                                                                             aCCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCtt     <  2:4445452/150‑1 (MQ=255)
                                                                                                                                                aGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTGTa  >  2:4314161/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: