Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,239 A→T 16.2% intergenic (+332/+40) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,2390AT16.2% 372.5 / 56.0 149intergenic (+332/+40)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (69/55);  new base T (15/9);  total (84/65)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 155 total aligned reads displayed.

TGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCG  >  NC_000913/1224140‑1224336
                                                                                                   |                                                                                                 
tgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAgaga                                                                                                   >  2:1013121/1‑100 (MQ=255)
 gtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                  >  2:748705/1‑100 (MQ=255)
 gtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAAAGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                  >  1:3290654/1‑100 (MQ=255)
  tAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATa                                                                                                 <  1:718474/100‑1 (MQ=255)
      gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                                             <  2:3137226/100‑1 (MQ=255)
      gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                                             >  1:2169859/1‑100 (MQ=255)
        gggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTg                                                                                           >  1:3729830/1‑100 (MQ=255)
          gAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                         >  2:3016418/1‑100 (MQ=255)
            ggTGGAAATAGCCATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCCAGGCAGAGATAACTCTGCCtt                                                                                       >  1:164702/1‑100 (MQ=255)
            ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAGTCTGCCtt                                                                                       >  2:1520000/1‑100 (MQ=255)
            ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                       <  2:857852/100‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAGTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                      <  1:1281009/100‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                      <  1:3584754/100‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                       >  1:3706684/1‑99 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                      <  2:3779506/100‑1 (MQ=255)
              tGGAAATAGCAATGATGACTATCAGCAAGAATACTCGCCTCTTTAGCATAACGTGGATGAGAGGGATGAAACACTCAAGGCAGAGATAACTCTGCCTTg                                                                                      <  1:2271451/99‑1 (MQ=255)
              tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCTCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGa                                                                                     <  2:2681753/100‑1 (MQ=39)
              tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGa                                                                                     >  1:3240207/1‑100 (MQ=255)
               ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGag                                                                                    <  2:2823330/100‑2 (MQ=255)
               ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGag                                                                                    <  1:124336/100‑2 (MQ=255)
               ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                    >  1:2812576/1‑100 (MQ=255)
                gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                   <  2:2859459/100‑1 (MQ=255)
                   ttccgatctaggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                <  2:1533155/91‑1 (MQ=255)
                   aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagtttt                                                                                <  1:3803401/100‑6 (MQ=38)
                   aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                <  2:2347896/100‑1 (MQ=255)
                    tAGCAATGAGGGGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa                                                                               <  2:2792995/100‑1 (MQ=255)
                    tAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa                                                                               >  2:2548101/1‑100 (MQ=255)
                      gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                             >  1:3275638/1‑100 (MQ=255)
                      gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                             >  2:505598/1‑100 (MQ=255)
                      gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                             <  2:1681698/100‑1 (MQ=255)
                            aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaagatcggaa                                                                       >  1:1533146/1‑92 (MQ=255)
                            aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                                                       >  2:2411453/1‑100 (MQ=255)
                             ggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                      <  2:31206/100‑1 (MQ=255)
                             ggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                      <  2:1426194/100‑1 (MQ=255)
                               agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                    <  2:2473579/100‑1 (MQ=255)
                               acgaTCAACAAGAATACTCGCCGCTTTAACACAACGTTGATGAGAAGGATGAAAAACTCAAGGCAGAGATAACTCTTCCTTTAAGATAAATgcg                                                                          <  1:3163888/91‑1 (MQ=255)
                                gTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCCAGGCAGAGATCACTCTGCCTTGAAGATAAATGCGCTTTTAc                                                                   >  1:2300992/1‑100 (MQ=255)
                                  aTCAGCAAGACTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAGTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                 <  1:2502416/100‑1 (MQ=255)
                                  aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAc                                                                   <  2:867290/98‑1 (MQ=255)
                                  aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                 >  1:867283/1‑100 (MQ=255)
                                  aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                 >  1:594181/1‑100 (MQ=255)
                                  aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                 >  1:2963626/1‑100 (MQ=255)
                                   tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                                                >  2:1536244/1‑100 (MQ=255)
                                     aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                              >  2:2010939/1‑100 (MQ=255)
                                      gCAAGAGTACTGGACGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                             <  1:3046462/100‑1 (MQ=255)
                                       cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                            <  2:3780658/100‑1 (MQ=255)
                                       cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                            <  1:2932474/100‑1 (MQ=255)
                                       cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                            <  2:3721633/100‑1 (MQ=255)
                                       cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGAATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                              <  1:436081/98‑1 (MQ=255)
                                        aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                           <  2:1459618/100‑1 (MQ=255)
                                          gAATACTCTCCGCTTTACCACACCGTGAATGAGAGGCATCAAAAACCCAAGGCAGAGATAACTCGGCGTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                           <  1:2509576/98‑1 (MQ=255)
                                          gAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTa                                                         >  2:1764215/1‑100 (MQ=255)
                                                tcgtggcTTAACCACAACATGGATGAGGGTGATGGAAGATTGAAGGCAGAGATAACTCTTCATTGAAGATAAATGCGCTTTTACAGCGGGGTTATTTc                                                     <  1:301082/93‑1 (MQ=255)
                                                 cgcctcTTCACCACAACTTGGATGAGAGTGATGAAAGACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                   <  1:2562736/99‑1 (MQ=255)
                                                  gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGct                                                 >  2:1654661/1‑100 (MQ=255)
                                                    cgcTTTACCACAACGTGGATGAGAGGGATGAAAATCTCAAAGCAAAGATAACTCTGCCTTGAAGATAAATGCGTTTTTACAGCTGGCTTATTTCAGc                                                  <  1:2096295/97‑1 (MQ=255)
                                                      ntttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGCCAGAGATAACACTGCATTGTAGATCAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                             >  1:1473337/2‑100 (MQ=255)
                                                       tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttcn                                            <  2:3293780/100‑2 (MQ=255)
                                                        ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                           <  2:3669290/100‑1 (MQ=255)
                                                           ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATGACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                        >  1:557602/1‑100 (MQ=255)
                                                           ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                        >  2:1164994/1‑100 (MQ=255)
                                                           ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                        >  1:3296865/1‑100 (MQ=255)
                                                           ccACAACATGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                        >  2:1997028/1‑100 (MQ=255)
                                                             cttccgatcttGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                      <  1:351443/90‑1 (MQ=255)
                                                             acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                      >  1:786775/1‑100 (MQ=255)
                                                             acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                      >  2:1793013/1‑100 (MQ=255)
                                                              caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg                                     >  1:2237719/1‑100 (MQ=255)
                                                                aCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                   >  1:993519/1‑100 (MQ=255)
                                                                aCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                   <  1:630733/100‑1 (MQ=255)
                                                                   tGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                >  1:1141264/1‑100 (MQ=255)
                                                                   tGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                >  2:1314963/1‑100 (MQ=255)
                                                                      aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTCCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc                             >  1:2043885/1‑100 (MQ=255)
                                                                      aTGAGAGGGATGAAAAACTCAAAGCAGAGATAACTCTGCCCTGAAGATAAATGCGCTTTTACAGCGGACTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc                             >  1:1694485/1‑100 (MQ=255)
                                                                       tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTccagatcggaag                            >  2:351448/1‑90 (MQ=255)
                                                                       tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCa                            <  2:1141445/100‑1 (MQ=255)
                                                                          gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGt                         >  2:2650591/1‑100 (MQ=255)
                                                                            gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                       <  2:3697808/100‑1 (MQ=255)
                                                                            gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTCTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                       >  1:227994/1‑100 (MQ=255)
                                                                              gATGAAAACCTCGAGGCAGAGGTATCTCTGTCTTTATGATAATTGCGCTGTTACAGCGGGCTTATTTCAGGTCTTTTGCTTCCGGTAAGGTCActt                         <  1:3089443/96‑3 (MQ=255)
                                                                               catttatcttCAAGTCAGAGATATCTCTGTTTTGAAGATAACTGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGtt                        <  1:803364/87‑1 (MQ=255)
                                                                                  aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                 >  2:3359485/1‑100 (MQ=255)
                                                                                   aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCaa                >  2:980358/1‑100 (MQ=255)
                                                                                    aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCaa                <  1:1064652/99‑1 (MQ=255)
                                                                                      aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa             >  2:2324660/1‑100 (MQ=255)
                                                                                        tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa           >  1:1353820/1‑100 (MQ=255)
                                                                                        tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa           >  2:840807/1‑100 (MQ=255)
                                                                                        tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa           >  1:3135998/1‑100 (MQ=255)
                                                                                         cAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCACGAATAg          >  1:120635/1‑100 (MQ=255)
                                                                                          aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa         >  2:3565802/1‑100 (MQ=255)
                                                                                          aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTCTTTCAGCTCTTCTGCTTCCGGTACGGTCACGTTCAGCTCAAGAATAGa         >  1:2296800/1‑100 (MQ=255)
                                                                                          aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa         <  1:2658502/100‑1 (MQ=255)
                                                                                          aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa         >  1:1946199/1‑100 (MQ=255)
                                                                                           aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa        >  1:3552809/1‑100 (MQ=255)
                                                                                           aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa        >  2:1341057/1‑100 (MQ=255)
                                                                                            ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa       >  2:3751770/1‑100 (MQ=255)
                                                                                              cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAata     >  1:3047926/1‑100 (MQ=255)
                                                                                                gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATc   >  1:2508665/1‑100 (MQ=255)
                                                                                                gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATc   >  2:564543/1‑100 (MQ=255)
                                                                                                 agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg  >  2:2761141/1‑100 (MQ=255)
                                                                                                 agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg  <  2:3263505/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
TGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCG  >  NC_000913/1224140‑1224336

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