Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,259,987 A→G 5.4% T178A (ACC→GCC)  fimV → motility protein FimV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,259,9870AG 5.4% 223.3 / 2.2 74T178A (ACC→GCC) fimVmotility protein FimV
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base G (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCGCCCAGGGCGAGGCGAACCGCGAGGTGG  >  NC_002947/2259840‑2260128
                                                                                                                                                   |                                                                                                                                             
ctCGACCAGGCCAAGGCACAGGGCCCGCGCCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGAcc                                                                                                                                             <  1:414084/150‑1 (MQ=255)
 tCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCa                                                                                                                                            >  1:703270/1‑150 (MQ=255)
 tCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCa                                                                                                                                            <  1:344355/150‑1 (MQ=255)
   gACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATg                                                                                                                                          >  2:73609/1‑150 (MQ=255)
      cAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATc                                                                                                                                       <  2:1071436/150‑1 (MQ=255)
           cAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGAt                                                                                                                                  >  1:1286204/1‑150 (MQ=255)
           cAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGAt                                                                                                                                  >  1:527439/1‑150 (MQ=255)
           cAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTAATCTGGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGAt                                                                                                                                  <  1:836058/150‑1 (MQ=255)
            aaGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTACCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATc                                                                                                                                 >  1:1333840/1‑150 (MQ=255)
             aGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATcc                                                                                                                                <  1:222960/150‑1 (MQ=255)
               gCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCCCCGTCGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAg                                                                                                                              <  1:514639/150‑1 (MQ=255)
                cacaGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAgg                                                                                                                             <  2:1333840/150‑1 (MQ=255)
                 acaGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGc                                                                                                                            >  1:190196/1‑150 (MQ=255)
                    gggCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCttt                                                                                                                         <  2:99100/150‑1 (MQ=255)
                    gggCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCttt                                                                                                                         >  2:1309969/1‑150 (MQ=255)
                          aCAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAAccc                                                                                                                   >  1:921985/1‑150 (MQ=255)
                          aCAGCCGGCGGCAGGAAATGTCACAACGAACAAAACCGGTGCCCGCCGCCACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAAccc                                                                                                                   <  2:1163291/150‑1 (MQ=255)
                                ggcggcAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGc                                                                                                             <  1:198307/150‑1 (MQ=255)
                                  cggcAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACAACACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCt                                                                                                           <  2:703270/150‑1 (MQ=255)
                                    gcaggcaATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTc                                                                                                         <  2:1596462/150‑1 (MQ=255)
                                      aggcaATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCAt                                                                                                       >  2:725897/1‑150 (MQ=255)
                                      aggcaATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCAt                                                                                                       >  2:847415/1‑150 (MQ=255)
                                              tCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGcaaca                                                                                               >  2:395992/1‑150 (MQ=255)
                                              tCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGcaaca                                                                                               >  2:894954/1‑150 (MQ=255)
                                                acacCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCaacaac                                                                                             >  1:212137/1‑150 (MQ=255)
                                                acacCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCAACGTCaacgac                                                                                             >  2:1350687/1‑150 (MQ=255)
                                                acacCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCACGTTGCCGCGCTTAACACACAGGGAGAAGCGATCCAGCACAGCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGcaacaaa                                                                                             >  2:1418290/1‑149 (MQ=255)
                                                  acCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGGCCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACAt                                                                                           >  2:305318/1‑150 (MQ=255)
                                                  acCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGGCCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACAt                                                                                           >  2:872247/1‑150 (MQ=255)
                                                   ccGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCACCAGACCACAATCGCGTTCCAGTCTTTGATCCCGGATGCCTTCATCGGCAACAACATc                                                                                          <  1:1309969/150‑1 (MQ=255)
                                                      gCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCTTTCATCGGCAACAACATCAAc                                                                                       <  2:338431/150‑1 (MQ=255)
                                                          tCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGc                                                                                   <  2:308175/150‑1 (MQ=255)
                                                              cGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGaa                                                                               >  2:1646926/1‑150 (MQ=255)
                                                                   ccAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGGCCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTg                                                                          <  1:872247/150‑1 (MQ=255)
                                                                   ccAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGGCCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTg                                                                          <  1:305318/150‑1 (MQ=255)
                                                                   ccAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTg                                                                          >  1:1311441/1‑150 (MQ=255)
                                                                      gcagcTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCATGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTg                                                                          <  2:1040822/147‑1 (MQ=255)
                                                                      gcagcTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTg                                                                          >  1:1040822/1‑147 (MQ=255)
                                                                         gcTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAgg                                                                    >  2:310565/1‑150 (MQ=255)
                                                                          cTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGt                                                                   >  2:1031314/1‑150 (MQ=255)
                                                                          cTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGt                                                                   <  1:894954/150‑1 (MQ=255)
                                                                           tACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCg                                                                                                                  >  2:1046206/1‑102 (MQ=255)
                                                                           tACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCg                                                                                                                  <  1:1046206/102‑1 (MQ=255)
                                                                           tACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTa                                                                  >  2:1239810/1‑150 (MQ=255)
                                                                             caccacCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAAACAATGATAGCGATACAGCGTTTGGACCGGGATGCCTTCATCGGAAACAACAGCACCCAGCTGACGGTTGGCCAggaaat                                                                >  2:179802/1‑146 (MQ=255)
                                                                               ccaccCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACCgc                                                              <  1:1646926/150‑1 (MQ=255)
                                                                                caccCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCNGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTggg                                                             >  1:719608/1‑148 (MQ=255)
                                                                                 accCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTgcgc                                                            <  1:73609/150‑1 (MQ=255)
                                                                                 accCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTgcgc                                                            >  1:623146/1‑150 (MQ=255)
                                                                                  cccAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGc                                                                                                             >  1:870612/1‑100 (MQ=255)
                                                                                  cccAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGc                                                                                                             <  2:870612/100‑1 (MQ=255)
                                                                                    cAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTg                                                         >  1:316611/1‑150 (MQ=255)
                                                                                     aGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGc                                                        <  2:190196/150‑1 (MQ=255)
                                                                                             gACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACcag                                                <  2:921985/150‑1 (MQ=255)
                                                                                                   ttGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACcagcagcag                                          >  2:1488283/1‑150 (MQ=255)
                                                                                                       ggCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGCACACGTGTGCTTTTATCGTCTAGGACATCATCCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATcc                                      <  1:1433667/150‑1 (MQ=255)
                                                                                                       ggCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGag                                                                               >  2:1121760/1‑108 (MQ=255)
                                                                                                       ggCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGaa                                                                               <  1:1121760/109‑1 (MQ=255)
                                                                                                        gCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCa                                     <  1:847415/150‑1 (MQ=255)
                                                                                                             tgtgcccgcgTTACCCCCCGGGTGAATCGATTCAACGGACCTTTATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCa                                <  2:719608/144‑1 (MQ=255)
                                                                                                                ccGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCg                             >  1:1644597/1‑150 (MQ=255)
                                                                                                                      gTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCGCCCAgg                       <  2:316611/150‑1 (MQ=255)
                                                                                                                       tAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCGCCCAggg                      <  1:310565/150‑1 (MQ=255)
                                                                                                                            cccAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCGCCCAGggcgagg                 <  2:527439/150‑1 (MQ=255)
                                                                                                                                   tGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCGCCCAGGGCGAGGCGAACcg          >  1:1116814/1‑150 (MQ=255)
                                                                                                                                     gatcgatcCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCag                                    >  1:189644/1‑122 (MQ=255)
                                                                                                                                     gatcgatcCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCag                                    <  2:189644/122‑1 (MQ=255)
                                                                                                                                          atcCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATc                              >  2:424054/1‑123 (MQ=255)
                                                                                                                                          atcCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATc                              <  1:424054/123‑1 (MQ=255)
                                                                                                                                           tcCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCGCCCAGGGCGAGGCGAACCGCGAGGTgg  >  2:1448082/1‑150 (MQ=255)
                                                                                                                                              agcagACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAg                                                                                    <  2:330480/65‑1 (MQ=255)
                                                                                                                                              agcagACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAg                                                                                    >  1:330480/1‑65 (MQ=255)
                                                                                                                                              agcagACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGccc                                                      <  1:1451899/95‑1 (MQ=255)
                                                                                                                                              agcagACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGccc                                                      >  2:1451899/1‑95 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATTGCCGCGCGTAACACCCAGGGTGGATCGATCCAGCAGACCATGATCGCGATCCAGGCTTTGAACCCGGATGCCTTCATCGGCAACAACATCAACCAGCTGAAGGTTGGCCAGGTACTGCGCCTGCCCGACCAGCAGCAGATCCAGAGCATCGCCCAGGGCGAGGCGAACCGCGAGGTGG  >  NC_002947/2259840‑2260128

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