Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,522,373 A→G 14.5% L126P (CTC→CCC)  hofM ← DNA catabolic putative pilus assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,522,3730AG14.5% 165.0 / 2.3 71L126P (CTC→CCC) hofMDNA catabolic putative pilus assembly protein
Reads supporting (aligned to +/- strand):  ref base A (11/47);  new base G (9/3);  total (21/50)
Fisher's exact test for biased strand distribution p-value = 3.32e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATCCGCCC  >  NC_000913/3522137‑3522611
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
ggcgaggcATGCCCACCACCAGTTTGCGCCCCCCGCTAAAGGCCGTCGGCAAACGCCACTGTTCGTTATCACGCAAGGCCAGACATTGCTGATGAAAAGGTAAAAAAGGCAGGAATCGCTGTAAGGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtgagt                                                                                                                                                                                                                                                     <  2:852028/232‑1 (MQ=255)
                  cccttttGCGCCCCCCGGCCTTGGCGGGGCGCCCCCAACCACTGGTCGTTATCACGCCAGGCCAGACATTGCTGATGGGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATc                                                                                                                                                                                                                                  <  1:959896/235‑1 (MQ=255)
                   cAGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGACGGGACTGGGGGAGTCTTCCGCTGAATCGa                                                                                                                                                                                                                                >  2:267282/1‑239 (MQ=255)
                    aGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGCCTGGGGAGTCTTTGCTGTAATCGaa                                                                                                                                                                                                                               >  2:245949/1‑239 (MQ=255)
                          gcccccccAGGTATAGGGCGTCGCCCACAGCCACTGTTCGTTAGCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCa                                                                                                                                                                                                                         <  2:731208/236‑1 (MQ=255)
                          cgcCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGGGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGAGTGAGTGAGTCTTCGCTGTACTCGAAGCGCac                                                                                                                                                                                                                        >  2:883803/1‑238 (MQ=255)
                           gcCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAgg                                                                                                                                                                                                                       <  2:377804/239‑1 (MQ=255)
                                ccagcttaaggcggTCGCCCAAAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTTATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGt                                                                                                                                                                                                                   <  2:626964/226‑1 (MQ=255)
                                 cAGCTATAGCGGGCGGACCACAGCTCATGTTCGTTGTCACGCCAGGCCAGACATTGCTGGTGAGAAGGTAAAGAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGAT‑‑CGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCGTTGCTGTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCg                                                                                                                                                                                                                <  2:236068/239‑1 (MQ=255)
                                 cAGCTATAGCGGCTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTc                                                                                                                                                                                                                  <  2:685858/239‑1 (MQ=255)
                                 cAGCTACAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCCGACATTGCCGATGAGAAGGTAAAAAAGGCAGGAATCGATGTAATGCACTGGCATCCGGTGGGATCGCGCTCACACAAACACGCACCCTTTCGGCCCTAGTAACCAGCTTTGCCAGCCCTTTACGTTGCCCGGCAGCCACTATACTGACCGAACTGCGACAGCGTTCGATGTAATCGAAacggaacaggga                                                                                                                                                                                                                  >  1:394041/1‑230 (MQ=255)
                                   gttctaGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGGCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCg                                                                                                                                                                                                                <  1:245949/235‑1 (MQ=255)
                                   ggtataGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCg                                                                                                                                                                                                                <  2:774552/237‑1 (MQ=255)
                                   gctcaaggggctCGCCCAAAGCCACTGTTCGTTATCCCGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCCGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCg                                                                                                                                                                                                                <  2:407482/228‑1 (MQ=255)
                                   gcctatgccggTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCTCTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCg                                                                                                                                                                                                                <  2:555877/229‑1 (MQ=255)
                                   gCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCg                                                                                                                                                                                                                <  2:604054/239‑1 (MQ=255)
                                       taGCGCGTCGGCACCAGCCACTGTTCGTTATCGCGCCAGGGCAGACATTGCTGATGAGAAGGTAAAAAAGGAAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATc                                                                                                                                                                                                            <  2:667918/239‑1 (MQ=255)
                                        aGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGGGAGTTTTCGCCGTAATccatggcccggggaaccggacc                                                                                                                                                                                                             >  2:1144541/1‑217 (MQ=255)
                                            cgTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGACGAGTGAGTCTTCTCTGTATTCGGAgggccaggggctcgggtccaaa                                                                                                                                                                                                        >  2:1225075/1‑220 (MQ=255)
                                              tCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGGGTAAGCAGCGTTGCCAGCTCTTTGCTTTCCGCGGCAGTCCCGTTTTAGCCGGGCCTAGGTGGCGTCTCCCTGTAATCCAAGCCCCGGGAGACCggggcccaaacc                                                                                                                                                                                                     >  1:592899/1‑229 (MQ=255)
                                              tCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGGAGGAGTCTTCGCCTTAATCCGACGGCAGGGGGTCCGGATcccaaatc                                                                                                                                                                                                     >  2:116714/1‑233 (MQ=255)
                                              tCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGGGAGGCTTCCCTGTAATCGAAGCGCAGGGAGGCCGGAtgccaatcc                                                                                                                                                                                                     >  2:175437/1‑231 (MQ=255)
                                                           cggTTCGTTAACACGCCCGGCCAGACATTGCTGAGGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCa                                                                                                                                                                                        <  2:875270/237‑1 (MQ=255)
                                                           cTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCa                                                                                                                                                                                        <  2:733016/239‑1 (MQ=255)
                                                                cGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGAGCGGGGGGAGTCTTTCCTGGTATAGCAGCGGCAGGGGGCCCGGACCCaatcccactccccggccccctc                                                                                                                                                                                     >  1:1024002/1‑217 (MQ=255)
                                                                cGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTAACGTTATAGGCGGGCCAGGGGGAGTTCTCTCTGTAATCGAAGCCCCGGGGGACCGGAACCCCATCCaaccccgggcgccacct                                                                                                                                                                                    >  2:666593/1‑222 (MQ=255)
                                                                      cACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCCTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGCATGAGGGAGTCTTTTCTTTTATCCAAACGCAGCGGGTCCGGGACCCTATCCCGCTCCCCGGCCCACGTCCCtca                                                                                                                                                                              >  2:712081/1‑236 (MQ=255)
                                                                        cGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGCCTGGGTGAGTCTTCCGTGTAATCTGAAAGCCAGGAGACCCGTTCCATATCCAACTCGCCGGCCCCCGGCCCTcaca                                                                                                                                                                           >  2:573976/1‑238 (MQ=255)
                                                                        cGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTTGCTGTAATCGTAACGCAGGGGGGCCGGAACCCTATCCCGCTCCCCGGCCATCCTCCCttcca                                                                                                                                                                            >  2:871081/1‑234 (MQ=255)
                                                                                  gACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTGCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTc                                                                                                                                                                                                <  2:869939/208‑1 (MQ=255)
                                                                                      ttGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCt                                                                                                                                                           <  2:1307057/239‑1 (MQ=255)
                                                                                       tGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGGGGGAGTTTTCGCTGTCATTGAAGGGCAGGGAGTCCGGGTCCCTATCCAGCTCCCGGGCCCACGTCCCTGACAGCCCAGCCGGTTTCt                                                                                                                                                           >  2:392862/1‑238 (MQ=255)
                                                                                              gagaAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGCTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGGTCCATATCCCGCCTGCCGGCCCTCCTCCCCGACACCCAGTCCGGTTGCTccccccc                                                                                                                                                    >  2:573283/1‑234 (MQ=255)
                                                                                                      aaaaaaGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTAACATGAACCCGCATCCTTTCTGCCAGAGTAAGCCGCGTTGCCAGCTCGTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAggg                                                                                                                                            <  2:517791/238‑1 (MQ=255)
                                                                                                              gatgtgccctcggaaTGCACTGGCCCACTGGGTGATCGCGGACCATGGATCCAGGAACCTTTCTGCCAGAGTAAGCGCGTTGTCGAGCTCTTTGATGTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCGCGCGGGCCATCGTCCCGGACTGCCGGGCCGTTGGCTCCCGCTCTCCAACGGACATCGAc                                                                                                                                   <  1:267282/225‑1 (MQ=255)
                                                                                                              cAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGCTGAGGGGGTCTTCGCTGTAATCGAAGCGCCGGGAGTCCGGATCCATATCCCGCTCCCGGGCCCTCGTCCCTGCACACCAGGCCGCTTGTTCCCCCTCACCCCGGGGCATCGACg                                                                                                                                  >  2:668786/1‑239 (MQ=255)
                                                                                                              cAGGAATCGCTGGAATGCACGGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGa                                                                                                                                    <  2:542600/237‑1 (MQ=255)
                                                                                                                ggAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACgg                                                                                                                                 <  2:263262/238‑1 (MQ=255)
                                                                                                                         gTAATGCTCTGGCAGCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAAt                                                                                                                       <  2:1238275/239‑1 (MQ=255)
                                                                                                                           aaTGCACTGGCATCCGGGTGGATCGCGCTCACATGAACCCGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGATTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGa                                                                                                                     <  2:59254/239‑1 (MQ=255)
                                                                                                                            aTGCAATGGCATCCGGGGTGATCGCGCTCACATGAACCCGCAACCTTTCTGACAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGAc                                                                                                                    <  2:790522/239‑1 (MQ=255)
                                                                                                                                    gCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGTCTGAGTGAGTCTTTCCTGTCATCGAAGCGGACGGGGTGCGGAGCCCTATTCAGCTCGCGGGCCATCGTCCCCTCCAGCCCCGCCGTTTTGCCCCCCCCCCCAAAGGAACACCCAGGGGCGGGGaaaaaacccgtgaa                                                                                                             >  1:1142378/1‑228 (MQ=255)
                                                                                                                                            ggttaccgcgCTCACATTAACACGCAACTTTTCGGCCAGAGGAAGCAGCGTTGCCAGCGCTGTGCTTTGCGCGGCAGTCTCGTTATAGGCGGGACTGAGTGAGTCTTCGCGGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGa                                                                                                     <  1:668786/232‑1 (MQ=255)
                                                                                                                                             ggtgagcgggcctactggAACACGCAACGTTTAGCCAAGGGGGAGCAGTGTTGCCGACGCTTTGCTTTGCGCGGCAGGCACGGTATAGGCGGGACGGAGGGAGTCTCCGCTGTAAGCGAAGCGCAGGGAGTCCGGATCCATATACAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGa                                                                                                     <  1:1144541/221‑1 (MQ=255)
                                                                                                                                                  cgcgCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGGTCCATATCCAGCTCGCGGGCCATCCTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCTACGGGGGCGGGAAAGACCGCTGTAATGTGCGACTGGCgg                                                                                             >  2:635952/1‑239 (MQ=255)
                                                                                                                                                    cgcTCAGATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGa                                                                                           <  2:1142378/239‑1 (MQ=255)
                                                                                                                                                    cgcTCACAAGAACCCGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAAGCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGa                                                                                           <  2:1317461/239‑1 (MQ=255)
                                                                                                                                                     cgtccccttggaccgaccaCCTTTCTGCCAAGGGAAGCAGGGGTGCCAGCTCTTGTCTGTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGa                                                                                           <  1:175437/221‑1 (MQ=255)
                                                                                                                                                      cTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCTCTGTAATCGAAGCGCAGGGGGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCCAGACGTTTGCGCCCGCTCACCAAGGGGCAACGCCGGGCGCGGGAATGACCGCTGGAAAGTGGGACTCGCggggga                                                                                         >  2:1325781/1‑236 (MQ=255)
                                                                                                                                                        cacaTGAACCCGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACg                                                                                       <  2:388255/239‑1 (MQ=255)
                                                                                                                                                        cacaTGAACACGCAACCTTTCGGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACg                                                                                       <  2:1321966/239‑1 (MQ=255)
                                                                                                                                                         acaTGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTGGCCGTCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGc                                                                                      <  2:458108/239‑1 (MQ=255)
                                                                                                                                                           aTGAACCCGCAACCTTTCTGCCAGGGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCa                                                                                    <  2:592899/239‑1 (MQ=255)
                                                                                                                                                              gAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTTTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAAc                                                                                  <  2:639641/239‑1 (MQ=255)
                                                                                                                                                                  ccGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCTGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGTTGTAATGTGCGACTGGCGGGAAACGCCAACATaa                                                                              <  2:1093469/238‑1 (MQ=255)
                                                                                                                                                                    gCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTTGGGAAGTTCTATCTGTAAACGAAGCGCAGGGGATGCGGGTCCCTATCCCACCGCCGGGGCACCTCCCCCCGCAGACAGGGCGGTAGCTCCCCCCCCCCACAGGGCATACGACGGGCGCGAAAAGAACGATGTGAAGTGCGACGCGCCGGGAAAGGCCaaaaaaaa                                                                             >  1:464446/1‑233 (MQ=255)
                                                                                                                                                                        ccTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCGCGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGtgat                                                                        <  2:1305787/239‑1 (MQ=255)
                                                                                                                                                                                           agcGTTGCCAGCTTGTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTcgc                                                      <  2:701994/238‑1 (MQ=255)
                                                                                                                                                                                             cGTTGCCAGCTCTTTGCTTTGCCGGGCAGTCCCGGTATAGGCGGGACTGAGTTAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGa                                                    <  2:576191/238‑1 (MQ=255)
                                                                                                                                                                                                tGCCAGCTCTTTGCTGTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTGCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTc                                                 <  2:875371/238‑1 (MQ=255)
                                                                                                                                                                                                 tgcggtctttttggttggcGCGACAGTATCGTTATACGGGGGCGTGGCTGAGTCCTCTCTGTAATGGAGACGAGGGGAGGCCGGATCCATATCCAGCTCGCGGGCCCTCGTCCCTGCCAGACAGTCCGTTTGCGCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGAGTCTGCGGCAGTGCGCGAcgc                                                 <  1:573283/229‑3 (MQ=255)
                                                                                                                                                                                                  tgccggcttttgttttggCGCGGAAGTCCCTTTACAGGCGGGACTGGGTGGGTCCTGGGTTGACTCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGGCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTcc                                                <  1:883803/231‑1 (MQ=255)
                                                                                                                                                                                                   cccggtctttggcggtgggcgGGAGTCACGGTATAGGCGGAATTTAGTGGGTCTTCGCTGGAATCGAAGCGGAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCTTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACAGCGACGGGCGCGGAAATGACCGCTTTAATGTGCGACTGGCGGGAAACGCCAACATAATGGGATGACGCTGCGGCAGTTCGCGACTCCa                                               <  1:1280389/221‑1 (MQ=255)
                                                                                                                                                                                                       ctctTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGt                                                                                                                                                                                    <  2:456906/105‑1 (MQ=255)
                                                                                                                                                                                                          tttgcttggcGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTaaca                                       <  2:495105/239‑1 (MQ=255)
                                                                                                                                                                                                          gtggctttggCCGGCAGTCACGTTATAGGGGGGACTGAGTTAGTCTTCGCTGGAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTaaca                                       <  2:57313/236‑1 (MQ=255)
                                                                                                                                                                                                           tgtgcgtggcgcGGGAGGCACGTTATAGGCGGGACTGAGTGAGTCTGCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCGCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTaac                                        <  2:1158053/230‑1 (MQ=255)
                                                                                                                                                                                                                         cAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAg                          <  2:8272/238‑1 (MQ=255)
                                                                                                                                                                                                                            tCACGTTATAGGCCGGACTGACGGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATAGCCAGCGCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAgctgc                      <  2:403175/239‑1 (MQ=255)
                                                                                                                                                                                                                                gTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGctgctg                    <  2:1008528/237‑1 (MQ=255)
                                                                                                                                                                                                                                  tataGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGAcc                                                                                                                   <  2:38713/141‑1 (MQ=255)
                                                                                                                                                                                                                                         ggggCTGGGTGAGTCTTCGGTGTAATCGGAGCGGAGGGGGTCCGGACCCATATCCCGCTCGCGGGCCATCGTCCCGGACAGCCAGGCCGGTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATcaaca          <  1:116714/238‑1 (MQ=255)
                                                                                                                                                                                                                                                agtgagtCTTCGCTGTAACCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATCCGccc  <  2:1024002/239‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
CGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATCCGCCC  >  NC_000913/3522137‑3522611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: