Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 114,506 G→T 5.1% intergenic (+19/+16) 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,5060GT 5.1% 236.7 / 3.3 79intergenic (+19/+16)guaC/hofCGMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base T (2/2);  total (43/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 = 3.02e-01

GCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCG  >  NC_000913/114358‑114655
                                                                                                                                                    |                                                                                                                                                     
gcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGa                                                                                                                                                      <  1:1312777/150‑1 (MQ=255)
 cgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAg                                                                                                                                                     >  2:1542978/1‑150 (MQ=255)
 cgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAg                                                                                                                                                     >  1:856185/1‑150 (MQ=255)
 cgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAg                                                                                                                                                     <  1:1405491/150‑1 (MQ=255)
   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                                 <  1:955770/150‑1 (MQ=255)
                                                                                                                     caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGt                                 >  2:1294709/1‑150 (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                     >  2:137666/1‑150 (MQ=255)
                                                                                                                                 aTCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCt                     <  1:1498497/150‑1 (MQ=255)
                                                                                                                                     cccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGa                 <  2:161448/150‑1 (MQ=255)
                                                                                                                                     cccAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGa                 >  2:1626179/1‑150 (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)
                                                                                                                                                    |                                                                                                                                                     
GCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCG  >  NC_000913/114358‑114655

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: