Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,537 G→C 5.6% A216G (GCC→GGC)  yjhP ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,5370GC 5.6% 385.6 / 11.2 107A216G (GCC→GGC) yjhPputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (53/48);  new base C (4/2);  total (57/50)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 115 total aligned reads displayed.

GTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCT  >  NC_000913/4532393‑4532680
                                                                                                                                                |                                                                                                                                               
gTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              >  3:498682/1‑148 (MQ=255)
            tatTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              >  3:285970/1‑136 (MQ=255)
            tatTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              <  4:285970/136‑1 (MQ=255)
                         aaaTACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccag                                                                                                                     >  2:140273/1‑148 (MQ=255)
                               cGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              <  2:217872/117‑1 (MQ=255)
                               cGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              >  1:737757/1‑117 (MQ=255)
                               cGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              <  2:737757/117‑1 (MQ=255)
                                gggAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcgg                                                                                                              >  3:1031331/1‑148 (MQ=255)
                                gggAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcgg                                                                                                              >  3:394881/1‑148 (MQ=255)
                                       tgttgtTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATgg                                                                                                       >  1:854412/1‑148 (MQ=255)
                                         ttgttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCAt                                                                                                                                        >  1:703341/1‑113 (MQ=255)
                                         ttgttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCAt                                                                                                                                        <  2:703341/113‑1 (MQ=255)
                                         ttgttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagc                                                                                                                    >  3:405532/1‑133 (MQ=255)
                                         ttgttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagc                                                                                                                    <  4:405532/133‑1 (MQ=255)
                                          tgttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCa                                                                                                    <  4:394881/148‑1 (MQ=255)
                                           gttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAg                                                                                                   <  4:1031331/148‑1 (MQ=255)
                                            ttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              >  2:441932/1‑104 (MQ=255)
                                            ttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              <  1:441932/104‑1 (MQ=255)
                                                cGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                                                                 >  1:640075/1‑97 (MQ=255)
                                                cGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                                                                 <  2:640075/97‑1 (MQ=255)
                                                    cGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                                                                         <  3:1104593/133‑1 (MQ=255)
                                                    cGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                                                                         >  4:1104593/1‑133 (MQ=255)
                                                       tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                     <  1:63186/102‑1 (MQ=255)
                                                       tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                     >  2:63186/1‑102 (MQ=255)
                                                       tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGc                                                                                       <  1:140273/148‑1 (MQ=255)
                                                               aaCACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAcc                                                                               >  4:115723/1‑148 (MQ=255)
                                                                   cGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTc                                                                                                                                          <  1:158043/85‑1 (MQ=255)
                                                                   cGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTc                                                                                                                                          >  2:158043/1‑85 (MQ=255)
                                                                          gCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                                     <  3:1033605/131‑1 (MQ=255)
                                                                            cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGcc                                                                     <  4:438137/145‑1 (MQ=255)
                                                                            cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGcc                                                                     >  3:438137/1‑145 (MQ=255)
                                                                             aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              >  1:873187/1‑71 (MQ=255)
                                                                             aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                                                              <  2:873187/71‑1 (MQ=255)
                                                                              aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                                <  3:309758/148‑1 (MQ=255)
                                                                              aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                                <  2:854412/148‑1 (MQ=255)
                                                                              aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                                >  2:1075962/1‑148 (MQ=255)
                                                                              aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                                >  3:82764/1‑148 (MQ=255)
                                                                                 cATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                                     >  3:50833/1‑124 (MQ=255)
                                                                                 cATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                                     <  4:50833/124‑1 (MQ=255)
                                                                                 cATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt                                                             <  4:746514/148‑1 (MQ=255)
                                                                                   ttCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCgg                                                                                        >  2:457028/1‑119 (MQ=255)
                                                                                   ttCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCgg                                                                                        <  1:457028/119‑1 (MQ=255)
                                                                                    tCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTg                                                          >  3:1090606/1‑148 (MQ=255)
                                                                                      ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATg                                                                                                        <  4:180327/100‑1 (MQ=255)
                                                                                      ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATg                                                                                                        >  3:180327/1‑100 (MQ=255)
                                                                                          gcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTa                                                    <  4:82764/148‑1 (MQ=255)
                                                                                             tAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCg                                                                                         >  4:945089/1‑108 (MQ=255)
                                                                                             tAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCg                                                                                         <  3:945089/108‑1 (MQ=255)
                                                                                               ggTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAtt                                               >  1:559166/1‑148 (MQ=255)
                                                                                                 tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                                                                    >  4:376375/1‑61 (MQ=255)
                                                                                                 tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                                                                    <  1:879583/61‑1 (MQ=255)
                                                                                                 tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                                                                    >  2:879583/1‑61 (MQ=255)
                                                                                                 tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                                                                    <  3:376375/61‑1 (MQ=255)
                                                                                                 tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTc                                             <  1:1075962/148‑1 (MQ=255)
                                                                                                  cACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCtt                                                                   >  3:977392/1‑125 (MQ=255)
                                                                                                  cACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTcc                                                                 >  4:856278/1‑127 (MQ=255)
                                                                                                  cACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTcc                                                                 <  3:856278/127‑1 (MQ=255)
                                                                                                   aCGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGCACTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGGACCATTTcca                                           <  4:100874/148‑1 (MQ=255)
                                                                                                      tAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAt                                                >  2:955188/1‑140 (MQ=255)
                                                                                                      tAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAt                                                <  1:766556/140‑1 (MQ=255)
                                                                                                      tAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAt                                                >  2:766556/1‑140 (MQ=255)
                                                                                                      tAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAt                                                <  1:955188/140‑1 (MQ=255)
                                                                                                       aGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTg                                                                                           >  2:91401/1‑96 (MQ=255)
                                                                                                       aGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTg                                                                                           <  1:91401/96‑1 (MQ=255)
                                                                                                         cGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                                                       >  1:1073291/1‑66 (MQ=255)
                                                                                                           tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTc                                             >  1:879361/1‑138 (MQ=255)
                                                                                                           tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTc                                             <  2:879361/138‑1 (MQ=255)
                                                                                                            ttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCAt                                  <  4:1090606/148‑1 (MQ=255)
                                                                                                               ggCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGc                               >  2:49720/1‑148 (MQ=255)
                                                                                                                 cgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                                                                   <  1:655763/94‑1 (MQ=255)
                                                                                                                 cgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                                                                   >  2:655763/1‑94 (MQ=255)
                                                                                                                            gTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGaa                  <  2:559166/148‑1 (MQ=255)
                                                                                                                             ttCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTAcc                                                  <  2:410584/115‑1 (MQ=255)
                                                                                                                             ttCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTAcc                                                  >  1:410584/1‑115 (MQ=255)
                                                                                                                             ttCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAAc                 >  1:461359/1‑148 (MQ=255)
                                                                                                                              tCGGCCCTGACTTCGGCGGCGAGGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCCGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACg                >  2:348988/1‑148 (MQ=255)
                                                                                                                              tCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcgg                                                                                                              >  4:1085324/1‑54 (MQ=255)
                                                                                                                              tCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcgg                                                                                                              <  3:1085324/54‑1 (MQ=255)
                                                                                                                                    cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTg                                                                                           <  2:1124615/67‑1 (MQ=255)
                                                                                                                                    cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTg                                                                                           >  1:1124615/1‑67 (MQ=255)
                                                                                                                                    cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCTACCACGTCATAGCCg                             <  4:256836/129‑1 (MQ=255)
                                                                                                                                    cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCg                             >  3:256836/1‑129 (MQ=255)
                                                                                                                                      gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTg                                                          >  2:528736/1‑98 (MQ=255)
                                                                                                                                      gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTg                                                          <  1:114059/98‑1 (MQ=255)
                                                                                                                                      gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTg                                                          >  2:114059/1‑98 (MQ=255)
                                                                                                                                       aCTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCAt                                  >  4:401649/1‑121 (MQ=255)
                                                                                                                                       aCTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACaa       >  3:468342/1‑148 (MQ=255)
                                                                                                                                       aCTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACaa       >  3:1045990/1‑148 (MQ=255)
                                                                                                                                        cTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcgg                                                                                                              <  2:446412/44‑1 (MQ=38)
                                                                                                                                        cTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcgg                                                                                                              >  1:446412/1‑44 (MQ=38)
                                                                                                                                         ttCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAttt                                              >  4:436608/1‑107 (MQ=255)
                                                                                                                                         ttCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAttt                                              <  3:436608/107‑1 (MQ=255)
                                                                                                                                          tCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATa                                                                                 <  1:151994/71‑1 (MQ=255)
                                                                                                                                          tCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATa                                                                                 >  2:151994/1‑71 (MQ=255)
                                                                                                                                          tCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTACCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTc    >  3:731948/1‑148 (MQ=255)
                                                                                                                                           cggcggCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTcc   <  1:519957/148‑1 (MQ=255)
                                                                                                                                            ggcggcGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCt  <  4:731948/148‑1 (MQ=255)
                                                                                                                                             gcggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCgag                           <  4:391275/122‑1 (MQ=255)
                                                                                                                                             gcggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCgag                           >  3:391275/1‑122 (MQ=255)
                                                                                                                                              cggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                                                               >  1:234855/1‑85 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCT  >  NC_000913/4532393‑4532680

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: