Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,242 A→T 17.7% intergenic (+335/+37) 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,2420AT17.7% 368.9 / 62.7 150intergenic (+335/+37)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (70/53);  new base T (15/12);  total (85/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 152 total aligned reads displayed.

AAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATC  >  NC_000913/1224143‑1224341
                                                                                                   |                                                                                                   
aaaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                                                     >  2:616601/1‑100 (MQ=255)
   gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                                                  <  2:3137226/100‑1 (MQ=255)
     gggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTg                                                                                                >  1:3729830/1‑100 (MQ=255)
       gAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                              >  2:3016418/1‑100 (MQ=255)
        aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                             <  2:3303061/100‑1 (MQ=255)
         ggTGGAAATAGCCATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCCAGGCAGAGATAACTCTGCCtt                                                                                            >  1:164702/1‑100 (MQ=255)
         ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAGTCTGCCtt                                                                                            >  2:1520000/1‑100 (MQ=255)
         ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                            <  2:857852/100‑1 (MQ=255)
         ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                            >  1:492271/1‑100 (MQ=255)
          gTGGAAATAGCAATGAGGAGTATCAGCAAGAGTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                           <  1:1281009/100‑1 (MQ=255)
          gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                           <  1:3584754/100‑1 (MQ=255)
          gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                           <  2:3779506/100‑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)
             gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                        <  2:2859459/100‑1 (MQ=255)
              aacaaGCAATGAGGAGTAACAGCAAGAATACAAGCCGCTTTACCACAACGTGTATAAAAAAAAAAAAAAACTCAAGGCAGAGTTATCTCTGCCTTGa                                                                                          <  1:3041759/93‑1 (MQ=16)
               aaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                      >  1:31204/1‑100 (MQ=255)
               aaTAGCAATGAGGAGTATCAGCAAGAATACTCGACGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                      >  1:2288680/1‑100 (MQ=255)
                aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                     <  2:2347896/100‑1 (MQ=255)
                 tAGCAATGAGGGGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa                                                                                    <  2:2792995/100‑1 (MQ=255)
                 tAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa                                                                                    >  2:2548101/1‑100 (MQ=255)
                   gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                  <  2:1681698/100‑1 (MQ=255)
                   gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                  >  2:505598/1‑100 (MQ=255)
                   gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                  >  1:2710164/1‑100 (MQ=255)
                   gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGACAGAGATAACTCTGCCTTGAAGATaaa                                                                                   >  1:991125/1‑99 (MQ=255)
                    cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTAGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                                                 <  1:859720/100‑1 (MQ=255)
                         aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaagatcggaa                                                                            >  1:1533146/1‑92 (MQ=255)
                         aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                                                            >  2:2411453/1‑100 (MQ=255)
                          ggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                           <  2:1426194/100‑1 (MQ=255)
                           gagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtttt                                                                          <  1:610385/100‑1 (MQ=255)
                            agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                         >  1:2947138/1‑100 (MQ=255)
                             gTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCCAGGCAGAGATCACTCTGCCTTGAAGATAAATGCGCTTTTAc                                                                        >  1:2300992/1‑100 (MQ=255)
                               aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAc                                                                        <  2:867290/98‑1 (MQ=255)
                               aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                      >  1:594181/1‑100 (MQ=255)
                                tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                                                     >  2:1536244/1‑100 (MQ=255)
                                  aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                   >  2:2010939/1‑100 (MQ=255)
                                    cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                                 <  2:3721633/100‑1 (MQ=255)
                                    cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                                 <  1:2932474/100‑1 (MQ=255)
                                    cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGAATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                   <  1:436081/98‑1 (MQ=255)
                                     aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                                <  2:1459618/100‑1 (MQ=255)
                                           aCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                          >  1:322731/1‑100 (MQ=255)
                                             tcgtggcTTAACCACAACATGGATGAGGGTGATGGAAGATTGAAGGCAGAGATAACTCTTCATTGAAGATAAATGCGCTTTTACAGCGGGGTTATTTc                                                          <  1:301082/93‑1 (MQ=255)
                                              cgcctcTTCACCACAACTTGGATGAGAGTGATGAAAGACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                        <  1:2562736/99‑1 (MQ=255)
                                               gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGct                                                      <  2:3549913/100‑1 (MQ=255)
                                                 cgcTTTACCACAACGTGGATGAGAGGGATGAAAATCTCAAAGCAAAGATAACTCTGCCTTGAAGATAAATGCGTTTTTACAGCTGGCTTATTTCAGc                                                       <  1:2096295/97‑1 (MQ=255)
                                                   ntttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGCCAGAGATAACACTGCATTGTAGATCAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                  >  1:1473337/2‑100 (MQ=255)
                                                    tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttcn                                                 <  2:3293780/100‑2 (MQ=255)
                                                     ttACCACAACGTGGATGAGAGGGATGAAAATCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATACGCTTTTACAGCAGGCTTATTTCAGCTGTtgtg                                                >  2:1754707/1‑100 (MQ=255)
                                                     ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                                <  2:3669290/100‑1 (MQ=255)
                                                        ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                             >  1:3296865/1‑100 (MQ=255)
                                                        ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                             >  2:1164994/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                                          <  2:3047712/100‑1 (MQ=255)
                                                           caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg                                          >  1:2237719/1‑100 (MQ=255)
                                                             aCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                        >  1:993519/1‑100 (MQ=255)
                                                               gTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTaa                                      <  1:2146256/100‑1 (MQ=255)
                                                               gTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTaa                                      <  1:2377896/100‑1 (MQ=255)
                                                                tGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                     >  1:1141264/1‑100 (MQ=255)
                                                                tGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                     >  2:1314963/1‑100 (MQ=255)
                                                                  gATGAGAGGTATGAAAAACTCAAGGCAGATATATCTCTGCCTTGAATATAAATGCTCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                     <  1:2506354/98‑1 (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:322736/100‑1 (MQ=255)
                                                                       gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGCCTTATTTCAGCTATTCTGCTTCCGGTAAGGTCACGt                              >  1:3285990/1‑100 (MQ=255)
                                                                         gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                            >  1:3085385/1‑100 (MQ=255)
                                                                         gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTCTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                            >  1:227994/1‑100 (MQ=255)
                                                                                aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCaa                     >  2:980358/1‑100 (MQ=255)
                                                                                 aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCaa                     <  1:1064652/99‑1 (MQ=255)
                                                                                 aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                    <  2:3729854/100‑1 (MQ=255)
                                                                                 aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                    <  2:2411882/100‑1 (MQ=255)
                                                                                  tcttCAAGGCAGAGTTATCTCTGCCTTGAGGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                   >  1:2546748/4‑100 (MQ=38)
                                                                                  aaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                   <  1:2005137/100‑1 (MQ=255)
                                                                                   aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa                  >  2:2324660/1‑100 (MQ=255)
                                                                                     tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa                >  2:264896/1‑100 (MQ=255)
                                                                                     tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa                >  2:840807/1‑100 (MQ=255)
                                                                                     tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa                >  1:1353820/1‑100 (MQ=255)
                                                                                     tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa                >  1:3135998/1‑100 (MQ=255)
                                                                                     tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa                <  2:1576829/100‑1 (MQ=255)
                                                                                      taaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg               <  2:1946209/99‑1 (MQ=255)
                                                                                      cAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCACGAATAg               >  1:120635/1‑100 (MQ=255)
                                                                                       aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa              >  2:3692538/1‑100 (MQ=255)
                                                                                       aaGGCAGAGTTATCTCTGCCCTGAAGATAAATGCGATTTTACCTCGGGCTTATTTCCGCTCTTCTGCTTCCGGTAAGGTCACGTTCCGCTCAACAATAg               >  1:224235/1‑99 (MQ=255)
                                                                                       aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACCGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTACAGCTCAAGAATAGa              >  2:3754740/1‑100 (MQ=255)
                                                                                       aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa              >  1:1946199/1‑100 (MQ=255)
                                                                                       aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa              >  2:395956/1‑100 (MQ=255)
                                                                                        aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa             >  2:1341057/1‑100 (MQ=255)
                                                                                        aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa             >  1:3552809/1‑100 (MQ=255)
                                                                                        aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa             >  2:1529729/1‑100 (MQ=255)
                                                                                        aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa             >  1:3762402/1‑100 (MQ=255)
                                                                                           cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAata          >  1:3047926/1‑100 (MQ=255)
                                                                                           cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAata          >  1:3667539/1‑100 (MQ=255)
                                                                                            agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAatat         >  1:766647/1‑100 (MQ=255)
                                                                                             gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATc        >  1:2508665/1‑100 (MQ=255)
                                                                                              agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg       >  2:2761141/1‑100 (MQ=255)
                                                                                                   aCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATc  <  2:2259392/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATC  >  NC_000913/1224143‑1224341

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