Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 776,394 A→T 8.2% intergenic (+20/‑130) lysW → / → valZ tRNA‑Lys/tRNA‑Val

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273776,3940AT 8.2% 560.7 / 8.2 223intergenic (+20/‑130)lysW/valZtRNA‑Lys/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base A (94/97);  new base T (2/21);  total (99/124)
Fisher's exact test for biased strand distribution p-value = 2.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

CACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTG  >  CP009273/776245‑776541
                                                                                                                                                     |                                                                                                                                                   
  ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTTGCTCAGTTTGTAGAGCAGTTTACTTTTAATATATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATg                                                                                                                                                   <  1:2449173/150‑1 (MQ=14)
  ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTTGCTCAGGTGGTTGAGCAGTTGTCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATg                                                                                                                                                   <  1:3896377/150‑1 (MQ=255)
    tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGtt                                                                                                                                                 <  1:1658908/150‑1 (MQ=255)
    tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGtt                                                                                                                                                 <  1:2595689/150‑1 (MQ=40)
     cccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTa                                                                                                                                                <  2:2361755/150‑1 (MQ=40)
         tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATAAATTTGTCGCCGGGTCGAATACAGCACCACCCACCACTTTTTACATCGAAGCCAGATGTTTAgag                                                                                                                                            >  1:4126250/1‑148 (MQ=11)
            aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTGAATCAAGTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACCGATGTTAAGCgagt                                                                                                                                         >  2:4200001/1‑150 (MQ=14)
             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                                                                                                                        <  2:1826441/150‑1 (MQ=255)
             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                                                                                                                        <  1:4577095/150‑1 (MQ=255)
              ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                                                                                                                       <  1:524479/150‑1 (MQ=255)
                 gggggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGAt                                                                                                                                    <  2:3790690/150‑1 (MQ=255)
                   gggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATaa                                                                                                                                  <  2:1498059/150‑1 (MQ=255)
                       cggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGtt                                                                                                                              <  1:1552831/150‑1 (MQ=255)
                       cggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGtt                                                                                                                              <  1:1553680/150‑1 (MQ=255)
                         gcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCCGTTGGTAGAGCTGTTGTCTTTTCATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGc                                                                                                                            <  1:2059514/150‑1 (MQ=255)
                           ggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTGAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTTTTAAAATAGAAGACAGATGTTAATCGTGTAGGATAAAGTTTCGt                                                                                                                          >  1:1779989/1‑150 (MQ=14)
                                gATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCa                                                                                                                     >  2:1727476/1‑150 (MQ=255)
                                gATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACCGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCa                                                                                                                     >  2:2630232/1‑150 (MQ=255)
                                gATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGGAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGTCCCACCAGTTTTAAAATCGAAGACAGCTGGTAAGCGTGTAGGATAACGTTGCGgcagca                                                                                                                     >  2:4024218/1‑150 (MQ=14)
                                        catcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCCGTTTTAAACTCGAAGACAGATGTTAAGCGTGTAGGATAACCTTGCGTCAGCAACGGCCCg                                                                                                             >  2:4483210/1‑150 (MQ=255)
                                         atcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGTGCTGTTTACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACTGATTTTAAGCGTGTAGGATAACGTTGCTTCAGCAACGGCCCGt                                                                                                            <  2:1111321/150‑1 (MQ=14)
                                              ccaccaCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgc                                                                                                       >  2:2823184/1‑150 (MQ=255)
                                                        gccGTTCGCTCCTTTGTTAGACCTTTTGACTTTTAATCAATTGGTCGCTGTTTCTAATCCTGCACCCCCCACCCTTTTTAACACCCAACACCGATTTTAAGCGTGTTGTATAAGTTGGCGTCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                             <  2:3247253/148‑1 (MQ=2)
                                                            ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                         <  2:2245110/150‑1 (MQ=255)
                                                             tcgCTCAGTTGGTAGAGCAGTTGACTTTTAACCAATTGGTCGCAGGTTCTAATCCGGCACGACCCACCAGTTTTAACATCGAACACAGATGTTAAGCGTGTAGTATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                        <  1:3647820/148‑1 (MQ=11)
                                                                cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATcc                                                                                     >  1:6477/1‑150 (MQ=255)
                                                                   tgtgggTAGAGCTGTTGCCTTTTAATCATTGGGTCGCTGGTTCGAATCCTGCCCGCCCCCCCCGTTTTACCATCGAAGACCGTTGTTAAGCGGGTCGGATAACGTTGCGTCCGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGc                                                                                  <  2:3665992/146‑1 (MQ=14)
                                                                     ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                                >  1:2202632/1‑150 (MQ=255)
                                                                               aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAGGTTATGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCcaccac                                                                      >  2:737059/1‑150 (MQ=255)
                                                                                gttcccTTTTCATCAATTGGTCGCAGTAGCTACTCCTGCACCACCCCCCCGTTTTAACATCTAACACAGTTTTTAAGCGTTTAGTATAACGTTGCGTCAGAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACt                                                                     <  1:624511/145‑1 (MQ=11)
                                                                                 ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTa                                                                                                                                                <  1:3837136/74‑1 (MQ=255)
                                                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAAt                                                                  >  1:279422/1‑150 (MQ=255)
                                                                                     cTTTTAATTAATTGGTCGCTGGTTCTAATCCTTCACGACCAAACTTTTTTAACAAGTACCACAGATGTTAATCGTGTAGTATATCTTTGTGTCAGAATCGGCCATTTGGGCGAGCGATTCGAGTCATCCTGCACGACCCACCACTAATGa                                                                <  2:67860/150‑1 (MQ=9)
                                                                                       tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACg                                                              <  2:3606813/150‑1 (MQ=255)
                                                                                       tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACg                                                              >  2:4772939/1‑150 (MQ=255)
                                                                                         tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGt                                                            >  1:1578437/1‑150 (MQ=255)
                                                                                         tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGt                                                            >  1:1401711/1‑150 (MQ=255)
                                                                                          aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTg                                                           >  2:2762493/1‑150 (MQ=255)
                                                                                          aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTg                                                           >  2:3193684/1‑150 (MQ=255)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGtt                                                       >  1:2533535/1‑150 (MQ=255)
                                                                                                 tGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCg                                                                                                                           >  1:2837272/1‑79 (MQ=255)
                                                                                                 tGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGAAGGGCTAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCgg                                                    >  1:404902/1‑150 (MQ=255)
                                                                                                 tGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTTCACCTCGAAGACCGAGGATGATCGTGGAGGATACCCTTGCGTCAGCGCCGGAGCCGAGGCCGAGCGAGGCGACAACTCCTGCACGCCCCACCACTAAAGACCGGGGGTTCgg                                                    >  2:2076342/1‑150 (MQ=9)
                                                                                                     cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGa                                                >  2:1850158/1‑150 (MQ=255)
                                                                                                      gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGaa                                               <  2:4039807/150‑1 (MQ=255)
                                                                                                        aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagt                                             >  2:1798471/1‑150 (MQ=255)
                                                                                                        aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagt                                             >  2:2907491/1‑150 (MQ=255)
                                                                                                         ggTTCGAATCCTGCACGACCCACCATTTTTAACATCGAAGACAGTTTTTATGGGTGTTGGATAACGTTGCTTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagta                                            <  1:1012492/150‑1 (MQ=18)
                                                                                                          gTTCGAATCCTGCACGACCCACCCGTTTTAACAGCGAAGCCAGAGGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagtag                                           <  2:4543744/150‑1 (MQ=255)
                                                                                                            tCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACCGATGTTAAGGGTGTAGTATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGtt                                         <  1:3512570/150‑1 (MQ=255)
                                                                                                             cGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGACGTAGttt                                        >  2:1178979/1‑150 (MQ=255)
                                                                                                              gAAGCCTGCACGCCCCACCCGTTTTAACATCGAAGACAGCTGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTACTGACGGTGGGTTCGGTGGAAGTAGTTTg                                       <  2:3891064/150‑1 (MQ=18)
                                                                                                               aaTCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgt                                      >  1:4379740/1‑150 (MQ=255)
                                                                                                               aaTCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGGCAGCAACGGCCCGTAGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCACTCTTCACGCTGGGTTCGGTGGACGTAGGTTgg                                      >  2:2440201/1‑149 (MQ=17)
                                                                                                                 tCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtag                                    >  2:2351416/1‑150 (MQ=255)
                                                                                                                 tCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtag                                    >  2:2283059/1‑150 (MQ=255)
                                                                                                                   cTGCACGACCCACCATTTTTACCTTCTCAGACTGTTGTTACGGGTGTAGTATACCGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtagta                                  <  1:1743676/150‑1 (MQ=21)
                                                                                                                      cACGACCCACCAGTTTTACCATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                               <  2:1433510/150‑1 (MQ=255)
                                                                                                                      cACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                               >  2:2133594/1‑150 (MQ=255)
                                                                                                                      cACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                               >  1:1183165/1‑150 (MQ=255)
                                                                                                                       aCGACCCACCTTTTTTAACACCGAACACAGATGTTAAGCGTGTAGTATACCGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCa                              <  2:1624536/150‑1 (MQ=21)
                                                                                                                        cGACCCACCTGTTTAAACATCGAAGAACGTTTTTATGGGTGTTGGATAACGTGGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAg                             <  1:4734991/150‑1 (MQ=25)
                                                                                                                        cGACCCACCCGTTGTACCATCGACGACCGCTGTTACGCGTGTAGGATACCGGGGCGGCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCGGCACGACCCACCCCGAATGACGGTGGGTTCGGTGGAAGGAGTTTGTAGTATCCAg                             <  2:1893085/150‑1 (MQ=21)
                                                                                                                        cGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGGGGGTTCGGTGGAAGTAGTTGGTAGTATcccg                             >  1:3347955/1‑148 (MQ=255)
                                                                                                                         gACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgc                            >  1:1621145/1‑150 (MQ=255)
                                                                                                                          aCCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcg                           <  2:2533535/150‑1 (MQ=255)
                                                                                                                          aCCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcg                           <  2:4077310/150‑1 (MQ=255)
                                                                                                                           cccACCAGTTTTAACATCGAAGACAGCTGTTAAGCGTGTAGGATACCGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcgc                          <  1:3576945/150‑1 (MQ=255)
                                                                                                                              acttttttaAAAATCTAAGACATTTTTTATGGGTGTTGTTTAACGTTGCGACTGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGt                       <  1:3704597/145‑1 (MQ=21)
                                                                                                                              accaGTTTTAACATCGAAGACAGATGTTAAGCGTGTTGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGt                       <  1:4330765/150‑1 (MQ=255)
                                                                                                                               cctgTTTTAACATCGAAGACCGTGGTTAAGGGTGTAGGATACGGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa                      <  1:784573/147‑1 (MQ=25)
                                                                                                                               ccaGTTTTAACTTGGATGACAGATGTTATGCGTGTTGGATATGGTTGCGTCAGAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa                      <  1:2416390/150‑1 (MQ=25)
                                                                                                                               ccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGTATAACTTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa                      <  1:375452/150‑1 (MQ=255)
                                                                                                                               ccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGTATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa                      <  1:1522962/150‑1 (MQ=255)
                                                                                                                                ctgTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTAt                     <  2:3860099/148‑1 (MQ=255)
                                                                                                                                 aGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATc                    >  1:2134981/1‑150 (MQ=255)
                                                                                                                                 aGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATc                    <  2:2511420/150‑1 (MQ=255)
                                                                                                                                  gTTTTAACATCGACGCCCGATGTTATGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGCCGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCg                   <  2:4509932/150‑1 (MQ=25)
                                                                                                                                    tttAACATCGAAGACCGATGTTATGCGTGTTGTATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCggg                 <  1:2630232/150‑1 (MQ=255)
                                                                                                                                    tttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCggg                 >  2:1320027/1‑150 (MQ=255)
                                                                                                                                      tAACATCGACGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTg               <  2:4375212/150‑1 (MQ=255)
                                                                                                                                      tAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTg               >  2:2622849/1‑150 (MQ=255)
                                                                                                                                      tAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTg               >  1:1763390/1‑150 (MQ=255)
                                                                                                                                       aaCATCGAAGACACATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa              >  2:159047/1‑150 (MQ=255)
                                                                                                                                       aaCAGCGAAGACAGATGTTATGGGTGTCGTATAACTTTGCGTCAGCCACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa              <  1:599789/150‑1 (MQ=37)
                                                                                                                                         cATCGAAGACAGATTTTAAGCGTGTAGGATACCGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGAtt            <  1:3071656/150‑1 (MQ=255)
                                                                                                                                         aaTCGAAGACAGTTGTTAAGCGTGTAGTATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGAtt            <  1:2164120/149‑1 (MQ=255)
                                                                                                                                          accgacgaCAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTa           <  2:4551065/148‑1 (MQ=255)
                                                                                                                                           tCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAg          >  2:518178/1‑150 (MQ=255)
                                                                                                                                             gAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGGGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCt        >  1:3715987/1‑150 (MQ=255)
                                                                                                                                              acgACAGCTGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTc       <  2:698231/148‑1 (MQ=255)
                                                                                                                                              aaGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTc       >  1:859737/1‑150 (MQ=255)
                                                                                                                                                gACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTGGTAGTATCCAGCGCAGTATCGGGTGATTATCTCAg     >  2:4292663/1‑150 (MQ=255)
                                                                                                                                                  cAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCt   <  2:2202699/150‑1 (MQ=255)
                                                                                                                                                  cAGATGTAAAGCGTGTAGGATACGGTTGCGTCACCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCt   <  1:1344423/150‑1 (MQ=255)
                                                                                                                                                   ataTGTTATGCGTGTTGGATAAGGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTg  <  1:174346/148‑1 (MQ=255)
                                                                                                                                                   aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACTGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTg  >  1:3242668/1‑150 (MQ=255)
                                                                                                                                                   aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTg  >  2:3317/1‑150 (MQ=255)
                                                                                                                                                   aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTg  >  2:891643/1‑150 (MQ=255)
                                                                                                                                                   aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGAGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTg  >  2:496437/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTG  >  CP009273/776245‑776541

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