Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,183,817 C→A 6.1% T193K (ACG→AAG)  yqiH → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,183,8170CA 6.1% 760.0 / 19.6 201T193K (ACG→AAG) yqiHputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (86/102);  new base A (6/7);  total (92/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

CATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTAATATTT  >  CP009273/3183668‑3183966
                                                                                                                                                     |                                                                                                                                                     
cATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAc                                                                                                                                                       >  1:4266393/1‑150 (MQ=255)
cATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAc                                                                                                                                                       <  2:4017682/150‑1 (MQ=255)
 aTTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACg                                                                                                                                                      >  1:4051078/1‑150 (MQ=255)
 aTTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCATAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACg                                                                                                                                                      <  1:1266757/150‑1 (MQ=255)
   tGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAAGCa                                                                                                                                                    <  2:3342593/150‑1 (MQ=255)
    gAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAAGCaa                                                                                                                                                   <  1:4394629/150‑1 (MQ=255)
      aTTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAg                                                                                                                                                 >  2:895353/1‑150 (MQ=255)
      aTTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAg                                                                                                                                                 <  1:4458211/150‑1 (MQ=255)
      aTTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAg                                                                                                                                                 <  1:3172287/150‑1 (MQ=255)
         agaagaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCtt                                                                                                                                              <  2:3756488/150‑1 (MQ=255)
          gaagaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTa                                                                                                                                             <  1:3492828/150‑1 (MQ=255)
             gaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATg                                                                                                                                          <  2:4510709/150‑1 (MQ=255)
             gaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAAGCAAAGCTTAATg                                                                                                                                          <  2:770563/150‑1 (MQ=255)
              aCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATgg                                                                                                                                         <  1:2243799/150‑1 (MQ=255)
               cGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGt                                                                                                                                        <  1:517160/150‑1 (MQ=255)
               cGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGt                                                                                                                                        <  1:2673151/150‑1 (MQ=255)
                 aTGGGGTTTTGACTTACGAAATAGTTGAAATCATTCGTAGTAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGTGTTATATCGGAACTAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTaa                                                                                                                                      <  2:1480073/150‑1 (MQ=255)
                 aTGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAAGCAAAGCTTAATGGTaa                                                                                                                                      >  1:615114/1‑150 (MQ=255)
                  tGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAc                                                                                                                                     >  1:48172/1‑150 (MQ=255)
                   ggggTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAcc                                                                                                                                    <  2:2524688/150‑1 (MQ=255)
                   ggggTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAcc                                                                                                                                    <  2:3568514/150‑1 (MQ=255)
                      gTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCt                                                                                                                                 <  2:3251466/150‑1 (MQ=255)
                       tAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCtt                                                                                                                                <  2:1282376/150‑1 (MQ=255)
                        aaTGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCttt                                                                                                                               >  2:2048562/1‑150 (MQ=255)
                        aaTGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCttt                                                                                                                               <  1:1571618/150‑1 (MQ=255)
                         aTGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCtttt                                                                                                                              <  1:4486589/150‑1 (MQ=255)
                         aTGACTTAAGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCtttt                                                                                                                              <  2:4432266/150‑1 (MQ=255)
                           gACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAg                                                                                                                            <  1:3100936/150‑1 (MQ=255)
                            aCTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGc                                                                                                                           >  2:1830551/1‑150 (MQ=255)
                             cTTACGAAAAAGTTGAAATCATTAGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCt                                                                                                                          <  1:2286628/150‑1 (MQ=255)
                                 cGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGc                                                                                                                      <  1:4390871/150‑1 (MQ=255)
                                   aaaatGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTCTCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAAGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCaa                                                                                                                    >  1:1242240/1‑150 (MQ=255)
                                   aaaaaTTTTAAATCTTTCGTAGAAATGACGGTTCGATCCTCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACTAATGGGGTAACGATGTTACCACAAACGCAAAGCTTTATGGTAACCCCTTTTAGCTATGCaa                                                                                                                    <  2:551228/150‑1 (MQ=255)
                                   aaaaaGTTGAAATCATTCGTAGAAGTGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCaa                                                                                                                    >  2:3818043/1‑150 (MQ=255)
                                     aaaGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAAt                                                                                                                  <  2:1880982/150‑1 (MQ=255)
                                      aaGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATa                                                                                                                 <  1:1395249/150‑1 (MQ=255)
                                          tGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGTTAACCCCTTTTAGCTATGCAAATACGCa                                                                                                             >  2:3358069/1‑150 (MQ=255)
                                               tCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTa                                                                                                        <  2:4142416/150‑1 (MQ=255)
                                                cATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTaa                                                                                                       >  1:1259501/1‑150 (MQ=255)
                                                 aTTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCACATACGCAGTTTaaa                                                                                                      >  1:3307/1‑150 (MQ=255)
                                                 aTTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTaaa                                                                                                      >  1:3155/1‑150 (MQ=255)
                                                         aaaTGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTcc                                                                                              >  2:3012122/1‑150 (MQ=255)
                                                           aTGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGt                                                                                            >  2:724419/1‑150 (MQ=255)
                                                           aTGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGt                                                                                            >  2:1638847/1‑150 (MQ=255)
                                                           aTGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGt                                                                                            >  2:2031791/1‑150 (MQ=255)
                                                           aTGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGt                                                                                            >  2:3917804/1‑150 (MQ=255)
                                                                ggTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAAGCAAAGCTTAATGGTAACCCCTTTTAGCTCTGCAAATACGCAGTTTAAAAATGTTCCGTCAACt                                                                                       >  1:2806048/1‑150 (MQ=255)
                                                                  ttCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACttt                                                                                     >  2:316668/1‑150 (MQ=255)
                                                                  ttCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACttt                                                                                     <  2:1748422/150‑1 (MQ=255)
                                                                   tCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACtttt                                                                                    <  2:482880/150‑1 (MQ=255)
                                                                     gATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCa                                                                                  >  2:3421118/1‑150 (MQ=255)
                                                                           cTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTggg                                                                            <  1:1115833/150‑1 (MQ=255)
                                                                           cTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGAggg                                                                            >  2:673513/1‑150 (MQ=255)
                                                                              tAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTta                                                                         >  1:963884/1‑150 (MQ=255)
                                                                               aaCAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTtat                                                                        <  1:2628374/150‑1 (MQ=255)
                                                                                aCAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTtata                                                                       >  2:1487034/1‑150 (MQ=255)
                                                                                aCAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTtata                                                                       <  2:3108461/150‑1 (MQ=255)
                                                                                  aaTAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATc                                                                     >  1:1642126/1‑150 (MQ=255)
                                                                                   aTAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCa                                                                    >  1:218566/1‑150 (MQ=255)
                                                                                     aaGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAAt                                                                  >  2:3037341/1‑150 (MQ=255)
                                                                                     aaGATGCCTTATCACGGTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAAt                                                                  <  2:4182497/150‑1 (MQ=255)
                                                                                      aGATGCCTTATCACGTGACGCGGGGTTATATCGGAACAAATGGGGTAACGATGTGACCACAACCGCAAAGCTTACGGGTAACCCCTTTTAGCTAGGCAACTACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATg                                                                 <  2:1187307/150‑1 (MQ=255)
                                                                                        aTGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGAc                                                               <  1:534099/150‑1 (MQ=255)
                                                                                        aTGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGAc                                                               <  1:826929/150‑1 (MQ=255)
                                                                                                   cGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGAtt                                                    >  2:4221419/1‑150 (MQ=255)
                                                                                                    gTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTa                                                   >  2:1011350/1‑150 (MQ=255)
                                                                                                    gTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTa                                                   >  2:1010988/1‑150 (MQ=255)
                                                                                                      tACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAg                                                 >  2:1871591/1‑150 (MQ=255)
                                                                                                       aCGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGt                                                <  1:1830551/150‑1 (MQ=255)
                                                                                                       aCGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAAGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGt                                                <  1:2901900/150‑1 (MQ=255)
                                                                                                        cGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGtt                                               <  1:2030502/150‑1 (MQ=255)
                                                                                                          cTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGtttt                                             <  1:1890372/150‑1 (MQ=255)
                                                                                                           tGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGttttt                                            >  1:3338296/1‑150 (MQ=255)
                                                                                                           tGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGttttt                                            <  1:862266/150‑1 (MQ=255)
                                                                                                            gggTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTa                                           <  2:2163046/150‑1 (MQ=255)
                                                                                                            gggTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTa                                           <  2:4435524/150‑1 (MQ=255)
                                                                                                                 atatCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGaaa                                      >  2:4475201/1‑150 (MQ=255)
                                                                                                                     cGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATaaa                                  >  2:852059/1‑150 (MQ=255)
                                                                                                                     cGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATaaa                                  >  2:1788504/1‑150 (MQ=255)
                                                                                                                       gAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACt                                >  2:4354182/1‑150 (MQ=255)
                                                                                                                       gAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACt                                <  2:2818818/150‑1 (MQ=255)
                                                                                                                         aCAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGc                              <  2:3593733/150‑1 (MQ=255)
                                                                                                                              tGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTgt                         <  2:2803980/150‑1 (MQ=255)
                                                                                                                              tGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGGGGATTAAGTTTTTATGAAATAAACTGCCCTgt                         >  2:2544474/1‑150 (MQ=255)
                                                                                                                                  gTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTgtagta                     <  1:4068867/150‑1 (MQ=255)
                                                                                                                                   tAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTaa                    >  1:1869644/1‑150 (MQ=255)
                                                                                                                                     aCGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAAt                  <  2:2289699/150‑1 (MQ=255)
                                                                                                                                      cGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGATTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATa                 <  2:3847132/150‑1 (MQ=255)
                                                                                                                                       gATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATaa                <  2:1586130/150‑1 (MQ=255)
                                                                                                                                         tGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACa              >  1:1945997/1‑150 (MQ=255)
                                                                                                                                          gTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACAt             <  2:2940977/150‑1 (MQ=255)
                                                                                                                                             aCCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTg          >  1:416414/1‑150 (MQ=255)
                                                                                                                                              ccACAAAAGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGt         >  1:4537531/1‑150 (MQ=255)
                                                                                                                                                acaAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTaa       >  2:1557312/1‑150 (MQ=255)
                                                                                                                                                acaAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTaa       >  1:4205028/1‑150 (MQ=255)
                                                                                                                                                 caAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTAat      >  1:4099252/1‑150 (MQ=255)
                                                                                                                                                  aaaCGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTAata     >  2:66320/1‑150 (MQ=255)
                                                                                                                                                  aaaCGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTAata     <  2:905250/150‑1 (MQ=255)
                                                                                                                                                   aaCGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTAatat    <  1:102377/150‑1 (MQ=255)
                                                                                                                                                     cGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTAATAttt  <  2:4051078/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTATGAAATAAACTGCCCTGTAGTAAATAACATTTGTAATATTT  >  CP009273/3183668‑3183966

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: