Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,571,921 (G)6→7 intergenic (+6/‑493) yjiS → / → yjiT DUF1127 domain‑containing protein YjiS/putative uncharacterized protein YjiT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,9151.G88.4% 1296.1 / 154.3 344coding (165/165 nt)yjiSDUF1127 domain‑containing protein YjiS
Reads supporting (aligned to +/- strand):  ref base . (10/30);  new base G (103/201);  total (113/231)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 383 total aligned reads displayed.

TGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGAGT  >  NC_000913/4571815‑4572010
                                                                                                     |                                                                                                
   ggCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggg                                                                                               <  2:8457209/101‑1 (MQ=255)
    gCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAgggg                                                                                             <  2:8457754/101‑1 (MQ=255)
    gCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAgggg                                                                                             >  2:8457812/1‑101 (MQ=255)
     cGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggggg                                                                                            <  2:8457280/101‑1 (MQ=255)
       ttGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggggggg                                                                                          <  2:8457342/101‑2 (MQ=255)
       ttGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggggggg                                                                                          <  2:8457427/101‑2 (MQ=255)
        tGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggggggga                                                                                         <  1:8457252/101‑3 (MQ=255)
        tGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggggggga                                                                                         <  1:8457150/101‑3 (MQ=255)
         ggcggcggcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAgggggggat                                                                                        <  2:8457453/101‑4 (MQ=255)
         ggcggcggcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAgggggggat                                                                                        <  2:8457634/101‑4 (MQ=255)
          gcggcggcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAgggggggata                                                                                       <  1:8457021/101‑5 (MQ=255)
           cggcggcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGAtata                                                                                      >  1:8457814/1‑101 (MQ=255)
             gcggcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAg                                                                                     >  2:8457819/1‑101 (MQ=255)
              cggcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGa                                                                                    >  1:8457821/1‑101 (MQ=255)
                gcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGAtt                                                                                  >  2:8457828/1‑101 (MQ=255)
                gcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGAtt                                                                                  <  1:8457144/101‑1 (MQ=255)
                gcAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGAtt                                                                                  <  1:8457296/101‑1 (MQ=255)
                 cAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGAttt                                                                                 <  2:8457238/101‑1 (MQ=255)
                  aaaTGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTta                                                                                >  1:8457831/1‑101 (MQ=255)
                  aaaTGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTta                                                                                <  1:8457076/101‑1 (MQ=255)
                  aaaTGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATCta                                                                                <  1:8456909/101‑1 (MQ=255)
                    aTGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtata                                                                              <  1:8457506/101‑1 (MQ=255)
                     tGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtatat                                                                             <  1:8457197/101‑1 (MQ=255)
                     tGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTtatata                                                                            <  2:8457059/101‑1 (MQ=255)
                      gCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtatata                                                                            >  1:8457836/1‑101 (MQ=255)
                      gCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtatata                                                                            <  1:8457085/101‑1 (MQ=255)
                      gCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtatata                                                                            <  2:8457173/101‑1 (MQ=255)
                       cAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtatatat                                                                           <  2:8456932/101‑1 (MQ=255)
                       cAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtatatat                                                                           <  1:8457366/101‑1 (MQ=255)
                        aGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTtatatata                                                                          <  2:8457306/101‑1 (MQ=255)
                         gACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATataa                                                                         <  1:8456995/101‑1 (MQ=255)
                          aCCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATataat                                                                        >  2:8457840/1‑101 (MQ=255)
                          aCCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATataat                                                                        <  2:8457438/101‑1 (MQ=255)
                            ccGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATAtaataa                                                                      >  1:8457845/1‑101 (MQ=255)
                            ccGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGGGAGGGGGGGATATAGATTTATATAtaataa                                                                      >  2:8457847/1‑101 (MQ=255)
                               aCGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAAc                                                                   <  2:8457565/101‑1 (MQ=255)
                               aCGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAAc                                                                   <  2:8457606/101‑1 (MQ=255)
                                 gcgTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGt                                                                 >  2:8457853/1‑101 (MQ=255)
                                 gcgTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGt                                                                 <  1:8457227/101‑1 (MQ=255)
                                  cgTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGtt                                                                <  1:8457101/101‑1 (MQ=255)
                                   gtgtTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGttt                                                               <  1:8457459/101‑1 (MQ=255)
                                   gtgtTACAGCAGATGAGGGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGttt                                                               >  2:8457858/1‑101 (MQ=255)
                                       tACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTAt                                                           >  2:8457862/1‑101 (MQ=255)
                                        aCAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATg                                                          <  1:8457063/101‑1 (MQ=255)
                                         cagcagATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGt                                                         <  2:8457132/101‑1 (MQ=255)
                                          agcagATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGtt                                                        <  2:8457664/101‑1 (MQ=255)
                                           gcagATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGttt                                                       >  1:8457866/1‑101 (MQ=255)
                                              gATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTa                                                    >  1:8457871/1‑101 (MQ=255)
                                              gATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTa                                                    <  2:8457105/101‑1 (MQ=255)
                                               aTGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTaa                                                   <  2:8457487/101‑1 (MQ=255)
                                               aTGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTaa                                                   <  2:8457321/101‑1 (MQ=255)
                                                tGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTaaa                                                  <  1:8457165/101‑1 (MQ=255)
                                                 gAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAAt                                                 <  1:8457270/101‑1 (MQ=255)
                                                 gAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAAt                                                 <  1:8457566/101‑1 (MQ=255)
                                                 gAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAAt                                                 >  2:8457872/1‑101 (MQ=255)
                                                    tgatgaGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTaa                                              >  2:8457876/1‑101 (MQ=255)
                                                       tgaGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGtt                                           <  1:8457479/101‑1 (MQ=255)
                                                       tgaGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGtt                                           <  1:8457315/101‑1 (MQ=255)
                                                        gaGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTta                                          <  2:8456957/101‑1 (MQ=255)
                                                         aGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTtat                                         <  1:8457324/101‑1 (MQ=255)
                                                           cGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaa                                       <  1:8457493/101‑1 (MQ=255)
                                                            ggTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaaa                                      <  1:8457204/101‑1 (MQ=255)
                                                            ggTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaaa                                      >  2:8457881/1‑101 (MQ=255)
                                                             gTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaaaa                                     <  1:8457540/101‑1 (MQ=255)
                                                             gTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGTGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaaaa                                     <  2:8457261/101‑1 (MQ=255)
                                                              ttGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaaaaa                                    <  1:8457048/101‑1 (MQ=255)
                                                              ttGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaaaaa                                    <  2:8457388/101‑1 (MQ=255)
                                                              ttGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATaaaaa                                    >  2:8457885/1‑101 (MQ=255)
                                                               tGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAAt                                   <  2:8457395/101‑1 (MQ=255)
                                                                gAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAAtt                                  <  2:8457020/101‑1 (MQ=255)
                                                                 aaGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAAttt                                 >  1:8457889/1‑101 (MQ=255)
                                                                 aaGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAAttt                                 <  1:8457290/101‑1 (MQ=255)
                                                                  aGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAAtttt                                <  1:8457175/101‑1 (MQ=255)
                                                                  aGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAAtttt                                >  2:8457893/1‑101 (MQ=255)
                                                                    gATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTcc                              >  2:8457896/1‑101 (MQ=255)
                                                                      tatCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCg                            <  1:8456835/101‑1 (MQ=255)
                                                                      tatCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCg                            <  2:8457220/101‑1 (MQ=255)
                                                                           ggTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACa                       >  2:8457902/1‑101 (MQ=255)
                                                                            gTTACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACaa                      >  1:8457903/1‑101 (MQ=255)
                                                                              tACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAAtt                    <  1:8457733/101‑1 (MQ=255)
                                                                              tACGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAAtt                    <  2:8457190/101‑1 (MQ=255)
                                                                               aCGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATtc                   >  2:8457907/1‑101 (MQ=255)
                                                                                cGCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATtca                  <  1:8457403/101‑1 (MQ=255)
                                                                                 gCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATtcat                 >  2:8457912/1‑101 (MQ=255)
                                                                                 gCAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATtcat                 <  1:8457519/101‑1 (MQ=255)
                                                                                  cAGGGAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATtcatc                >  1:8457915/1‑101 (MQ=255)
                                                                                    gggAGGATGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCaa              <  2:8457550/101‑1 (MQ=255)
                                                                                     ggaggaTGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAAt             >  2:8457921/1‑101 (MQ=255)
                                                                                     ggaggaTGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAAt             >  1:8457922/1‑101 (MQ=255)
                                                                                     ggaggaTGTGGAGTGAGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATa            <  1:8457631/101‑1 (MQ=255)
                                                                                        ggaTGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAgg          <  1:8457425/101‑1 (MQ=255)
                                                                                          aTGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAgg          >  1:8457926/1‑99 (MQ=255)
                                                                                          aTGTGGAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTa        <  1:8457217/101‑1 (MQ=255)
                                                                                          aTGTGGAGTGAGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAg           <  1:8457441/97‑1 (MQ=255)
                                                                                              ggAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGa    <  2:8457374/101‑1 (MQ=255)
                                                                                               gAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGAg   <  2:8457153/101‑1 (MQ=255)
                                                                                               gAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGAg   <  1:8457348/101‑1 (MQ=255)
                                                                                               gAGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGAg   <  1:8457121/101‑1 (MQ=255)
                                                                                                aGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGAgt  >  1:8457932/1‑101 (MQ=255)
                                                                                                aGTGAGGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGAgt  <  1:8457134/101‑1 (MQ=255)
                                                                                                     |                                                                                                
TGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTTATATATAATAAAACGTTTTTATGTTTTTAAATTAAGTTATAAAAATTTTCCCGAGACAATTCATCAATAGGTATGGAGT  >  NC_000913/4571815‑4572010

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