Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,056,373 A→G K129K (AAA→AAG)  yihK → putative GTP‑binding factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,3730AG99.8% 1507.7 / ‑5.0 468K129K (AAA→AAG) yihKputative GTP‑binding factor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (221/245);  total (221/247)
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.99e-01

Only 100 of 478 total aligned reads displayed.

TGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGA  >  NC_000913/4056283‑4056462
                                                                                          |                                                                                         
tGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGt                                                                                         <  2:6847714/93‑1 (MQ=255)
 gACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGtt                                                                                        <  1:5876447/93‑1 (MQ=255)
  aCGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTg                                                                                       >  2:3775795/1‑93 (MQ=255)
  aCGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTg                                                                                       <  2:2202932/93‑1 (MQ=255)
    gCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGAc                                                                                     <  1:6643998/93‑1 (MQ=255)
    gCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGAc                                                                                     <  2:1626937/93‑1 (MQ=255)
     cATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGAcc                                                                                    <  2:9734049/93‑1 (MQ=255)
     cATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGAcc                                                                                    >  2:7631077/1‑93 (MQ=255)
     cATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGAcc                                                                                    <  1:7989381/93‑1 (MQ=255)
      aTTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCg                                                                                   >  2:6433855/1‑93 (MQ=255)
      aTTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCg                                                                                   <  2:3408173/93‑1 (MQ=255)
       tttGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGc                                                                                  <  1:6529328/93‑1 (MQ=255)
       tttGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGc                                                                                  >  1:8728324/1‑93 (MQ=255)
       tttGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGc                                                                                  >  2:1505461/1‑93 (MQ=255)
       tttGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGc                                                                                  >  2:7285444/1‑93 (MQ=255)
          ggcGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCt                                                                               <  1:4232043/91‑1 (MQ=255)
          gACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCt                                                                               >  2:5453297/1‑93 (MQ=255)
           aCGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTg                                                                              >  1:3384274/1‑93 (MQ=255)
           aCGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTg                                                                              <  1:4173151/93‑1 (MQ=255)
           aCGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTg                                                                              <  1:9192642/93‑1 (MQ=255)
             ggCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGgc                                                                            <  1:5655310/93‑1 (MQ=255)
                ccGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGgcgcg                                                                         <  1:2011418/93‑1 (MQ=255)
                 cGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGgcgcgc                                                                        >  2:1975429/1‑93 (MQ=255)
                 cGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGgcgcgc                                                                        <  2:3437027/93‑1 (MQ=255)
                   aTGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGt                                                                      >  2:4562434/1‑93 (MQ=255)
                    tGCCGCGAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTc                                                                     >  2:6805243/1‑93 (MQ=255)
                    tGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTc                                                                     <  2:9330719/93‑1 (MQ=255)
                       cGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTg                                                                  <  2:8348515/93‑1 (MQ=255)
                       cGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTg                                                                  >  2:8741059/1‑93 (MQ=255)
                        gCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTgt                                                                 >  1:7582667/1‑92 (MQ=255)
                        gCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGa                                                                 >  1:704740/1‑93 (MQ=255)
                         cAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGAt                                                                >  1:4121934/1‑93 (MQ=255)
                           aaCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTg                                                              >  2:8273734/1‑93 (MQ=255)
                            aCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTgg                                                             <  2:7135912/93‑1 (MQ=255)
                              gcgcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGt                                                           <  2:8007332/93‑1 (MQ=255)
                               ctcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGtt                                                          <  2:5131602/91‑1 (MQ=255)
                                gcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTtg                                                         >  1:5358327/1‑93 (MQ=255)
                                  ttCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTtgtg                                                       <  1:8567925/93‑1 (MQ=255)
                                  ttCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTtgtg                                                       >  2:2818409/1‑93 (MQ=255)
                                   tCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTgg                                                      >  1:2994632/1‑93 (MQ=255)
                                   tCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTgg                                                      <  2:960756/93‑1 (MQ=255)
                                   tCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTgg                                                      <  2:4353705/93‑1 (MQ=255)
                                    cGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGa                                                     <  1:7097595/93‑1 (MQ=255)
                                     gTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCTCGCGTCCTGTTTGGGTTGTGGaa                                                    >  2:7771985/1‑92 (MQ=255)
                                      tAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATc                                                   <  2:240469/93‑1 (MQ=255)
                                      tAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATc                                                   >  2:396873/1‑93 (MQ=255)
                                      tAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATc                                                   >  2:6340726/1‑93 (MQ=255)
                                       aaCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCa                                                  >  1:9925805/1‑93 (MQ=255)
                                       aaCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCa                                                  >  1:8883113/1‑93 (MQ=255)
                                       aaCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCa                                                  <  1:8178174/93‑1 (MQ=255)
                                       aaCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCa                                                  <  2:5938678/93‑1 (MQ=255)
                                        aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAg                                                 <  2:7955814/93‑1 (MQ=255)
                                        aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAg                                                 <  2:2593354/93‑1 (MQ=255)
                                        aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAg                                                 >  1:4755999/1‑93 (MQ=255)
                                        aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAg                                                 <  1:3615945/93‑1 (MQ=255)
                                         ccAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAgg                                                <  2:1078358/93‑1 (MQ=255)
                                          cAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGt                                               >  2:6223129/1‑93 (MQ=255)
                                          cAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGt                                               >  2:1157620/1‑93 (MQ=255)
                                           aaaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTa                                              >  2:4694343/1‑93 (MQ=255)
                                            aaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTAt                                             >  1:6834639/1‑93 (MQ=255)
                                               aaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCg                                          >  2:3090038/1‑93 (MQ=255)
                                               aaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCg                                          >  2:6688573/1‑93 (MQ=255)
                                               aaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCg                                          <  1:8411023/93‑1 (MQ=255)
                                                aGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGa                                         <  2:9214724/93‑1 (MQ=255)
                                                 gCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGAt                                        >  1:1264031/1‑93 (MQ=255)
                                                 gCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGAt                                        <  2:2030698/93‑1 (MQ=255)
                                                 gCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGAt                                        <  2:7420164/93‑1 (MQ=255)
                                                  cGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATc                                       <  1:2656507/93‑1 (MQ=255)
                                                  cGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATc                                       <  2:3848100/93‑1 (MQ=255)
                                                    tttGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTg                                     >  1:724438/1‑93 (MQ=255)
                                                    tttGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTg                                     <  2:4895731/93‑1 (MQ=255)
                                                     ttGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGt                                    <  2:9000222/93‑1 (MQ=255)
                                                      tGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGtt                                   <  2:952473/93‑1 (MQ=255)
                                                       gCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTc                                  <  2:2133290/93‑1 (MQ=255)
                                                       gCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTc                                  >  1:1142730/1‑93 (MQ=255)
                                                         ttACGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGt                                <  2:5625776/93‑1 (MQ=255)
                                                            cGGCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTaa                             <  2:3670099/93‑1 (MQ=255)
                                                             ggCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAAc                            >  1:1843344/1‑93 (MQ=255)
                                                              gCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCGTGTATTCGATCTGTTCGTTAAcc                           >  2:5590702/1‑93 (MQ=255)
                                                              gCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcc                           <  1:2093684/93‑1 (MQ=255)
                                                              gCCTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcc                           >  2:4267616/1‑93 (MQ=255)
                                                               ccTGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTCACCt                          >  2:7057340/1‑93 (MQ=255)
                                                                 tGAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCg                        <  1:241588/93‑1 (MQ=255)
                                                                  gAAGCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGa                       <  2:2466377/93‑1 (MQ=255)
                                                                     gCCGATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGAcgc                    <  1:1482248/93‑1 (MQ=255)
                                                                        gATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgac                 >  1:5525530/1‑93 (MQ=255)
                                                                        gATTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgac                 <  1:4399564/93‑1 (MQ=255)
                                                                         aTTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgacc                <  2:5882752/93‑1 (MQ=255)
                                                                         aTTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgacc                >  2:9113087/1‑93 (MQ=255)
                                                                         aTTGTTGTTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgacc                >  1:3827916/1‑93 (MQ=255)
                                                                           tgttgtTATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccga              <  1:7813972/93‑1 (MQ=255)
                                                                            gttgttATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccgac             <  1:9601371/93‑1 (MQ=255)
                                                                              tgttATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACcgacga           <  2:1888471/93‑1 (MQ=255)
                                                                              tgttATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACcgacga           <  1:1081947/93‑1 (MQ=255)
                                                                               gttATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGag          >  1:1646573/1‑93 (MQ=255)
                                                                                ttATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATTAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGagc         >  1:1369120/1‑93 (MQ=255)
                                                                                ttATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGagc         <  1:977720/93‑1 (MQ=255)
                                                                                ttATCAACAAGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGagc         >  1:6386978/1‑93 (MQ=255)
                                                                                     aacaaGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTg    >  1:4974555/1‑93 (MQ=255)
                                                                                      acaaGGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGATGCGACCGACGAGCCGCTgg   >  1:6033669/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGA  >  NC_000913/4056283‑4056462

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