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)
aagaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCC‑CCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttc                                                                                                                    >  2:1117691/1‑138 (MQ=255)
  gaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   >  2:1637365/1‑137 (MQ=255)
  gaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   <  2:60037/139‑3 (MQ=255)
      aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAggg                                                                                                              <  2:1373062/139‑1 (MQ=255)
          aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTc                                                                                                          >  1:1184987/1‑139 (MQ=255)
          aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTc                                                                                                          >  1:1729972/1‑139 (MQ=255)
            cGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                        >  2:527029/1‑139 (MQ=255)
             gTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc                                                                                                       >  1:1823296/1‑139 (MQ=255)
              ttGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCC‑CCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACaa                                                                                                     <  1:1117691/139‑1 (MQ=255)
               tGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACaa                                                                                                     <  2:764486/139‑1 (MQ=255)
                   gTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCtt                                                                                                 <  2:1729972/139‑1 (MQ=255)
                        aaaCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAgg                                                                                                               <  1:620485/120‑1 (MQ=255)
                        aaaCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAgg                                                                                                               >  2:620485/1‑120 (MQ=255)
                        aaaCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGAt                                                                                            >  2:1890962/1‑139 (MQ=255)
                           cGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                                                                         >  1:359018/1‑139 (MQ=255)
                              cATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAggg                                                                                      >  1:73043/1‑139 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc                                                                                                       >  1:467696/1‑121 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc                                                                                                       <  2:467696/121‑1 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGc                                                                                     <  2:1823296/139‑1 (MQ=255)
                               aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGc                                                                                     >  1:1236549/1‑139 (MQ=255)
                                 cGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGACCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGa                                                                                   >  1:1494411/1‑139 (MQ=255)
                                  gTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                                                                                  <  2:359018/139‑1 (MQ=255)
                                         cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           <  2:1700118/139‑1 (MQ=255)
                                                 gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   <  1:1997286/92‑3 (MQ=255)
                                                 gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   >  2:1997286/1‑90 (MQ=255)
                                                 gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg                                                                   >  1:690775/1‑139 (MQ=255)
                                                  tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  2:73043/139‑1 (MQ=255)
                                                  tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  1:527029/139‑1 (MQ=255)
                                                  tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  2:2032327/139‑1 (MQ=255)
                                                   ttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 <  1:1890962/139‑1 (MQ=255)
                                                   ttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 <  1:876409/139‑1 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   >  2:115810/1‑85 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttca                                                                                                                   <  1:115810/87‑3 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc                                                                                                       <  2:522096/98‑1 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc                                                                                                       >  1:522096/1‑98 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg                                                              >  2:1371784/1‑139 (MQ=255)
                                                      cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg                                                              >  2:1373243/1‑139 (MQ=255)
                                                       gcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTtctgttctgt                                                             <  1:1243309/139‑1 (MQ=255)
                                                            cccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAt                                                        <  1:1637365/139‑1 (MQ=255)
                                                               ggCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgtg                                                     <  1:1373243/139‑1 (MQ=255)
                                                                gCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATgtgt                                                    <  2:1184987/139‑1 (MQ=255)
                                                                gCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATgtgt                                                    <  2:593952/139‑1 (MQ=255)
                                                                  tACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGttt                                                  <  2:690775/139‑1 (MQ=255)
                                                                     gTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 >  1:859249/1‑121 (MQ=255)
                                                                     gTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 <  2:859249/121‑1 (MQ=255)
                                                                      tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                        <  1:1445462/81‑1 (MQ=255)
                                                                      tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                        >  2:1445462/1‑81 (MQ=255)
                                                                                ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           <  1:1516599/100‑1 (MQ=255)
                                                                                ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           >  2:1516599/1‑100 (MQ=255)
                                                                                 cAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt                                   <  2:1236549/139‑1 (MQ=255)
                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTACATGTGtt                                                   <  1:1345842/114‑1 (MQ=255)
                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCa                          >  2:1444423/1‑139 (MQ=255)
                                                                                          cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑CTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                   >  2:1345842/1‑114 (MQ=255)
                                                                                               gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                       >  2:1843517/1‑121 (MQ=255)
                                                                                               gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                       <  1:1843517/121‑1 (MQ=255)
                                                                                                    aGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCTAAGTCGGAAAACttctgtt                                                                 <  1:1519788/90‑1 (MQ=255)
                                                                                                    aGTGCCTGAATGGGCGTACAGGTATG‑AACCCCTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                 >  2:1519788/1‑90 (MQ=255)
                                                                                                      tGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTACATGTGTTTTGCTCa                                            >  1:1798168/1‑109 (MQ=255)
                                                                                                      tGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCa                                            <  2:1798168/109‑1 (MQ=255)
                                                                                                       gCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCt             <  2:1494411/139‑1 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  2:1719733/82‑1 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  >  1:1719733/1‑82 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTa                  >  1:870565/1‑130 (MQ=255)
                                                                                                           gAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTa                  <  2:870565/130‑1 (MQ=255)
                                                                                                            aaTGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGTAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTgg                                    <  1:584760/111‑1 (MQ=255)
                                                                                                            aaTGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTgg                                    >  2:584760/1‑111 (MQ=255)
                                                                                                            aaTGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACg        <  1:1899232/139‑1 (MQ=255)
                                                                                                                ggCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                         >  1:735671/1‑102 (MQ=255)
                                                                                                                ggCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                         <  2:735671/102‑1 (MQ=255)
                                                                                                                  cGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCg  >  1:489416/1‑139 (MQ=255)
                                                                                                                    tACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtgt                                      >  2:44549/1‑101 (MQ=255)
                                                                                                                    tACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtgt                                      <  1:44549/101‑1 (MQ=255)
                                                                                                                          tATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           >  1:1248442/1‑58 (MQ=255)
                                                                                                                          tATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                           <  2:1248442/58‑1 (MQ=255)
                                                                                                                          tATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  <  1:1758412/67‑1 (MQ=255)
                                                                                                                          tATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                  >  2:1758412/1‑67 (MQ=255)
                                                                                                                           aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAatat                            >  2:642369/1‑104 (MQ=255)
                                                                                                                           aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAatat                            <  1:642369/104‑1 (MQ=255)
                                                                                                                           aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                   >  2:1945358/1‑113 (MQ=255)
                                                                                                                           aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                   <  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: ‑