Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,485,560 T→G E131A (GAG→GCG)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,485,5600TG97.3% 637.8 / 2.1 184E131A (GAG→GCG) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (80/99);  total (82/102)
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.81e-01

Only 100 of 186 total aligned reads displayed.

GGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCG  >  NC_000913/4485463‑4485656
                                                                                                 |                                                                                                
ggTCGGCACGTTGAACACCATCGTACGCAGCGGACGATGCGGTTCGCTTTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTcgcg                                                                                                <  2:3720428/100‑1 (MQ=255)
  tCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGAc                                                                                              <  1:3687710/100‑1 (MQ=255)
  tCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGAc                                                                                              <  1:3441163/100‑1 (MQ=255)
    ggCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAg                                                                                            >  1:2714284/1‑100 (MQ=255)
    ggCACGCTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAg                                                                                            <  2:4410552/100‑1 (MQ=255)
      cACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGcc                                                                                          <  1:4012190/100‑1 (MQ=255)
      cACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGcc                                                                                          <  1:3113297/100‑1 (MQ=255)
       aCGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGccn                                                                                         <  1:110919/100‑2 (MQ=255)
        cGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTg                                                                                        >  1:467069/1‑100 (MQ=255)
         gTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGa                                                                                       >  1:3733761/1‑100 (MQ=255)
         gTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGa                                                                                       <  1:1347390/100‑1 (MQ=255)
          ttGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGAc                                                                                      <  2:2673331/100‑1 (MQ=255)
          ttGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGAc                                                                                      >  2:2482289/1‑100 (MQ=255)
            gAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGc                                                                                    >  1:1943812/1‑100 (MQ=255)
              acacCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCAc                                                                                  >  1:3380364/1‑100 (MQ=255)
              acacCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCAc                                                                                  >  1:4502557/1‑100 (MQ=255)
               cacCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACt                                                                                 <  1:2110581/100‑1 (MQ=255)
               cacCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACt                                                                                 <  1:4368054/100‑1 (MQ=255)
                acCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACtt                                                                                <  2:1271029/100‑1 (MQ=255)
                 ccATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACttt                                                                               >  1:3881291/1‑100 (MQ=255)
                 ccATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACttt                                                                               >  1:3509519/1‑100 (MQ=255)
                 ccATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACttt                                                                               >  1:833098/1‑100 (MQ=255)
                   aTCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTcc                                                                             <  1:5431550/100‑1 (MQ=255)
                   aTCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTcc                                                                             >  1:3987248/1‑100 (MQ=255)
                   aTCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTcc                                                                             <  1:4870542/100‑1 (MQ=255)
                     cTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCag                                                                           <  2:3476987/100‑1 (MQ=255)
                       tACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGCCAGCCTGACGCACTTTCCagta                                                                         >  1:1881760/1‑100 (MQ=255)
                        aCGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCagtag                                                                        <  1:1610017/100‑1 (MQ=255)
                        aCGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCagtag                                                                        <  2:2574414/100‑1 (MQ=255)
                        aCGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCagtag                                                                        <  2:700734/100‑1 (MQ=255)
                         cGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCagtagt                                                                       >  2:4924372/1‑100 (MQ=255)
                         cGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCagtagt                                                                       >  2:2371488/1‑100 (MQ=255)
                          gCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAgtt                                                                      >  2:1843366/1‑100 (MQ=255)
                           cAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAgttg                                                                     >  2:4104517/1‑100 (MQ=255)
                           cAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAgttg                                                                     <  2:873264/100‑1 (MQ=255)
                           cAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAgttg                                                                     <  2:4135173/100‑1 (MQ=255)
                           cAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAgttg                                                                     <  1:3347084/100‑1 (MQ=255)
                            aGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAgttgt                                                                    <  1:1928607/100‑1 (MQ=255)
                             gCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAgttgtt                                                                   <  1:2357568/100‑1 (MQ=255)
                               ggACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTAc                                                                 <  2:2065278/100‑1 (MQ=255)
                               ggACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTAc                                                                 <  1:5455102/100‑1 (MQ=255)
                                gacgacGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACg                                                                >  1:4285544/1‑100 (MQ=255)
                                  cgacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACgag                                                              >  2:3068757/1‑98 (MQ=255)
                                  cgacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGc                                                              >  2:318788/1‑100 (MQ=255)
                                  cgacgCGGTTCGCTCTTGGTCGTTTTCAGCTGATCGAAACTGTAGGGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTCGTTGTTACGGc                                                              <  2:2494456/100‑1 (MQ=255)
                                   gacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGcc                                                             >  1:3905746/1‑100 (MQ=255)
                                    acgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCt                                                            >  2:4485656/1‑100 (MQ=255)
                                     cgcgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCtt                                                           <  1:2882943/100‑1 (MQ=255)
                                     cgcgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCtt                                                           >  2:2584985/1‑100 (MQ=255)
                                     cgcgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCtt                                                           >  2:1241959/1‑100 (MQ=255)
                                     cgcgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCtt                                                           <  1:12565/100‑1 (MQ=255)
                                      gcgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCttt                                                          <  1:4745573/100‑1 (MQ=255)
                                         gTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAAc                                                       >  1:146610/1‑100 (MQ=255)
                                         gTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAAc                                                       <  2:577489/100‑1 (MQ=255)
                                          ttCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACg                                                      <  1:5310550/100‑1 (MQ=255)
                                          ttCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACg                                                      <  1:5545526/100‑1 (MQ=255)
                                           tCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGt                                                     <  1:1787053/100‑1 (MQ=255)
                                           tCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGt                                                     <  1:2914317/100‑1 (MQ=255)
                                             gCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGc                                                   <  2:4537503/100‑1 (MQ=255)
                                             gCTCTTGGTCGGTTTCAGCTGATCGAAACTGTAGAGCGACTCTTTTGCCGTCGCGACAGACTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACgcgc                                                   >  2:1907804/1‑100 (MQ=255)
                                                ctTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGc                                                >  2:1342505/1‑100 (MQ=255)
                                                    ttCGCTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCtcag                                            <  1:3085118/100‑1 (MQ=255)
                                                      cGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCtcagtc                                          <  2:2187595/100‑1 (MQ=255)
                                                        ttttCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCtcagtcag                                        >  1:2450081/1‑100 (MQ=255)
                                                           tCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGaaa                                     <  2:1968757/100‑1 (MQ=255)
                                                             aGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGc                                   <  2:4242905/100‑1 (MQ=255)
                                                             aGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGc                                   >  1:1737608/1‑100 (MQ=255)
                                                             aGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGc                                   >  1:3820676/1‑100 (MQ=255)
                                                              gCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCa                                  <  2:5276187/100‑1 (MQ=255)
                                                              gCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCa                                  >  1:1072617/1‑100 (MQ=255)
                                                               cTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAg                                 <  1:4155027/100‑1 (MQ=255)
                                                               cTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAg                                 <  2:5191489/100‑1 (MQ=255)
                                                                 gATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGAc                               <  1:260389/100‑1 (MQ=255)
                                                                 gATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGAc                               >  2:1707299/1‑100 (MQ=255)
                                                                 gATCGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGAc                               <  1:4541999/100‑1 (MQ=255)
                                                                    cGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGc                            >  2:2455693/1‑100 (MQ=255)
                                                                    cGAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGc                            >  2:3744237/1‑100 (MQ=255)
                                                                     gAAACTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCGCTTTCCGGTAGTTGTTACGGCCTTTAACGTGCAGCTCCGTCAGAAAGCAGACCGCt                           >  1:5395163/1‑100 (MQ=255)
                                                                       aaCTGTAGAGCGTCTCTTTTGACGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTACCGTGCAGCTCAGTCAGAAAGCAGACCGCTTc                         <  2:1097950/100‑1 (MQ=255)
                                                                        aCTGTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTcc                        >  1:2215064/1‑100 (MQ=255)
                                                                           gTAGAGCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCAtt                     <  2:1620437/100‑1 (MQ=255)
                                                                             agagCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGa                   <  2:2430346/100‑1 (MQ=255)
                                                                              gagCGTCTCTTTTGCCGTCGCGACAGCTTGACGCAGTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTGTGGAATGCGACCGCTTCCATTGAg                  >  2:3769827/1‑100 (MQ=255)
                                                                               agCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTTTTACGGCCATTAACGTGCAGCTCAGTTAGAAAGTCGATCGCTTCCTTTGAGc                 >  2:1405841/1‑100 (MQ=255)
                                                                               agCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGc                 <  1:3513388/100‑1 (MQ=255)
                                                                               agCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGc                 <  1:3224647/100‑1 (MQ=255)
                                                                                gCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGcc                <  1:32045/100‑1 (MQ=255)
                                                                                gCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGcc                >  2:4608715/1‑100 (MQ=255)
                                                                                gCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGcc                <  1:2988591/100‑1 (MQ=255)
                                                                                gCGTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGcc                <  1:2417494/100‑1 (MQ=255)
                                                                                 cGTCTCTTTTGCCGTCGCGACAGCCTGACGCATTTTCCAGTAGTCGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCa               <  1:5107834/100‑1 (MQ=255)
                                                                                  gTCTCTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAg              >  1:1787702/1‑100 (MQ=255)
                                                                                   tctcTTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGt             <  2:3630154/100‑1 (MQ=255)
                                                                                    ctctTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCATTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTa            >  2:4981431/1‑100 (MQ=255)
                                                                                     tctTTTGCCGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTTCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTAt           >  1:3362401/1‑100 (MQ=255)
                                                                                      ctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATc          >  2:3380215/1‑100 (MQ=255)
                                                                                        tttGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCAt        >  1:5234480/1‑100 (MQ=255)
                                                                                            ccGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAg    <  1:599294/100‑1 (MQ=255)
                                                                                            ccGTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAg    >  2:463949/1‑100 (MQ=255)
                                                                                              gTCGCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCg  <  2:4918671/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCG  >  NC_000913/4485463‑4485656

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