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,0280CA97.6% 774.8 / 8.2 249D476Y (GAC→TAC) ravAhexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base A (133/110);  total (134/115)
Fisher's exact test for biased strand distribution p-value = 9.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01
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 249 total aligned reads displayed.

CAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTG  >  CP009273/3922883‑3923172
                                                                                                                                                 |                                                                                                                                                
cagcagTGGGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGc                                                                                                                                                <  1:807818/148‑1 (MQ=255)
  gcagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCa                                                                                                                                              >  2:2557030/1‑148 (MQ=255)
   cagTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAg                                                                                                                                             >  2:1931897/1‑148 (MQ=255)
     gTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAggg                                                                                                                                           <  2:184846/148‑1 (MQ=255)
     gTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAggg                                                                                                                                           <  1:2840452/148‑1 (MQ=255)
           tgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                     <  2:2526344/148‑1 (MQ=255)
           tgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                     >  2:2490222/1‑148 (MQ=255)
           tgatCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAAc                                                                                                                                     >  2:2763292/1‑148 (MQ=255)
              tcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGa                                                                                                                                  >  1:2298921/1‑148 (MQ=255)
              tcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGa                                                                                                                                  >  1:2297568/1‑148 (MQ=255)
              tcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGa                                                                                                                                  >  2:2278757/1‑148 (MQ=255)
              tcatcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGa                                                                                                                                  >  2:307447/1‑148 (MQ=255)
                atcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACa                                                                                                                                <  2:1974492/148‑1 (MQ=255)
                atcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACa                                                                                                                                <  1:1905460/148‑1 (MQ=255)
                atcTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACa                                                                                                                                >  2:1060101/1‑148 (MQ=255)
                 tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAt                                                                                                                               <  2:1956679/148‑1 (MQ=255)
                 tctcTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAt                                                                                                                               >  2:2993163/1‑148 (MQ=255)
                   tctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAttt                                                                                                                             <  2:1754820/148‑1 (MQ=255)
                   tctTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACAttt                                                                                                                             <  1:3245921/148‑1 (MQ=255)
                      tCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTg                                                                                                                          >  1:2654101/1‑148 (MQ=255)
                      tCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTg                                                                                                                          <  2:2567151/148‑1 (MQ=255)
                         aTCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                       <  2:1197526/148‑1 (MQ=255)
                         aTCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTg                                                                                                                       <  1:783534/148‑1 (MQ=255)
                          tCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGt                                                                                                                      <  1:2052422/148‑1 (MQ=255)
                           cAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGtt                                                                                                                     >  1:1340827/1‑148 (MQ=255)
                            aaTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTc                                                                                                                    >  2:2483993/1‑148 (MQ=255)
                             aTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCg                                                                                                                   <  2:1346518/148‑1 (MQ=255)
                             aTCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCg                                                                                                                   <  1:3457738/148‑1 (MQ=255)
                               cccTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCt                                                                                                                 >  2:144842/1‑148 (MQ=255)
                                 ctctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGa                                                                                                               >  2:3309073/1‑148 (MQ=255)
                                  tctTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATTCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGaa                                                                                                              >  1:2437053/1‑148 (MQ=255)
                                     tCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATaa                                                                                                           >  1:3554177/1‑148 (MQ=255)
                                        cTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGACTAAAGc                                                                                                        >  1:2200008/1‑148 (MQ=255)
                                        cTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGc                                                                                                        >  1:1108502/1‑148 (MQ=255)
                                         tGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCg                                                                                                       <  1:1498921/148‑1 (MQ=255)
                                           aCGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgtg                                                                                                     >  2:2145210/1‑148 (MQ=255)
                                            cGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCgtgt                                                                                                    >  2:803354/1‑148 (MQ=255)
                                             ggCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGtt                                                                                                   >  1:3170570/1‑148 (MQ=255)
                                               ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgc                                                                                                 >  2:3225460/1‑148 (MQ=255)
                                               ccAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGttgc                                                                                                 <  2:3403701/148‑1 (MQ=255)
                                                     tCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTgcgc                                                                                           >  1:1457093/1‑148 (MQ=255)
                                                        caTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAg                                                                                        >  1:407885/1‑148 (MQ=255)
                                                         aTTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGt                                                                                       <  2:2175480/148‑1 (MQ=255)
                                                          ttAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTc                                                                                      >  2:112805/1‑148 (MQ=255)
                                                          ttAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTc                                                                                      >  2:571525/1‑148 (MQ=255)
                                                            aaGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTCACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGc                                                                                    <  2:354293/148‑1 (MQ=255)
                                                             aGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCt                                                                                   >  2:1966166/1‑148 (MQ=255)
                                                              gCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCtt                                                                                  <  2:594341/148‑1 (MQ=255)
                                                               cGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCttt                                                                                 <  1:1750126/148‑1 (MQ=255)
                                                               cGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCttt                                                                                 <  2:2951185/148‑1 (MQ=255)
                                                                gTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTc                                                                                <  1:578958/148‑1 (MQ=255)
                                                                  aTCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCaa                                                                              >  1:2569317/1‑148 (MQ=255)
                                                                  aTCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCaa                                                                              >  1:3201683/1‑148 (MQ=255)
                                                                   tCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAg                                                                             <  1:518496/148‑1 (MQ=255)
                                                                    ccAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCTACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGc                                                                            <  2:3517231/148‑1 (MQ=255)
                                                                          ttAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTcc                                                                      <  2:1515679/148‑1 (MQ=255)
                                                                              cATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgtt                                                                  <  2:3323107/148‑1 (MQ=255)
                                                                                ttgttgTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgt                                                                >  1:2667528/1‑148 (MQ=255)
                                                                                ttgttgTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgt                                                                >  1:1322841/1‑148 (MQ=255)
                                                                                ttgttgTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgt                                                                >  2:151005/1‑148 (MQ=255)
                                                                                ttgttgTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgt                                                                >  1:3128287/1‑148 (MQ=255)
                                                                                ttgttgTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgt                                                                >  2:2320082/1‑148 (MQ=255)
                                                                                   ttgtGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAANCAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttgt                                                             <  1:1805535/148‑1 (MQ=255)
                                                                                    tgtgCCTGGCGAATCTGTGCACCTACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAgttgttgtt                                                            <  2:1293546/148‑1 (MQ=255)
                                                                                      tgCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGtttt                                                          >  2:626134/1‑148 (MQ=255)
                                                                                       gCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTa                                                         <  1:840894/148‑1 (MQ=255)
                                                                                        ccTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAt                                                        <  1:1166924/148‑1 (MQ=255)
                                                                                         cTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAtt                                                       <  1:2908925/148‑1 (MQ=255)
                                                                                          tGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTAttt                                                      <  2:1698388/148‑1 (MQ=255)
                                                                                           ggCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTc                                                     <  1:2161887/148‑1 (MQ=255)
                                                                                            gCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCa                                                    >  1:592259/1‑148 (MQ=255)
                                                                                            gCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCa                                                    >  1:2380861/1‑148 (MQ=255)
                                                                                             cGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCAc                                                   <  1:3334871/148‑1 (MQ=255)
                                                                                               aaTCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCa                                                 >  2:1031361/1‑148 (MQ=255)
                                                                                                aTCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAg                                                >  2:2680153/1‑148 (MQ=255)
                                                                                                 tCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAgg                                               <  1:1512674/148‑1 (MQ=255)
                                                                                                         ccGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCtt                                       <  2:2190397/148‑1 (MQ=255)
                                                                                                         ccGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCtt                                       <  2:394507/148‑1 (MQ=255)
                                                                                                            aCATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAg                                    >  1:641206/1‑148 (MQ=255)
                                                                                                              aTCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGcc                                  >  1:398585/1‑148 (MQ=255)
                                                                                                                 cTGTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGAt                               <  2:234804/148‑1 (MQ=255)
                                                                                                                   gTAGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGa                             >  2:2434444/1‑148 (MQ=255)
                                                                                                                     aGGCTGGCTTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAc                           <  1:2815616/148‑1 (MQ=255)
                                                                                                                      ggctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAAcc                          <  2:93402/148‑1 (MQ=255)
                                                                                                                       gctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCCGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCg                         >  1:2360335/1‑148 (MQ=255)
                                                                                                                       gctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCg                         >  1:364963/1‑148 (MQ=255)
                                                                                                                       gctggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCg                         <  1:103564/148‑1 (MQ=255)
                                                                                                                         tggctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCggg                       <  1:2246009/148‑1 (MQ=255)
                                                                                                                           gctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGaa                     <  1:3309073/148‑1 (MQ=255)
                                                                                                                            ctTCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAg                    <  1:1588781/148‑1 (MQ=255)
                                                                                                                              tCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTg                  >  2:1606579/1‑148 (MQ=255)
                                                                                                                              tCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTg                  <  1:3391253/148‑1 (MQ=255)
                                                                                                                              tCAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTg                  >  1:1084394/1‑148 (MQ=255)
                                                                                                                               cAATGCGACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGc                 <  1:3260888/148‑1 (MQ=255)
                                                                                                                                     gACCTAACCCGTAGCCAGGGATAAACAGACATTTTTGCTGGTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATGTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGt           <  2:2202137/148‑1 (MQ=255)
                                                                                                                                     gACCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGt           <  1:2731355/148‑1 (MQ=255)
                                                                                                                                      aCCTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTa          <  1:1931897/148‑1 (MQ=255)
                                                                                                                                        cTAACCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCATCAGCCGATGAACCGGGAAGTGCCCGCGTACt        >  1:1260905/1‑148 (MQ=255)
                                                                                                                                           aCCAGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGcc     <  1:3006055/148‑1 (MQ=255)
                                                                                                                                              aGTAGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTg  >  2:1706053/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAGCAGTGCGATGATCATCTCTTCGATCAATCCCTCTTCGCTGACGGCCAGCATCACATTAAGCGTATCCAGCGTTAGCATTGTTGTGCCTGGCGAATCTGTGCACCGACATCCTGTAGGCTGGCTTCAATGCGACCTAACCAGTCGCCAGGGATAAACAGACATTTTTGCTGTTCGCTGAATAAAGCGTGTTGCTTGCGCCAGTCGCTTTCAAGCTCTTCCAGTTGTTGTTTTATTTCACCAGGCAGACCTTCAGCCGATGAACCGGGAAGTGCCAGCGTACTGCCTTG  >  CP009273/3922883‑3923172

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: