Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,683,685 C→A 6.1% L268L (CTC→CTA)  rstB → sensory histidine kinase of RstAB two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,683,6850CA 6.1% 214.6 / 3.5 66L268L (CTC→CTA) rstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base A (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

GCTCACGAACTGCGAACACCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGAT  >  NC_000913/1683539‑1683829
                                                                                                                                                  |                                                                                                                                                
gCTCACGAACTGCGAACACCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCg                                                                                                                                                 <  7:161797/148‑1 (MQ=255)
           gCGAACACCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCAGCTGACTTATGCCCGACTCGATCGCCcaca                                                                                                                                       >  7:359479/1‑147 (MQ=255)
               acacCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACg                                                                                                                                  <  8:42314/148‑1 (MQ=255)
                    gTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCAAAAAACGAGCtt                                                                                                                             >  3:182918/1‑148 (MQ=255)
                         tGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACg                                                                                                                                  <  3:188724/138‑1 (MQ=255)
                         tGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACg                                                                                                                                  >  4:188724/1‑138 (MQ=255)
                          gcgcCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCt                                                                                                                        >  1:149955/1‑147 (MQ=255)
                                    tATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATAGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAAc                                                                                                                >  5:96013/1‑145 (MQ=255)
                                     aTCGACGTGAGATGAGCGATAACCTGAGCGCCGCAGAATCCCAGGCGTTGAATCGTGATCTCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTTACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGa                                                                                                            <  4:276908/148‑1 (MQ=255)
                                               gaTGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTtg                                                                                                  <  5:29043/148‑1 (MQ=255)
                                                        tAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATATTAGCGAACCAGCCCTGCCGTAGTGGCTGTc                                                                                          >  8:381083/1‑147 (MQ=255)
                                                        tAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTAGATCGCCCACaaaa                                                                                                                                    >  4:70288/1‑105 (MQ=255)
                                                        tAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTAGATCGCCCACaaaa                                                                                                                                    <  3:75346/105‑1 (MQ=255)
                                                        tAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTAGATCGCCCACaaaa                                                                                                                                    <  3:70288/105‑1 (MQ=255)
                                                        tAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTAGATCGCCCACaaaa                                                                                                                                    >  4:75346/1‑105 (MQ=255)
                                                           ccTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACg                                                                                      <  4:182918/148‑1 (MQ=255)
                                                             tGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCa                                                                                    >  4:281397/1‑148 (MQ=255)
                                                               aGCGCCGCCGAATCCCCGGCGTTGAATCGTGCTATCAGTCAACTTGAAGCTTTAATTGAAGCGCTGCTGACTTATGCCCGACTCGCTCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCAt                                                                                   <  6:22401/147‑1 (MQ=255)
                                                                      ccgAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGa                                                                           <  8:233669/148‑1 (MQ=255)
                                                                          aTCCCAGGCGTTGAATCGTGATATCAGTCAACATGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATAtt                                                                       >  2:362866/1‑148 (MQ=255)
                                                                               aGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAgg                                                                   >  8:30741/1‑147 (MQ=255)
                                                                                   ggtGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTa                                                              <  1:116597/146‑1 (MQ=255)
                                                                                   gTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTa                                                              <  2:106866/148‑1 (MQ=255)
                                                                                   gTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTa                                                              <  2:149955/148‑1 (MQ=255)
                                                                                     tGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAAc                                                            >  1:194932/1‑148 (MQ=255)
                                                                                       aaTCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTc                                                                                                                            >  3:153008/1‑82 (MQ=255)
                                                                                       aaTCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTc                                                                                                                            <  4:153008/82‑1 (MQ=255)
                                                                                         tCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGACCCAGCCCTGCCGGGGTGGCTGGCCACGCATCTGGCAGATATTCAGGCAGTACCGccc                                                        >  6:33329/1‑148 (MQ=255)
                                                                                              atatCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATa                                                    >  1:90977/1‑147 (MQ=255)
                                                                                                atCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATaaaa                                                 >  6:152229/1‑148 (MQ=255)
                                                                                                  cAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGt                                                                                                    >  7:377024/1‑95 (MQ=255)
                                                                                                  cAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGt                                                                                                    <  8:377024/95‑1 (MQ=255)
                                                                                                  cAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGc                                                                  >  5:125586/1‑129 (MQ=255)
                                                                                                  cAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGc                                                                  <  6:125586/129‑1 (MQ=255)
                                                                                                  cAGTCAACTAGAAGCTTTAATTGAAGAGCTGCTGACTTATGCACGACTCGATCGCCCACAAAACCAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCAGGCAGATATTCAGGCAGTAACGCCCGCTAAACCg                                               >  1:33371/1‑148 (MQ=255)
                                                                                                    gTCAACTTGAAGCTTTAATTGTAGTGGTGCTGACTTATGCCCGACTCGATCTCCCCCAAAACGATCTTGATCTTAGCGATCCAGACCTGCCGGTGTGGCTGTCATCGCATCTGGCTGATATTCAGGCAGTAACGCCCGTTAAAACGGt                                             <  2:33371/148‑1 (MQ=255)
                                                                                                    gTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGt                                             <  3:111881/148‑1 (MQ=255)
                                                                                                      cAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGt                                             >  2:83301/1‑146 (MQ=255)
                                                                                                      cAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGt                                             <  1:83301/146‑1 (MQ=255)
                                                                                                      cAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTAc                                           >  1:141588/1‑148 (MQ=255)
                                                                                                         cTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCATCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGa                                        <  2:194932/148‑1 (MQ=255)
                                                                                                           tGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGat                                                                          <  4:293357/112‑1 (MQ=255)
                                                                                                           tGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGAt                                       >  3:50738/1‑147 (MQ=255)
                                                                                                           tGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGAAGat                                                                          >  3:293357/1‑112 (MQ=255)
                                                                                                                  ttAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGAGAACCAGACCTGCCGTAGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGc                               >  6:207828/1‑148 (MQ=255)
                                                                                                                     aTTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCg                            <  1:362866/148‑1 (MQ=255)
                                                                                                                        gAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTg                          <  8:84426/147‑1 (MQ=255)
                                                                                                                         aaGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTTTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCa                        <  6:96013/148‑1 (MQ=255)
                                                                                                                         aaGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCa                        <  2:90977/148‑1 (MQ=255)
                                                                                                                          agagCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGc                                                                                                                               <  6:143569/44‑1 (MQ=255)
                                                                                                                          agagCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGc                                                                                                                               >  5:143569/1‑44 (MQ=255)
                                                                                                                          agagCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGc                         >  3:267528/1‑146 (MQ=255)
                                                                                                                             gctgctGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTc                                                                      <  7:108263/98‑1 (MQ=255)
                                                                                                                             gctgctGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTc                                                                      >  8:108263/1‑98 (MQ=255)
                                                                                                                                gctgACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                 >  2:113026/1‑148 (MQ=255)
                                                                                                                                 ctgACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                 >  4:353583/1‑147 (MQ=255)
                                                                                                                                 ctgAATTATGCCCGACTCGATCGCCCACAAAACGATCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCa                  <  5:22755/146‑1 (MQ=255)
                                                                                                                                   gACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGAt                                       >  7:10267/1‑123 (MQ=255)
                                                                                                                                   gACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGAt                                       <  8:10267/123‑1 (MQ=255)
                                                                                                                                   atgttATGCCCGACTCGATCGCCCATATAACGTGCTTCATCTTAGCGATCCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAtt                <  2:126151/143‑1 (MQ=255)
                                                                                                                                       tATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATgcgg          >  7:155711/1‑148 (MQ=255)
                                                                                                                                        aTGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATgcggc         >  2:405795/1‑148 (MQ=255)
                                                                                                                                         tGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTCTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAtaa               >  1:126151/1‑139 (MQ=255)
                                                                                                                                               aCTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat  <  8:359479/148‑1 (MQ=255)
                                                                                                                                               aCTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat  <  4:50738/148‑1 (MQ=255)
                                                                                                                                               aCTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa   >  6:131993/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCTCACGAACTGCGAACACCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGAGCTGCTGACTTATGCCCGACTCGATCGCCCACAAAACGAGCTTCATCTTAGCGAACCAGACCTGCCGTTGTGGCTGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGAT  >  NC_000913/1683539‑1683829

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑