Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,181,645 A→G 88.7% N83N (AAT→AAC yehA ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,181,6450AG88.7% 613.9 / 52.9 194N83N (AAT→AACyehAputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (14/8);  new base G (86/86);  total (100/94)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 200 total aligned reads displayed.

AATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTG  >  CP009273/2181498‑2181791
                                                                                                                                                   |                                                                                                                                                  
  tcatcaCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAagttta                                                                                                                                                <  2:1045463/150‑1 (MQ=255)
       aCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTTAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAAtttt                                                                                                                                           <  1:2067758/150‑1 (MQ=255)
       aCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAAtttt                                                                                                                                           <  2:2787292/150‑1 (MQ=255)
       aCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAAtttt                                                                                                                                           <  2:3011833/150‑1 (MQ=255)
       aCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAAtttt                                                                                                                                           <  2:618271/150‑1 (MQ=255)
       aCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAAtttt                                                                                                                                           >  2:181457/1‑150 (MQ=255)
         tGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTAt                                                                                                                                         <  1:3967208/150‑1 (MQ=255)
          gAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATg                                                                                                                                        >  2:2891834/1‑150 (MQ=255)
          gAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCAGTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTCTGGGTCCATTCAAAATGTAAGTTTAATTTTATg                                                                                                                                        >  1:2373346/1‑150 (MQ=255)
           aaCCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTAGTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATgg                                                                                                                                       >  1:150438/1‑150 (MQ=255)
            aCCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGt                                                                                                                                      <  2:823602/150‑1 (MQ=255)
              cgccgcCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTAt                                                                                                                                    <  2:393796/150‑1 (MQ=255)
              cgccgcCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTAt                                                                                                                                    >  2:2763321/1‑150 (MQ=255)
              cgccgcCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTAt                                                                                                                                    <  2:393830/150‑1 (MQ=255)
              cgccgcCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTAt                                                                                                                                    >  1:530656/1‑150 (MQ=255)
              cgccgcCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTAt                                                                                                                                    >  1:1851557/1‑150 (MQ=255)
                 cgccCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGt                                                                                                                                 >  2:841624/1‑150 (MQ=255)
                  gccCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGtt                                                                                                                                <  1:1554364/150‑1 (MQ=255)
                   ccccccGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTta                                                                                                                               <  2:801395/150‑1 (MQ=255)
                   ccccccGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTta                                                                                                                               <  2:2974864/150‑1 (MQ=255)
                        cGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAg                                                                                                                          <  1:3758598/150‑1 (MQ=255)
                           aGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGTTATCGTTATAAATTcc                                                                                                                       >  1:425895/1‑150 (MQ=255)
                           aGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTcc                                                                                                                       >  2:820053/1‑150 (MQ=255)
                           aGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTcc                                                                                                                       >  1:1971363/1‑150 (MQ=255)
                            gAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCca                                                                                                                      <  2:2192065/150‑1 (MQ=255)
                             aaaCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCcaa                                                                                                                     <  2:2224806/150‑1 (MQ=255)
                              aaCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCaaac                                                                                                                    >  2:3676967/1‑150 (MQ=255)
                               aCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCcaaca                                                                                                                   >  2:86420/1‑150 (MQ=255)
                                   tCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACaataa                                                                                                               >  2:1444864/1‑150 (MQ=255)
                                   tCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACaataa                                                                                                               >  2:55901/1‑150 (MQ=255)
                                     aCTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAataata                                                                                                             >  2:1836088/1‑150 (MQ=255)
                                      cTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCCCTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAatat                                                                                                            <  2:2782176/150‑1 (MQ=255)
                                      cTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAatat                                                                                                            <  2:2856587/150‑1 (MQ=255)
                                       ttttGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATc                                                                                                           >  1:1815887/1‑150 (MQ=255)
                                        tttGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATcc                                                                                                          <  2:1215986/150‑1 (MQ=255)
                                        tttGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATcc                                                                                                          >  1:3962417/1‑150 (MQ=255)
                                        tttGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATcc                                                                                                          <  2:3229228/150‑1 (MQ=255)
                                        tttGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATcc                                                                                                          <  1:3693960/150‑1 (MQ=255)
                                         ttGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTTAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCg                                                                                                         <  1:3643076/150‑1 (MQ=255)
                                                aCTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCt                                                                                                  <  1:643960/150‑1 (MQ=255)
                                                 cTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTc                                                                                                 <  2:4239759/150‑1 (MQ=255)
                                                 cTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTc                                                                                                 <  2:2565121/150‑1 (MQ=255)
                                                 cTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTc                                                                                                 >  2:1543932/1‑150 (MQ=255)
                                                   ggCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCAt                                                                                               >  1:278452/1‑150 (MQ=255)
                                                   ggCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCAt                                                                                               >  2:3119308/1‑150 (MQ=255)
                                                    gCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATa                                                                                              >  1:4166762/1‑150 (MQ=255)
                                                         aGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGt                                                                                         >  2:3920453/1‑150 (MQ=255)
                                                          ggTGTCACCGTAACTGAATAACCCCTGGTGAACCGGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTg                                                                                        <  2:1411083/150‑1 (MQ=255)
                                                            tgtCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGAt                                                                                      <  1:775994/150‑1 (MQ=255)
                                                            tgtCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGAt                                                                                      >  2:382545/1‑150 (MQ=255)
                                                                aCCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTa                                                                                  >  2:3604043/1‑150 (MQ=255)
                                                                aCCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTa                                                                                  >  1:4192983/1‑150 (MQ=255)
                                                                  cGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTAtt                                                                                >  2:1213968/1‑150 (MQ=255)
                                                                  cGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTAtt                                                                                >  2:2843343/1‑150 (MQ=255)
                                                                       cTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTAc                                                                           <  2:1349926/150‑1 (MQ=255)
                                                                           aTAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTAc                                                                       <  1:2707503/150‑1 (MQ=255)
                                                                            tAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACg                                                                      <  1:492647/150‑1 (MQ=255)
                                                                            tAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACg                                                                      <  2:3437838/150‑1 (MQ=255)
                                                                             aaCCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGc                                                                     >  1:756810/1‑150 (MQ=255)
                                                                              aCCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCa                                                                    <  1:2939945/150‑1 (MQ=255)
                                                                                 aCTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGa                                                                 <  1:1105358/150‑1 (MQ=255)
                                                                                      tGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAAt                                                            >  1:1173823/1‑150 (MQ=255)
                                                                                      tGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAAt                                                            <  2:450012/150‑1 (MQ=255)
                                                                                        aaaCTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCa                                                          >  2:3492228/1‑150 (MQ=255)
                                                                                        aaaCTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCa                                                          >  1:3846931/1‑150 (MQ=255)
                                                                                          aCTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCAtt                                                        <  1:1196251/150‑1 (MQ=255)
                                                                                          aCTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCAtt                                                        <  1:2188491/150‑1 (MQ=255)
                                                                                          aCTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCAtt                                                        <  2:2064825/150‑1 (MQ=255)
                                                                                                  gATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTAt                                                <  1:3378480/150‑1 (MQ=255)
                                                                                                   attattTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATg                                               <  1:2891834/150‑1 (MQ=255)
                                                                                                     tattTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCt                                             <  1:735320/150‑1 (MQ=255)
                                                                                                        ttGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCat                                          >  2:499182/1‑150 (MQ=255)
                                                                                                         tGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATGGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCata                                         >  1:2571035/1‑150 (MQ=255)
                                                                                                         tGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCata                                         >  1:2381227/1‑150 (MQ=255)
                                                                                                         tGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCata                                         >  1:2568870/1‑150 (MQ=255)
                                                                                                            cAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGt                                      <  2:2007690/150‑1 (MQ=255)
                                                                                                              aCGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCg                                    >  1:3327259/1‑150 (MQ=255)
                                                                                                                  gATGTTTTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAt                                <  2:3441983/150‑1 (MQ=255)
                                                                                                                  gATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAt                                >  2:1765252/1‑150 (MQ=255)
                                                                                                                  gATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAt                                >  1:2212215/1‑150 (MQ=255)
                                                                                                                   aTGTTGTTTTTATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAtt                               >  1:3245180/1‑150 (MQ=255)
                                                                                                                           tttATTGGTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTg                       <  1:247808/150‑1 (MQ=255)
                                                                                                                             tATTGGTCCATTCAAAATTTAAGTTTAATTTTATGGTATCGTTATAAAGTCCATCTATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAg                     <  2:1700336/150‑1 (MQ=255)
                                                                                                                             tATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAg                     >  2:3157234/1‑150 (MQ=255)
                                                                                                                                 ggTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAg                 <  1:2103346/150‑1 (MQ=255)
                                                                                                                                  gTCCATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGt                <  1:2858124/150‑1 (MQ=255)
                                                                                                                                    ccATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGc              >  1:135108/1‑150 (MQ=255)
                                                                                                                                     cATTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCa             <  2:4289461/150‑1 (MQ=255)
                                                                                                                                      aTTCAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCaa            >  2:2814613/1‑150 (MQ=255)
                                                                                                                                         cAAAATGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAAt         >  2:4018090/1‑150 (MQ=255)
                                                                                                                                          aaaaTGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATc        >  2:2650965/1‑150 (MQ=255)
                                                                                                                                          aaaaTGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATc        <  1:1319522/150‑1 (MQ=255)
                                                                                                                                          aaaaTGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATc        >  1:3022703/1‑150 (MQ=255)
                                                                                                                                          aaaaTGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATc        >  1:3419229/1‑150 (MQ=255)
                                                                                                                                          aaaaTGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATc        >  1:3420307/1‑150 (MQ=255)
                                                                                                                                           aaaTGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcc       >  2:1612310/1‑150 (MQ=255)
                                                                                                                                             aTGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccc     >  2:3953316/1‑150 (MQ=255)
                                                                                                                                             aTGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccc     >  2:3861603/1‑150 (MQ=255)
                                                                                                                                             aTGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccc     >  2:1840704/1‑150 (MQ=255)
                                                                                                                                             aTGTAAGTTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccc     >  1:3753845/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTG  >  CP009273/2181498‑2181791

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: