Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,616,606 G→C 5.2% A171A (GCG→GCC bepA → OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,616,6060GC 5.2% 384.9 / 7.9 112A171A (GCG→GCCbepAOM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (53/53);  new base C (3/3);  total (56/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 112 total aligned reads displayed.

TTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGC  >  NC_000913/2616459‑2616751
                                                                                                                                                   |                                                                                                                                                 
ttATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGa                                                                                                                                                  >  2:18913/1‑149 (MQ=255)
ttATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGa                                                                                                                                                  <  3:655757/149‑1 (MQ=255)
 tATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGAt                                                                                                                                                 <  1:894044/149‑1 (MQ=255)
  aTTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                                                                <  3:100368/149‑1 (MQ=255)
    tCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGTCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCGGTCCTTGGTGGGTGTATGAGGTTTTAGTTTACGGGCGATGGc                                                                                                                                              <  4:558236/149‑1 (MQ=255)
    tCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                              <  4:823001/149‑1 (MQ=255)
       gATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAg                                                                                                                                           >  2:643524/1‑149 (MQ=255)
            cGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGc                                                                                                                                      >  2:159442/1‑149 (MQ=255)
            cGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                      <  1:643524/149‑5 (MQ=255)
              aaaGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCGACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaata                                                                                                                                    <  3:889498/149‑7 (MQ=255)
                       tGGCTTCAGTTATGGCGCACGAAATCTCTCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATg                                                                                                                           <  3:86993/149‑1 (MQ=255)
                       tGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATg                                                                                                                           >  1:812314/1‑149 (MQ=255)
                          cTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcg                                                                                                                        >  3:649534/1‑149 (MQ=255)
                          cTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcc                                                                                                                        >  4:1058142/1‑148 (MQ=255)
                           ttCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                                                                >  4:1004187/1‑124 (MQ=255)
                           ttCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                                                                <  3:1004187/124‑1 (MQ=255)
                                     gcgcACGAAATCTCCCACGTCACCCAACGCCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCgg                                                                                                             >  4:628632/1‑149 (MQ=255)
                                      cgcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGt                                                                                                            <  1:725175/149‑1 (MQ=255)
                                       gcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTa                                                                                                           >  3:369256/1‑149 (MQ=255)
                                         aCGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTaca                                                                                                         >  1:258592/1‑149 (MQ=255)
                                              aTCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGc                                                                                                    <  1:980942/149‑1 (MQ=255)
                                               tctcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGAc                                                                                                                <  1:759904/136‑1 (MQ=255)
                                               tctcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGAc                                                                                                                >  2:759904/1‑136 (MQ=255)
                                               tctcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCg                                                                                                   >  1:164795/1‑149 (MQ=255)
                                               tctcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCg                                                                                                   >  1:373867/1‑149 (MQ=255)
                                                 tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCggg                                                                                                 >  4:723059/1‑149 (MQ=255)
                                                 tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCggg                                                                                                 >  2:112066/1‑149 (MQ=255)
                                                      cGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAcgc                                                                                            <  4:666416/149‑1 (MQ=255)
                                                      cGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAcgc                                                                                            <  2:1068064/149‑1 (MQ=255)
                                                          aCCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCa                                                                                        >  1:967079/1‑149 (MQ=255)
                                                            ccAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGTCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgg                                                                                      >  4:741058/1‑149 (MQ=255)
                                                            ccAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgg                                                                                      >  4:834309/1‑149 (MQ=255)
                                                            ccAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgg                                                                                      <  4:1069904/149‑1 (MQ=255)
                                                              aaCGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgggg                                                                                    <  1:161167/149‑1 (MQ=255)
                                                                cGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgat                                                                                  >  4:643847/1‑149 (MQ=255)
                                                                 gTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatg                                                                                 <  1:1084625/149‑1 (MQ=255)
                                                                  tCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatga                                                                                <  2:967079/149‑1 (MQ=255)
                                                                   cACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatgat                                                                               <  3:643847/149‑1 (MQ=255)
                                                                   cACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatgat                                                                               <  2:812314/149‑1 (MQ=255)
                                                                    aCCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATc                                                                              >  2:671157/1‑149 (MQ=255)
                                                                      cTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAg                                                                            >  2:706365/1‑149 (MQ=255)
                                                                      cTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAg                                                                            >  3:747194/1‑149 (MQ=255)
                                                                      cTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAg                                                                            >  2:788296/1‑149 (MQ=255)
                                                                       tGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCTGGGa                                                                                                                             >  3:932818/1‑99 (MQ=255)
                                                                       tGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGa                                                                                                                             <  4:932818/99‑1 (MQ=255)
                                                                       tGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGt                                                                           <  2:223469/149‑1 (MQ=255)
                                                                       tGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGt                                                                           <  2:1044935/149‑1 (MQ=255)
                                                                           gcgagcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCa                                                                       <  3:511414/149‑1 (MQ=255)
                                                                            cgagcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCAc                                                                      <  1:319106/149‑1 (MQ=255)
                                                                                 gaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGc                                                                 >  2:378304/1‑149 (MQ=255)
                                                                                    ggAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCaaa                                                              >  1:87088/1‑149 (MQ=255)
                                                                                      aaGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAAt                                                            >  3:440985/1‑149 (MQ=255)
                                                                                       aGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATg                                                           >  3:643257/1‑149 (MQ=255)
                                                                                        gATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGa                                                          <  1:706365/149‑1 (MQ=255)
                                                                                         aTCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                                                              >  3:1017003/1‑80 (MQ=255)
                                                                                         aTCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                                                              <  4:1017003/80‑1 (MQ=255)
                                                                                         aTCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATgg                                                                                                                          >  1:414219/1‑84 (MQ=255)
                                                                                         aTCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATgg                                                                                                                          <  2:414219/84‑1 (MQ=255)
                                                                                         aTCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGaa                                                         <  3:834309/149‑1 (MQ=255)
                                                                                         aTCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGaa                                                         <  1:639553/149‑1 (MQ=255)
                                                                                          tCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcggc                                                                                                                      >  4:128782/1‑87 (MQ=255)
                                                                                          tCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcggc                                                                                                                      <  3:128782/87‑1 (MQ=255)
                                                                                            agcagcGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCa                                                                                        <  4:893445/115‑1 (MQ=255)
                                                                                            agcagcGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCa                                                                                        >  3:893445/1‑115 (MQ=255)
                                                                                               agcgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGaa                                                   <  1:981084/149‑1 (MQ=255)
                                                                                                gcgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGGTTCACCCAGCAAAATGAACAGGAAg                                                  >  3:95667/1‑149 (MQ=255)
                                                                                                      cgcgcCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGAcc                                            <  4:163512/149‑1 (MQ=255)
                                                                                                         gccgcTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACGCTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCa                                         <  1:183346/149‑1 (MQ=255)
                                                                                                           cgcTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAtt                                       <  1:162950/149‑1 (MQ=255)
                                                                                                                  cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCgg                                               >  3:370873/1‑134 (MQ=255)
                                                                                                                  cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCgg                                               <  4:370873/134‑1 (MQ=255)
                                                                                                                  cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                                <  4:95667/149‑1 (MQ=255)
                                                                                                                  cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                                >  3:580511/1‑149 (MQ=255)
                                                                                                                  cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAATTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                                >  2:533281/1‑149 (MQ=255)
                                                                                                                   tGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCa                               <  3:164084/149‑1 (MQ=255)
                                                                                                                    gggTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAg                              <  2:164795/149‑1 (MQ=255)
                                                                                                                    gggTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAg                              <  4:649534/149‑1 (MQ=255)
                                                                                                                      gTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGt                            <  3:640101/149‑1 (MQ=255)
                                                                                                                      gTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGt                            <  1:206328/149‑1 (MQ=255)
                                                                                                                          gcgcGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGtgctg                        <  2:22870/149‑1 (MQ=255)
                                                                                                                              gTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAAc                    <  3:911915/149‑1 (MQ=255)
                                                                                                                               ttAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACg                   >  1:752254/1‑149 (MQ=255)
                                                                                                                               ttAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACg                   >  1:553018/1‑149 (MQ=255)
                                                                                                                                   gTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCg               <  2:786515/149‑1 (MQ=255)
                                                                                                                                    ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg              >  3:969526/1‑149 (MQ=255)
                                                                                                                                    ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg              <  2:608197/149‑1 (MQ=255)
                                                                                                                                     tCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCggg             <  3:723059/149‑1 (MQ=255)
                                                                                                                                       tATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGCTGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCCGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAAAGCTGGGGAt           >  2:343077/1‑149 (MQ=255)
                                                                                                                                         ttttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGt                                                                                                            <  4:198025/50‑1 (MQ=255)
                                                                                                                                         ttttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGt                                                                                                            >  3:198025/1‑50 (MQ=255)
                                                                                                                                         ttttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTc         >  3:611108/1‑149 (MQ=255)
                                                                                                                                         ttttACTGGCGATGGCCAGACCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTc         >  1:859497/1‑149 (MQ=255)
                                                                                                                                           ttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGATGCTGCAACGCTCGGGATTCGa       <  1:671157/149‑1 (MQ=255)
                                                                                                                                           ggaCTGGCCATCGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGa       >  4:578923/3‑149 (MQ=255)
                                                                                                                                           ggaCTGGCCATCGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGa       >  4:912846/3‑149 (MQ=255)
                                                                                                                                               tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCTCCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCg   >  1:356112/1‑149 (MQ=255)
                                                                                                                                               tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCg   >  3:457801/1‑149 (MQ=255)
                                                                                                                                               tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCg   >  3:489853/1‑149 (MQ=255)
                                                                                                                                               tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCg   >  3:573396/1‑149 (MQ=255)
                                                                                                                                               tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGGTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCg   >  4:455152/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGC  >  NC_000913/2616459‑2616751

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