Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,114,362 T→C 6.1% R377G (AGA→GGA)  mdtG ← putative drug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,114,3620TC 6.1% 200.9 / 3.4 65R377G (AGA→GGA) mdtGputative drug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

ATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGT  >  NC_000913/1114223‑1114507
                                                                                                                                           |                                                                                                                                                 
aTAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGt                                                                                                                                           <  2:327499/148‑1 (MQ=255)
            ccGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGaaa                                                                                                                               >  3:68942/1‑148 (MQ=255)
             cGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGaa                                                                                                                                <  3:315544/146‑1 (MQ=255)
                 tttGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAACAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCgct                                                                                                                          >  3:82250/1‑148 (MQ=255)
                         tATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCAt                                                                                                                  >  5:42504/1‑148 (MQ=255)
                         tATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCACCCAt                                                                                                                  >  7:33085/1‑148 (MQ=255)
                           tGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGGTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCa                                                                                                                <  6:342820/148‑1 (MQ=255)
                            gAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGt                                                                                                                                           >  7:114627/1‑120 (MQ=255)
                            gAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGt                                                                                                                                           <  8:114627/120‑1 (MQ=255)
                              aaGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATg                                                                                                             <  2:311695/148‑1 (MQ=255)
                                  cGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACTGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTcc                                                                                                         <  3:85052/148‑1 (MQ=255)
                                      aaaTCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTa                                                                                                     <  8:112589/148‑1 (MQ=255)
                                         tCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCTTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACg                                                                                                  <  6:234754/148‑1 (MQ=255)
                                         tCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACg                                                                                                  <  6:241060/148‑1 (MQ=255)
                                         tCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACg                                                                                                  >  8:173628/1‑148 (MQ=255)
                                                gATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCaa                                                                                           <  6:221068/148‑1 (MQ=255)
                                                      tGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGCAATACCGCTCTGAACCCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCAc                                                                                     >  2:185408/1‑148 (MQ=255)
                                                              tCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGAt                                                                             >  3:84313/1‑148 (MQ=255)
                                                              tCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGAt                                                                             >  8:45322/1‑148 (MQ=255)
                                                              tCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGAt                                                                             >  2:393884/1‑148 (MQ=255)
                                                                gacgacgacGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTg                                                                           >  8:230009/1‑148 (MQ=255)
                                                                 acgacgacgTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTgg                                                                          <  7:436813/148‑1 (MQ=255)
                                                                   gacgacgTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGtt                                                                        <  7:45322/148‑1 (MQ=255)
                                                                   gacgacgTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGtt                                                                        <  7:173628/148‑1 (MQ=255)
                                                                    acgacgTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTta                                                                       >  5:166547/1‑148 (MQ=255)
                                                                            gACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGa                                                                <  4:84313/147‑1 (MQ=255)
                                                                            gACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACAGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTACCGTAGCCAATATCACGAAACGATTGGGTATAGCTGaa                                                               >  1:66409/1‑148 (MQ=255)
                                                                              cTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTg                                                                 >  3:78021/1‑144 (MQ=255)
                                                                              cTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCCGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGaa                                                                                  <  2:41971/127‑1 (MQ=255)
                                                                              cTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCCGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGaa                                                                                  >  1:41971/1‑127 (MQ=255)
                                                                              cTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCCGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGaa                                                                                  >  1:45753/1‑127 (MQ=255)
                                                                              cTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCCGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGaa                                                                                  <  2:45753/127‑1 (MQ=255)
                                                                                gTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATc                                                                                       <  4:400362/120‑1 (MQ=255)
                                                                                gTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATc                                                                                       >  3:400362/1‑120 (MQ=255)
                                                                                 ttCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGAt                                                            <  6:312535/146‑1 (MQ=255)
                                                                                 ttCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGAt                                                            >  5:312535/1‑146 (MQ=255)
                                                                                  tCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTgg                                                                          >  1:201957/1‑131 (MQ=255)
                                                                                  tCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTgg                                                                          <  2:201957/131‑1 (MQ=255)
                                                                                   ccATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACg                                                         <  6:42504/147‑1 (MQ=255)
                                                                                    cATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGcc                                                       >  3:40170/1‑148 (MQ=255)
                                                                                          tAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGa                                                 >  2:195073/1‑148 (MQ=255)
                                                                                            gACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATc                                               >  3:335963/1‑148 (MQ=255)
                                                                                                   ttGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGc                                        >  4:129218/1‑148 (MQ=255)
                                                                                                   ttGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGc                                        <  7:230009/148‑1 (MQ=255)
                                                                                                    tGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCt                                       >  4:125385/1‑148 (MQ=255)
                                                                                                     gAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCt                                       <  1:195073/147‑1 (MQ=255)
                                                                                                       aTAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAATCGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGa                                    <  8:33085/148‑1 (MQ=255)
                                                                                                          acacTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTg                                             >  1:208453/1‑136 (MQ=255)
                                                                                                          acacTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTg                                             <  2:208453/136‑1 (MQ=255)
                                                                                                            acTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTg                                >  7:247258/1‑147 (MQ=255)
                                                                                                             cTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTa                              <  1:393884/148‑1 (MQ=255)
                                                                                                              tACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTtat                                                                      >  3:88865/1‑107 (MQ=255)
                                                                                                              tACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTtat                                                                      <  4:88865/107‑1 (MQ=255)
                                                                                                               ccGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAACTCTCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATCCGCCCGGCGATCTGGTTGCTCGAGTTGTaaa                            <  2:66409/147‑1 (MQ=255)
                                                                                                                cGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTaa                             <  4:82250/146‑1 (MQ=255)
                                                                                                                  ccAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAAcc                          >  7:316791/1‑147 (MQ=255)
                                                                                                                  ccAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCAGGAGATCTGGTTGCTCTAGTTGTAAAcc                          >  1:97496/1‑147 (MQ=255)
                                                                                                                    aGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGATCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAAc                       <  4:68942/148‑1 (MQ=255)
                                                                                                                            ggtgaAATACCTCGCTGAAACCTTAGTTCGCTGACATCGCTGCTAACATCAATGGTCCGGTAACGTTGCCAATTTCACGAAACGTTTGGTTATTGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCtt               <  2:97496/145‑2 (MQ=255)
                                                                                                                            cGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTg               <  8:316791/148‑1 (MQ=255)
                                                                                                                              agaAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTa             >  1:194204/1‑148 (MQ=255)
                                                                                                                                     ccGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTAcggcgg       <  8:247258/147‑1 (MQ=255)
                                                                                                                                     ccGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCggg      >  7:245928/1‑148 (MQ=255)
                                                                                                                                      cGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCgggg     <  3:125385/148‑1 (MQ=255)
                                                                                                                                       gCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCggg      >  6:159464/1‑146 (MQ=255)
                                                                                                                                         tctGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGagt  <  4:40170/148‑1 (MQ=255)
                                                                                                                                           tGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGagt  <  6:166547/146‑1 (MQ=255)
                                                                                                                                           tGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGagt  <  5:159464/146‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
ATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGT  >  NC_000913/1114223‑1114507

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: ‑