Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 114,510 A→T 5.3% intergenic (+23/+12) guaC → / ← hofC GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913114,5100AT 5.3% 220.0 / 2.6 76intergenic (+23/+12)guaC/hofCGMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/31);  new base T (2/2);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

AGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCAT  >  NC_000913/114361‑114659
                                                                                                                                                     |                                                                                                                                                     
agaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCa                                                                                                                                                       <  2:1028465/150‑1 (MQ=255)
 gaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCAACTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCaa                                                                                                                                                      >  2:1119068/1‑150 (MQ=255)
       ttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGcc                                                                                                                                                <  2:582190/150‑1 (MQ=255)
         gggCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCAg                                                                                                                                              <  2:483929/150‑1 (MQ=255)
          ggcggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCAGc                                                                                                                                             <  2:1549025/150‑1 (MQ=255)
           gcggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCg                                                                                                                                            <  1:949490/150‑1 (MQ=255)
             ggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGtt                                                                                                                                          >  1:522406/1‑150 (MQ=255)
              gcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTa                                                                                                                                         >  1:1665343/1‑150 (MQ=255)
                cTGCGTTCAGCTTGTACATACGGTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCAGCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatc                                                                                                                                       <  1:389491/150‑1 (MQ=255)
                  gCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatccc                                                                                                                                     <  1:1542978/150‑1 (MQ=255)
                  gCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCAGAAAGAGCTGACCAAGCGCACCACGATTATTCCTGTGCAGGAACAAGAAAACCGCATCTTCAACAACATGTAATCTCCCAACGCTATCGTGGAGCAACACGACAGCTTTatccc                                                                                                                                     >  2:1408068/1‑150 (MQ=255)
                      tCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatcc                                                                                                                                 >  2:165023/1‑150 (MQ=255)
                       cAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatccc                                                                                                                                <  1:162141/150‑1 (MQ=255)
                       cAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGGGGAGCAACACGCCAGCGTTatcccatccc                                                                                                                                <  1:755589/150‑1 (MQ=255)
                         gCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATATCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCAc                                                                                                                              >  1:1436080/1‑150 (MQ=255)
                          cTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCCCGTTTATTCGTGTGCAGGAACAAGACAACCGCATCGTCAACCACCTCTAATCTCCCAACCCTGGCGTGGAGCAACACGCCAGCGTAAACCCATCCAACt                                                                                                                             >  1:1123213/1‑150 (MQ=255)
                          cTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACt                                                                                                                             <  1:1161349/150‑1 (MQ=255)
                               aCATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCACCAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACCCGCCAGCGATATCCCATCCCACTcatcg                                                                                                                        >  1:333320/1‑150 (MQ=255)
                               aCATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcg                                                                                                                        >  1:675674/1‑150 (MQ=255)
                                cATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgc                                                                                                                       >  1:161448/1‑150 (MQ=255)
                                cATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgc                                                                                                                       >  1:218475/1‑150 (MQ=255)
                                   aCGTTGGGGCTTCACGCCTGAAAGAGCTGACGAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcatc                                                                                                                    <  1:254070/150‑1 (MQ=255)
                                    cGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTTCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcatcg                                                                                                                   <  2:126630/150‑1 (MQ=255)
                                          ggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTaaa                                                                                                             >  1:1392035/1‑150 (MQ=255)
                                           gCTTCACGCCTGAAAGAGCTGCCCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCCCTCATCGCATCGCCTAAAt                                                                                                            >  2:1412102/1‑150 (MQ=255)
                                           gCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                            >  2:955770/1‑150 (MQ=255)
                                           gCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                            <  2:675674/150‑1 (MQ=255)
                                           gCTTCACGCCTGAAAGAGCTGACCAAGAGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTTGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                            <  1:94141/150‑1 (MQ=255)
                                            cTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTTATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCCTCCCACTCATCGCATCGCCTaaaag                                                                                                           >  1:245540/1‑148 (MQ=255)
                                            cTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATg                                                                                                           <  2:844037/150‑1 (MQ=255)
                                                 cGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaaaa                                                                                                      <  1:846924/150‑1 (MQ=255)
                                                   ccTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAAtt                                                                                                    <  1:1310219/150‑1 (MQ=255)
                                                   ccTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGAGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAAtt                                                                                                    <  1:1296957/150‑1 (MQ=255)
                                                      gAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGc                                                                                                 >  2:1375978/1‑150 (MQ=255)
                                                           agCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATa                                                                                            <  2:218475/150‑1 (MQ=255)
                                                            gCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATAc                                                                                           <  2:856185/150‑1 (MQ=255)
                                                               gACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACAtt                                                                                        >  2:1154307/1‑150 (MQ=255)
                                                                  cAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGcc                                                                                     >  2:1488950/1‑150 (MQ=255)
                                                                     gcgcACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccacc                                                                                  <  2:392315/150‑1 (MQ=255)
                                                                      cgcACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccacca                                                                                 >  2:1498497/1‑150 (MQ=255)
                                                                       gcACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccaccac                                                                                <  1:1481522/150‑1 (MQ=255)
                                                                            acGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCg                                                                           <  1:165023/150‑1 (MQ=255)
                                                                                      gtgtGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataat                                                                 >  2:6324/1‑150 (MQ=255)
                                                                                      gtgtGAAGTAACAAGAAAACCGAAGCTTCAACAACCTGTAATCTACCAACGCTGGCGTGGAGCAAAACGCCAGCGTTATACCATCCCAATCATCGCTTCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataat                                                                 <  1:1408068/150‑1 (MQ=255)
                                                                                         tgCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtcc                                                              <  1:1563469/150‑1 (MQ=255)
                                                                                         tgCACGACCAAGAAAACCGCATCTTCAACAACCTGTATTCTCCCAACGCTGGCGTGCAGCAACACGCCCGCGTTATCCCATCCCACTCATCGCCTCGCCTAAATGGAACATTGGCAGATTCATAGCCACCACCCGCGGTCCAATAATtcc                                                              >  2:1416985/1‑150 (MQ=255)
                                                                                             ggAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTccc                                                          >  2:1234160/1‑150 (MQ=255)
                                                                                             ggAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTACTccc                                                          >  1:1518594/1‑150 (MQ=255)
                                                                                              gAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCg                                                         <  1:1154307/150‑1 (MQ=255)
                                                                                               aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGt                                                        >  2:1187548/1‑150 (MQ=255)
                                                                                                 cAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCAACGCCAAAATGGAAAAATGGCACATACATTGCCACCACCAGCGAACCAAAAATTCCTCCCGtaa                                                      >  1:1217526/1‑148 (MQ=255)
                                                                                                   aGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATg                                                    <  2:1392035/150‑1 (MQ=255)
                                                                                                    gAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCCGATACATTGCCACCACCAGCGTCCCAATAATTCCTCCCGTTATGa                                                   >  1:338360/1‑150 (MQ=255)
                                                                                                        aCCGCAACTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAg                                               <  1:1119068/150‑1 (MQ=255)
                                                                                                                  caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCCCTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGt                                     <  2:411641/150‑1 (MQ=255)
                                                                                                                  caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGt                                     >  2:1294709/1‑150 (MQ=255)
                                                                                                                  caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGt                                     <  1:955770/150‑1 (MQ=255)
                                                                                                                   aacaacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTAGCAATAATTCCTCCCGTTATGATCAGCAACGCCGGtt                                    >  2:696030/1‑150 (MQ=255)
                                                                                                                    acaacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTc                                   >  2:1128544/1‑150 (MQ=255)
                                                                                                                        ccTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGt                               >  2:567684/1‑150 (MQ=255)
                                                                                                                          tGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTaa                             >  2:1672379/1‑150 (MQ=255)
                                                                                                                              aTCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCTTCGCCTAATTGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCt                         <  1:1156772/150‑1 (MQ=255)
                                                                                                                              aTCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCt                         <  1:1498497/150‑1 (MQ=255)
                                                                                                                              aTCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCt                         >  2:137666/1‑150 (MQ=255)
                                                                                                                                  cccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGa                     >  2:1626179/1‑150 (MQ=255)
                                                                                                                                  cccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGa                     <  2:161448/150‑1 (MQ=255)
                                                                                                                                    cAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCCTCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAgg                   >  1:150010/1‑150 (MQ=255)
                                                                                                                                       cGCTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTa                >  1:1410750/1‑150 (MQ=255)
                                                                                                                                       cGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTa                >  2:1517437/1‑150 (MQ=255)
                                                                                                                                        gCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTAt               >  1:1095170/1‑150 (MQ=255)
                                                                                                                                         cTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATc              >  1:531091/1‑150 (MQ=255)
                                                                                                                                         cTGGCGAGCAGCAACACGCCAGCGTTATCCCCTCCCACTCAACGCATAGCCTAAATGGAAATTTGGGAGATACATTGCCACAACCAGAGAACAAATAATTCCTCCCGTTTTGATCAGCAACTCCGGTTCCAGTAAAGCTTCGAGGTTATc              >  2:1130152/1‑150 (MQ=255)
                                                                                                                                           ggCGTGGAGCAACACGCCAGCGCTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCg            >  2:995063/1‑150 (MQ=255)
                                                                                                                                             cGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAgg                   <  2:627147/141‑1 (MQ=255)
                                                                                                                                             cGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAgg                   >  1:627147/1‑141 (MQ=255)
                                                                                                                                                 gAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAgcg      >  2:1684213/1‑150 (MQ=255)
                                                                                                                                                 gAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGCAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAgcg      >  2:1717985/1‑150 (MQ=255)
                                                                                                                                                    cAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCa   >  2:619239/1‑150 (MQ=255)
                                                                                                                                                     aaCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCAt  >  2:699706/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGCCAT  >  NC_000913/114361‑114659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑