Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,469 +G intergenic (+36/‑197) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4661.G98.3% 755.5 / 2.9 237intergenic (+33/‑200)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (109/123);  total (112/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGAC  >  NC_000913/4008399‑4008551
                                                                        |                                                                                          
actgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg                                                                     >  1:7427623/5‑93 (MQ=255)
actgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg                                                                     >  2:6918389/5‑93 (MQ=255)
 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGCAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATCGGc                                                                    >  2:7442846/4‑93 (MQ=255)
 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc                                                                    <  1:1016132/90‑1 (MQ=255)
 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc                                                                    <  1:7346012/90‑1 (MQ=255)
 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc                                                                    <  1:3790479/90‑1 (MQ=255)
  tgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa                                                                   <  2:5431870/91‑1 (MQ=255)
   gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAATGCTGCATTATGGGCAg                                                                  <  2:9560886/92‑1 (MQ=255)
   gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg                                                                  <  2:3520406/92‑1 (MQ=255)
    ggCTACAGCCTGAAGGATTGTATTGAGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt                                                                 <  2:5907352/93‑1 (MQ=255)
     gCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg                                                                <  1:1490113/93‑1 (MQ=255)
     gCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg                                                                <  2:9885562/93‑1 (MQ=255)
     gCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg                                                                <  1:5093849/93‑1 (MQ=255)
     gCTACAGCCTGAAGGATTGTATTGCGTTTGGTGAAGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTg                                                                <  1:6885365/93‑1 (MQ=255)
      cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc                                                               >  2:3591843/1‑92 (MQ=255)
      cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc                                                               <  1:6556074/93‑2 (MQ=255)
       tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg                                                              >  1:2286645/1‑91 (MQ=255)
       tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg                                                              <  1:6646566/93‑3 (MQ=255)
           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacc                                                          <  2:9121825/91‑1 (MQ=255)
            ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacca                                                         >  2:2557917/2‑93 (MQ=255)
             cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg                                                        <  1:8986337/93‑1 (MQ=255)
              tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                       >  1:4424522/1‑93 (MQ=255)
              tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                       >  1:4394163/1‑93 (MQ=255)
               gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGTTGCATTATGGGCAGTGCGCACCAGCg                                                      >  1:2268480/1‑93 (MQ=255)
               gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                      >  2:6486777/1‑93 (MQ=255)
               gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                      >  1:338086/1‑93 (MQ=255)
                aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                                     >  1:8016447/1‑93 (MQ=255)
                aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                                     <  2:1114875/93‑1 (MQ=255)
                 aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                    >  1:760100/1‑93 (MQ=255)
                 aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                    <  1:2050830/93‑1 (MQ=255)
                 aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                    <  1:2250235/93‑1 (MQ=255)
                 aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                    >  2:7801316/1‑93 (MQ=255)
                  ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATTGCGGGGCAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                   >  2:4046463/1‑93 (MQ=255)
                  ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                   >  1:7201063/1‑93 (MQ=255)
                               tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                      >  1:7667105/5‑93 (MQ=255)
                               tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                      <  1:9899297/89‑1 (MQ=255)
                                ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                                     <  2:2733140/90‑1 (MQ=255)
                                ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                                     <  1:4602986/90‑1 (MQ=255)
                                ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                                     <  1:4274339/90‑1 (MQ=255)
                                ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                                     <  2:7428799/90‑1 (MQ=255)
                                 tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                    <  1:5598663/91‑1 (MQ=255)
                                   gtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                  <  1:6171928/93‑1 (MQ=255)
                                    tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                 <  2:212095/93‑1 (MQ=255)
                                    tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                 <  2:2159170/93‑1 (MQ=255)
                                    tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                 >  2:9042858/1‑93 (MQ=255)
                                     gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                <  1:5840377/93‑1 (MQ=255)
                                      aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                               <  1:110867/93‑1 (MQ=255)
                                       cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                              >  1:9448758/1‑93 (MQ=255)
                                         ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                           >  2:4201697/1‑92 (MQ=255)
                                         ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                           <  2:825032/93‑2 (MQ=255)
                                           aTGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                         <  1:2792070/93‑4 (MQ=255)
                                            tGAGCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                        <  1:4050853/93‑5 (MQ=255)
                                            tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCCTCCCGAGCTGGAgg                        >  1:3267306/1‑92 (MQ=255)
                                            tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                        >  2:8102031/1‑89 (MQ=255)
                                            tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                        <  2:1293452/93‑5 (MQ=255)
                                            tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                        >  1:52866/1‑89 (MQ=255)
                                            tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                        <  2:7292451/93‑5 (MQ=255)
                                            tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                        <  2:2070297/93‑5 (MQ=255)
                                            tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                        >  1:6916590/1‑89 (MQ=255)
                                              aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                       >  1:5440669/1‑93 (MQ=255)
                                              aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                       <  2:3069874/93‑1 (MQ=255)
                                               acgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                      <  2:5182193/93‑1 (MQ=255)
                                               acgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                      <  1:6076560/93‑1 (MQ=255)
                                                cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                     >  1:3068276/1‑93 (MQ=255)
                                                cgacgCAGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                     <  1:4798201/93‑1 (MQ=255)
                                                 gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGGGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtg                    >  1:3489390/1‑92 (MQ=255)
                                                 gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                    <  1:966285/93‑1 (MQ=255)
                                                 gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGGGAtt                    >  2:5700810/1‑93 (MQ=255)
                                                  acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                   >  1:4709955/1‑93 (MQ=255)
                                                  acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                   <  2:7867718/93‑1 (MQ=255)
                                                  acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                   <  1:1935044/93‑1 (MQ=255)
                                                   cgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                  <  1:6428980/93‑1 (MQ=255)
                                                    gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                 <  1:8571730/93‑1 (MQ=255)
                                                    gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                 <  2:7233090/93‑1 (MQ=255)
                                                    gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                 >  2:1897124/1‑93 (MQ=255)
                                                     ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                <  2:4507461/93‑1 (MQ=255)
                                                     ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                >  2:8626148/1‑93 (MQ=255)
                                                     ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                <  2:2655949/93‑1 (MQ=255)
                                                     ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                <  2:2365581/93‑1 (MQ=255)
                                                     ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                >  1:7642186/1‑93 (MQ=255)
                                                      cGAGATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc               <  2:61351/93‑1 (MQ=255)
                                                      cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc               >  1:9836448/1‑93 (MQ=255)
                                                      cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc               <  2:6240468/93‑1 (MQ=255)
                                                       gAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt              >  1:9007313/1‑93 (MQ=255)
                                                           tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg          >  1:9096104/1‑93 (MQ=255)
                                                           tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg          >  2:1447203/1‑93 (MQ=255)
                                                           tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg          >  1:4621631/1‑93 (MQ=255)
                                                           tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg          <  1:4827014/93‑1 (MQ=255)
                                                            gCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc         <  2:5554618/93‑1 (MQ=255)
                                                             cTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc        >  2:1753246/1‑93 (MQ=255)
                                                             cTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc        >  1:2428242/1‑93 (MQ=255)
                                                              tGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg       >  2:4562025/1‑93 (MQ=255)
                                                               ggcGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga      <  2:8691627/91‑1 (MQ=255)
                                                               gTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga      >  2:4889679/1‑93 (MQ=255)
                                                                tCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac     >  1:3519909/1‑93 (MQ=255)
                                                                tCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac     >  1:9472067/1‑93 (MQ=255)
                                                                  gATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga   <  2:8553131/93‑1 (MQ=255)
                                                                  gATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga   >  2:8910148/1‑93 (MQ=255)
                                                                  gATGGCGGAG‑‑AAAGGCTTCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga   <  2:1592122/93‑1 (MQ=255)
                                                                   aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac  >  2:6387962/1‑93 (MQ=255)
                                                                        |                                                                                          
AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGAC  >  NC_000913/4008399‑4008551

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