Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,965,557 (T)8→7 coding (74/1266 nt) rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,5500T.98.6% 295.7 / ‑2.3 73coding (81/1266 nt)rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (32/40);  minor base C (1/0);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

AAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCG  >  NC_000913/3965418‑3965671
                                                                                                                                    |                                                                                                                         
aagaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttttt                                                                                                                     >  1:60037/1‑139 (MQ=255)
aagaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttc                                                                                                                    >  2:1117691/1‑138 (MQ=255)
  gaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   >  2:1637365/1‑137 (MQ=255)
  gaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   <  2:60037/139‑3 (MQ=255)
      aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAggg                                                                                                              <  2:1373062/139‑1 (MQ=255)
          aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTc                                                                                                          >  1:1184987/1‑139 (MQ=255)
          aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTc                                                                                                          >  1:1729972/1‑139 (MQ=255)
            cGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTa                                                                                                        >  2:527029/1‑139 (MQ=255)
             gTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTAc                                                                                                       >  1:1823296/1‑139 (MQ=255)
              ttGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACaa                                                                                                     <  1:1117691/139‑1 (MQ=255)
               tGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACaa                                                                                                     <  2:764486/139‑1 (MQ=255)
                   gTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCtt                                                                                                 <  2:1729972/139‑1 (MQ=255)
                        aaaCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAgg                                                                                                               <  1:620485/120‑1 (MQ=255)
                        aaaCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAgg                                                                                                               >  2:620485/1‑120 (MQ=255)
                        aaaCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGAt                                                                                            >  2:1890962/1‑139 (MQ=255)
                           cGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                                                                         >  1:359018/1‑139 (MQ=255)
                              cATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAggg                                                                                      >  1:73043/1‑139 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTAc                                                                                                       >  1:467696/1‑121 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTAc                                                                                                       <  2:467696/121‑1 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGc                                                                                     <  2:1823296/139‑1 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGc                                                                                     >  1:1236549/1‑139 (MQ=255)
                                 cGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGACCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGa                                                                                   >  1:1494411/1‑139 (MQ=255)
                                  gTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                                                                                  <  2:359018/139‑1 (MQ=255)
                                         cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           <  2:1700118/139‑1 (MQ=255)
                                                 gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   <  1:1997286/92‑3 (MQ=255)
                                                 gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   >  2:1997286/1‑90 (MQ=255)
                                                 gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg                                                                   >  1:690775/1‑139 (MQ=255)
                                                  tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  2:73043/139‑1 (MQ=255)
                                                  tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  1:527029/139‑1 (MQ=255)
                                                  tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  2:2032327/139‑1 (MQ=255)
                                                   ttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 <  1:1890962/139‑1 (MQ=255)
                                                   ttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 <  1:876409/139‑1 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   >  2:115810/1‑85 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   <  1:115810/87‑3 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTAc                                                                                                       <  2:522096/98‑1 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTAc                                                                                                       >  1:522096/1‑98 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg                                                              >  2:1371784/1‑139 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg                                                              >  2:1373243/1‑139 (MQ=255)
                                                       gcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTtctgttctgt                                                             <  1:1243309/139‑1 (MQ=255)
                                                            cccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAt                                                        <  1:1637365/139‑1 (MQ=255)
                                                               ggCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgtg                                                     <  1:1373243/139‑1 (MQ=255)
                                                                gCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATgtgt                                                    <  2:1184987/139‑1 (MQ=255)
                                                                gCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATgtgt                                                    <  2:593952/139‑1 (MQ=255)
                                                                  tACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGttt                                                  <  2:690775/139‑1 (MQ=255)
                                                                     gTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 >  1:859249/1‑121 (MQ=255)
                                                                     gTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 <  2:859249/121‑1 (MQ=255)
                                                                      tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTa                                                                                                        <  1:1445462/81‑1 (MQ=255)
                                                                      tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTa                                                                                                        >  2:1445462/1‑81 (MQ=255)
                                                                                ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           <  1:1516599/100‑1 (MQ=255)
                                                                                ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           >  2:1516599/1‑100 (MQ=255)
                                                                                 cAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt                                   <  2:1236549/139‑1 (MQ=255)
                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTACATGTGtt                                                   <  1:1345842/114‑1 (MQ=255)
                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCa                          >  2:1444423/1‑139 (MQ=255)
                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCCTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                   >  2:1345842/1‑114 (MQ=255)
                                                                                               gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                       >  2:1843517/1‑121 (MQ=255)
                                                                                               gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                       <  1:1843517/121‑1 (MQ=255)
                                                                                                    aGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCTAAGTCGGAAAACttctgtt                                                                 <  1:1519788/90‑1 (MQ=255)
                                                                                                    aGTGCCTGAATGGGCGTACAGGTATGAACCCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 >  2:1519788/1‑90 (MQ=255)
                                                                                                      tGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTACATGTGTTTTGCTCa                                            >  1:1798168/1‑109 (MQ=255)
                                                                                                      tGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCa                                            <  2:1798168/109‑1 (MQ=255)
                                                                                                       gCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCt             <  2:1494411/139‑1 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  2:1719733/82‑1 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  >  1:1719733/1‑82 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTa                  >  1:870565/1‑130 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTa                  <  2:870565/130‑1 (MQ=255)
                                                                                                            aaTGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGTAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTgg                                    <  1:584760/111‑1 (MQ=255)
                                                                                                            aaTGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTgg                                    >  2:584760/1‑111 (MQ=255)
                                                                                                            aaTGGGCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACg        <  1:1899232/139‑1 (MQ=255)
                                                                                                                ggCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                         >  1:735671/1‑102 (MQ=255)
                                                                                                                ggCGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                         <  2:735671/102‑1 (MQ=255)
                                                                                                                  cGTACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCg  >  1:489416/1‑139 (MQ=255)
                                                                                                                    tACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtgt                                      >  2:44549/1‑101 (MQ=255)
                                                                                                                    tACAGTTATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtgt                                      <  1:44549/101‑1 (MQ=255)
                                                                                                                          tATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           >  1:1248442/1‑58 (MQ=255)
                                                                                                                          tATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           <  2:1248442/58‑1 (MQ=255)
                                                                                                                          tATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  1:1758412/67‑1 (MQ=255)
                                                                                                                          tATGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  >  2:1758412/1‑67 (MQ=255)
                                                                                                                           aTGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAatat                            >  2:642369/1‑104 (MQ=255)
                                                                                                                           aTGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAatat                            <  1:642369/104‑1 (MQ=255)
                                                                                                                           aTGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                   >  2:1945358/1‑113 (MQ=255)
                                                                                                                           aTGAAACCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                   <  1:1945358/113‑1 (MQ=255)
                                                                                                                                    |                                                                                                                         
AAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCG  >  NC_000913/3965418‑3965671

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