Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 776,389 G→C 6.4% intergenic (+15/‑135) lysW → / → valZ tRNA‑Lys/tRNA‑Val

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273776,3890GC 6.4% 646.8 / 4.7 219intergenic (+15/‑135)lysW/valZtRNA‑Lys/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base G (94/105);  new base C (2/13);  total (97/122)
Fisher's exact test for biased strand distribution p-value = 1.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

AGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAG  >  CP009273/776243‑776538
                                                                                                                                                  |                                                                                                                                                     
agCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCCGTTTTAACATCGAAGACa                                                                                                                                                    >  2:1522962/1‑150 (MQ=255)
  cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGa                                                                                                                                                  <  2:2237200/150‑1 (MQ=40)
   aCCTCCCTTACACGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTTGCTGTGTTGTTAGAGGTGTTGACTTTTAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                                                                                                                 <  1:2161553/150‑1 (MQ=14)
    ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCCGGTGGTTGAGCAGGTGACTTTTAATAAATTGGTCGCAGGGTCGAAACCTGCACGCCCCCCCAGTTTTAACAAAGAATAAAGATg                                                                                                                                                >  1:2511420/1‑150 (MQ=11)
      tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGtt                                                                                                                                              <  1:1658908/150‑1 (MQ=255)
       cccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTa                                                                                                                                             <  2:2361918/150‑1 (MQ=40)
         cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAg                                                                                                                                           <  1:3337264/150‑1 (MQ=40)
           tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAATCg                                                                                                                                         >  2:32905/1‑150 (MQ=255)
              aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTATCGCCGGTGGGAGAGCCGTTTACTTGGAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAGTTTTAACCTCGGAGACCGATGTTAAGAgtgt                                                                                                                                      >  2:4401230/1‑150 (MQ=14)
              aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTGAATCAAGTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACCGATGTTAAGCgagt                                                                                                                                      >  2:4200001/1‑150 (MQ=14)
               aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                                                                                                                     <  1:4577095/150‑1 (MQ=255)
               aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                                                                                                                     <  2:1826441/150‑1 (MQ=255)
                ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                                                                                                                    <  1:524479/150‑1 (MQ=255)
                 gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAgg                                                                                                                                   <  2:1239185/150‑1 (MQ=255)
                    ggggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATa                                                                                                                                <  2:2600350/150‑1 (MQ=255)
                    ggggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATa                                                                                                                                <  2:170296/150‑1 (MQ=255)
                     gggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATaa                                                                                                                               <  2:1498059/150‑1 (MQ=255)
                         cggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGtt                                                                                                                           <  1:1552928/150‑1 (MQ=255)
                             ggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCCATTGGTCGCAGGTTAGAATACTGCACGACCCACCAGTTTTAACATAGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCgg                                                                                                                       >  1:1779984/1‑149 (MQ=255)
                               ttCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                                                                                                                     <  1:2106047/150‑1 (MQ=255)
                                  gATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACCGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCa                                                                                                                  >  2:2630232/1‑150 (MQ=255)
                                  gATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGGAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGTCCCACCAGTTTTAAAATCGAAGACAGCTGGTAAGCGTGTAGGATAACGTTGCGgcagca                                                                                                                  >  2:4024218/1‑150 (MQ=14)
                                   aTCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa                                                                                                                 <  2:3870647/150‑1 (MQ=255)
                                       cGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGc                                                                                                             <  1:634553/150‑1 (MQ=255)
                                           atcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGTGCTGTTTACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACTGATTTTAAGCGTGTAGGATAACGTTGCTTCAGCAACGGCCCGt                                                                                                         <  2:1111321/150‑1 (MQ=14)
                                                ccaccaCCGGGTCGTTAGCTCAGTTGGGAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAGTTTTAACATCGCAGACAGTTGTTAAGCGTGTAGGATAACGTGGCGTCAGCAACGGCCCGTAGGgc                                                                                                    >  2:3008518/1‑150 (MQ=14)
                                                  accaccGGGTCGTTCGCTCCGTTGGTAGAGCATTTGACTTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCTGTTTTAACATCGAAGACAGTTTTTATGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcga                                                                                                  <  1:2771906/150‑1 (MQ=14)
                                                    caccGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAGTGGTCGCACGATCGAATCCTGCACGACCCACCAGTTTTAACATCGCACAAACAAGGTAAGCGTGGAGGATAACGAAGCGTCAGCAACGGCCCGGCGGgcgagc                                                                                                >  2:3760338/1‑150 (MQ=11)
                                                           tctttAGCACAGTTGTTACAGATTTTTAATTTTAAAAAATTGGTCGCAGTTTCGAATCCTGCACTACCCACAATTTTTAACATCTAAGACAGATTTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                         <  1:2822705/147‑1 (MQ=14)
                                                              ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                      <  2:2245328/150‑1 (MQ=255)
                                                               tAGCTCATTTGGTAGAGCAGTTTACTTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCATTTTTAACATCGAACACAGATTTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                     <  1:3646619/150‑1 (MQ=11)
                                                                     tgtgggTAGAGCTGTTGCCTTTTAATCATTGGGTCGCTGGTTCGAATCCTGCCCGCCCCCCCCGTTTTACCATCGAAGACCGTTGTTAAGCGGGTCGGATAACGTTGCGTCCGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGc                                                                               <  2:3665992/146‑1 (MQ=14)
                                                                               gctgTTGACTTTTAATATTTGGGCCGCTGGTTCCAAACCTGCACGACCCACAAGTTTAAAAATCGAACACTTATTTTAAGGGTTTTGTATAAGGTTGCGTCACCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                                     <  1:3856279/147‑1 (MQ=2)
                                                                                  gttcccTTTTCATCAATTGGTCGCAGTAGCTACTCCTGCACCACCCCCCCGTTTTAACATCTAACACAGTTTTTAAGCGTTTAGTATAACGTTGCGTCAGAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACt                                                                  <  1:624511/145‑1 (MQ=11)
                                                                                   ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTa                                                                                                                                             >  2:3837136/1‑74 (MQ=255)
                                                                                      aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACACACCAGTTTTAACATCGAAGACATATTTTATGCTTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATg                                                              <  2:2520227/150‑1 (MQ=14)
                                                                                       cTTTTAATTAATTGGTCGCTGGTTCTAATCCTTCACGACCAAACTTTTTTAACAAGTACCACAGATGTTAATCGTGTAGTATATCTTTGTGTCAGAATCGGCCATTTGGGCGAGCGATTCGAGTCATCCTGCACGACCCACCACTAATGa                                                             <  2:67860/150‑1 (MQ=9)
                                                                                        ttttAAACTTTTGGTCGCTTGTTCTAATCCTGCACGACACACCATTTTTAACTTCGACGAGTTTTTTTATGGGTTTTTGATATGGTTGGGTCAGAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGAc                                                            <  1:3126343/150‑1 (MQ=12)
                                                                                         tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACg                                                           <  2:3606813/150‑1 (MQ=255)
                                                                                         tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACg                                                           >  2:2763646/1‑150 (MQ=255)
                                                                                         tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACg                                                           >  2:4773781/1‑150 (MQ=255)
                                                                                           tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGt                                                         >  1:1578381/1‑150 (MQ=255)
                                                                                            aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTg                                                        >  2:3193684/1‑150 (MQ=255)
                                                                                             aacacttGGTCGCATGTTCGAATCCTGCCCGACCCACCAGTTTTCACATCGAATACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTgg                                                       <  2:1416968/145‑1 (MQ=17)
                                                                                                 tttGGTCTCTGTTTCTAATCGTGCACGACCAACATTTTTTAAAATCTAATAAAGTTGTTATGGGTTTTGTATAAGTTTGCGTCCCCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTc                                                   <  1:4120184/149‑1 (MQ=12)
                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                         >  1:1828975/1‑128 (MQ=255)
                                                                                                   tGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCg                                                                                                                        <  2:2837272/79‑1 (MQ=255)
                                                                                                   tGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCgg                                                 >  2:1474489/1‑150 (MQ=255)
                                                                                                   tGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTTCACCTCGAAGACCGAGGATGATCGTGGAGGATACCCTTGCGTCAGCGCCGGAGCCGAGGCCGAGCGAGGCGACAACTCCTGCACGCCCCACCACTAAAGACCGGGGGTTCgg                                                 >  2:2076342/1‑150 (MQ=9)
                                                                                                    ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGAAAACGTTGCGACATCAACTCACAAAACCCCTAACCAAGCCAGTAATCCTGCAAGACCCACCACTAATGACGGTGGGTGCggg                                                >  1:3674139/1‑149 (MQ=12)
                                                                                                       cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGa                                             >  2:1850158/1‑150 (MQ=255)
                                                                                                        gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGaa                                            <  2:4039807/150‑1 (MQ=255)
                                                                                                        gCAGGTTCGAATCCTGCACGACCCACCAGTTACAACATCGAAGACAGATGTTATGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGaa                                            <  2:3923476/150‑1 (MQ=255)
                                                                                                        gCAGGTCCTAAGCCGGCACGCCCCACCCGTTTTAACATCGAAGACTGCTGTTAAGCGTGTTGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGaa                                            <  2:842162/150‑1 (MQ=18)
                                                                                                          aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagt                                          >  2:2343522/1‑150 (MQ=255)
                                                                                                          aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagt                                          >  2:1798471/1‑150 (MQ=255)
                                                                                                           gtttCGAAACCGGCACGCCCCACCATTTTTAACATCGAAGACAGATGTTACGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagta                                         <  1:3463046/148‑1 (MQ=18)
                                                                                                           ggTTCGAATCCTGCACGACCCACCATTTTTAACATCGAAGACAGTTTTTATGGGTGTTGGATAACGTTGCTTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagta                                         <  1:1012492/150‑1 (MQ=18)
                                                                                                            gTTCGAATCCTGCACGACCCACCCGTTTTAACAGCGAAGCCAGAGGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagtag                                        <  2:4543744/150‑1 (MQ=255)
                                                                                                             ttCGAATCCTGCACGCCCCACCAGTTTAAACTTCGAAGACAGATTTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagtagt                                       <  1:2198550/150‑1 (MQ=255)
                                                                                                               ctattCCGGCACGACCCCCCAGTTTAAACATCGAAGACAGATGTTTAGCTTGTAGGATAACGTTGCGACCGCTACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGttt                                     <  1:394784/146‑1 (MQ=18)
                                                                                                               cGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGttt                                     >  1:70139/1‑150 (MQ=255)
                                                                                                                 aaTCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGGCAGCAACGGCCCGTAGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCACTCTTCACGCTGGGTTCGGTGGACGTAGGTTgg                                   >  2:2440201/1‑149 (MQ=17)
                                                                                                                   tCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtag                                 >  2:2351659/1‑150 (MQ=255)
                                                                                                                   tCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCACGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGTCGTAGTTTgttg                                 >  1:1485116/1‑148 (MQ=255)
                                                                                                                   tCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAAAGTTCCGTCGTCAACGCCACTTAGCGCGAGACAACCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTGGTGGgtag                                 >  1:4322726/1‑150 (MQ=18)
                                                                                                                    ccTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCCCTAATGACGGGGGGTTCGGTGGAAGTAGTTTgtagt                                >  1:2352459/1‑150 (MQ=255)
                                                                                                                     cTGCACGACCCACCATTTTTACCTTCTCAGACTGTTGTTACGGGTGTAGTATACCGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtagta                               <  1:1743676/150‑1 (MQ=21)
                                                                                                                       gCACGCCCCACCAGTTTTAACATCGACGACAGAGGTTAAGGTTGTAGGATAAGGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATc                             <  2:1741498/150‑1 (MQ=21)
                                                                                                                        cACGACCCACCATTTTTAACTGTGAACCCGGTTTTTACGGGTGTTGGATAAGGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                            <  1:1763131/150‑1 (MQ=21)
                                                                                                                        cACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                            >  1:1183165/1‑150 (MQ=255)
                                                                                                                        cACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                            >  2:4583769/1‑150 (MQ=255)
                                                                                                                        cACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                            >  2:1152180/1‑150 (MQ=255)
                                                                                                                        cACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                            >  2:2133594/1‑150 (MQ=255)
                                                                                                                        cACGACCAACCATTTTTAAAATCTAAGACAGTTTTTATGGGTTTTGTATATGGTTGCTTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc                            <  1:3787269/150‑1 (MQ=21)
                                                                                                                           gACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgc                         >  1:1621145/1‑150 (MQ=255)
                                                                                                                           gACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgc                         >  1:4578516/1‑150 (MQ=255)
                                                                                                                            aCCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcg                        <  2:4077310/150‑1 (MQ=255)
                                                                                                                             cccACCAGTTTTAACATCGAAGACAGCTGTTAAGCGTGTAGGATACCGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcgc                       <  1:3576945/150‑1 (MQ=255)
                                                                                                                             cccACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcgc                       >  1:2428892/1‑150 (MQ=255)
                                                                                                                                accaGTTTTACAATCGAAGACAGATGTTATGCGTGTTGGATTACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGt                    <  1:4331547/150‑1 (MQ=25)
                                                                                                                                acatttttaAAAATCTAAGACTTTTGTTATGGGTTTTGTATAAGGTTGCGTCTGCACCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGt                    <  1:3715443/145‑1 (MQ=21)
                                                                                                                                 ccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGTATAACTTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa                   <  1:375452/150‑1 (MQ=255)
                                                                                                                                 ccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGTATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa                   <  1:1522962/150‑1 (MQ=255)
                                                                                                                                    gTTTTAACATCGACGCCCGATGTTATGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGCCGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCg                <  2:4509932/150‑1 (MQ=25)
                                                                                                                                     ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGAGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg               >  1:1347130/1‑150 (MQ=255)
                                                                                                                                      tttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCggg              <  1:2630192/150‑1 (MQ=255)
                                                                                                                                       ttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGt             >  1:1595366/1‑150 (MQ=255)
                                                                                                                                        tAACATCGACGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTg            <  2:4375212/150‑1 (MQ=255)
                                                                                                                                        tAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTg            >  2:2622879/1‑150 (MQ=255)
                                                                                                                                         aaCATCGAAGACAGATGTTATGCGTGTTGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa           <  1:3032915/150‑1 (MQ=255)
                                                                                                                                         aaCATCGAAGACACATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa           >  2:159047/1‑150 (MQ=255)
                                                                                                                                         aaCAGCGAAGACAGATGTTATGGGTGTCGTATAACTTTGCGTCAGCCACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa           <  1:599789/150‑1 (MQ=37)
                                                                                                                                           cATCGAAGACAGATGTTAAGCGTTTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTACTGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGAtt         <  2:1925418/150‑1 (MQ=255)
                                                                                                                                            accgacgaCAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTa        <  2:4551065/148‑1 (MQ=255)
                                                                                                                                                aaGACAGATGTTAAGCGTGTAGGATACCGGTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTc    <  2:697790/150‑1 (MQ=255)
                                                                                                                                                aaGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTc    >  1:2588185/1‑150 (MQ=255)
                                                                                                                                                aaGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTc    >  1:859737/1‑150 (MQ=255)
                                                                                                                                                aaGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTc    <  2:698905/150‑1 (MQ=255)
                                                                                                                                                  gACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTGGTAGTATCCAGCGCAGTATCGGGTGATTATCTCAg  >  2:4292663/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAG  >  CP009273/776243‑776538

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: