Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,451,714 C→T 5.4% R219R (CGC→CGT ytfR → putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,451,7140CT 5.4% 237.5 / 3.0 74R219R (CGC→CGTytfRputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base T (2/2);  total (38/36)
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.27e-01

AACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGACAAACCCGTTGCCGCGTTCAA  >  NC_000913/4451569‑4451860
                                                                                                                                                 |                                                                                                                                                  
aaCCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCaa                                                                                                                                                  <  3:159321/148‑1 (MQ=255)
 aCCCACCTCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACg                                                                                                                                                <  7:105613/149‑1 (MQ=255)
 aCCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACg                                                                                                                                                >  4:66884/1‑149 (MQ=255)
    cACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGGCCGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCTGCTCACCGTCTTACGCAAAgg                                                                                                                                               >  2:313758/1‑147 (MQ=255)
     aCCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAg                                                                                                                                            >  2:90749/1‑149 (MQ=255)
       cGCCAGTCTCGACACTCAGGAAGTAGAGTTACTGTTTGGCCTGATGCGTCAGCTTCGCGATCTCGGCGTCAGCCTGATCTTCGTCACTCATTTTCTCGATCAGGTTTATCAGGTTAGCGATCGGATCACCGTCTTACGTAACGGCAGtt                                                                                                                                          >  8:409909/1‑149 (MQ=255)
       cGCCAGTCTCGACACTCAGGAAGAAGAGTTACTGTTTGGCCTGATGCGTCAGCTTCGCGATCGCGGCGTCAGCCTCATCTTCGTCACTCTTTTTCTCGATCAGGTTTATCAGGATAGCGATCGGATCACCGTCTTACGTAACTGCAGtt                                                                                                                                          >  8:415401/1‑149 (MQ=255)
          cAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCg                                                                                                                                       <  5:165314/149‑1 (MQ=255)
            gTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGt                                                                                                                                      <  7:14375/148‑1 (MQ=255)
            gTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGt                                                                                                                                      >  8:14375/1‑148 (MQ=255)
            gTCTAGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCTGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTa                                                                                                                                     <  3:192651/149‑1 (MQ=255)
               tcGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGc                                                                                                                                  >  3:373747/1‑149 (MQ=255)
               tcGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGc                                                                                                                                  >  8:219764/1‑149 (MQ=255)
                 gACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAgg                                                                                                                                   <  2:352497/146‑1 (MQ=255)
                  acacCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTTTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGt                                                                                                                               >  8:348881/1‑149 (MQ=255)
                     ccccGGAAGTGGAGGTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGGCAGCGATCGGATCCCCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCgg                                                                                                                            <  1:313758/149‑1 (MQ=255)
                     cccAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTa                                                                                                                                     <  4:225769/140‑1 (MQ=255)
                     cccAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTa                                                                                                                                     >  3:225769/1‑140 (MQ=255)
                          gAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAAc                                                                                                                       >  3:19715/1‑149 (MQ=255)
                             ttGGAGTTACTGTTTGACCTTATGCGTCAGTTTCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGc                                                                                                                                              <  8:150897/122‑1 (MQ=255)
                             gTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGc                                                                                                                                              >  7:150897/1‑123 (MQ=255)
                              tGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGc                                                                                                                   >  7:272039/1‑149 (MQ=255)
                              tGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGAGGCGTCAGCCTGATTTTTGTAACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGc                                                                                                                   >  3:143547/1‑149 (MQ=255)
                                   ttACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGc                                                                                                               >  4:80967/1‑148 (MQ=255)
                                      cTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTAc                                                                                                            <  1:207099/148‑1 (MQ=255)
                                              ccTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGtt                                                                                                                                          >  7:350243/1‑110 (MQ=255)
                                              ccTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGtt                                                                                                                                          <  8:350243/110‑1 (MQ=255)
                                                  aTGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAAc                                                                                                <  4:19715/148‑1 (MQ=255)
                                                   tGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTg                                                                                              <  3:66884/149‑1 (MQ=255)
                                                       tCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACg                                                                                                                                                <  1:267681/95‑1 (MQ=255)
                                                       tCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACg                                                                                                                                                >  2:267681/1‑95 (MQ=255)
                                                       tCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTaa                                                                                          <  4:373747/149‑1 (MQ=255)
                                                          gTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCTATCGGATCACCGTCTTACGCAACGTCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTaaaaa                                                                                       <  7:348881/149‑1 (MQ=255)
                                                             gcgcGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAatg                                                                                     >  2:413417/1‑148 (MQ=255)
                                                                           tCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAg                                                                      >  7:199226/1‑149 (MQ=255)
                                                                            cAGCCTGATCTTCGTCACTCATTTTCTCGATCAGGTTTATCAGGTTAGCGATCGGATCACCGTCTTACGTAACGGCAGTTTTGTGGGCTGTCGGGAAACCCGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGTGaac                                                                     <  7:409909/149‑3 (MQ=255)
                                                                            cAGCCTGATCTTCGTCACTCATTTTCTCGATCAGGTTTATCAGGTTAGCGATCGGAGCACCGTCTTACGTAACGGCATTTTTGTGGGCTGTCGGGAAACCCGCGAGCTAGCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGTGaac                                                                     <  7:415401/149‑3 (MQ=255)
                                                                               ccTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTg                                                                   >  4:286004/1‑148 (MQ=255)
                                                                               ccTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTg                                                                   >  6:220093/1‑148 (MQ=255)
                                                                                cTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGa                                                                 <  8:220982/149‑1 (MQ=255)
                                                                                         gTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCAcgc                                                        >  5:310420/1‑149 (MQ=255)
                                                                                            actcactTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGaa                                                                                                 <  2:388431/105‑1 (MQ=255)
                                                                                            actcactTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGaa                                                                                                 >  1:388431/1‑105 (MQ=255)
                                                                                               cactTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTa                                                    >  5:252136/1‑147 (MQ=255)
                                                                                                     ctcGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTg                                             <  3:74553/148‑1 (MQ=255)
                                                                                                      tcGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCAc                                                          <  2:370319/134‑1 (MQ=255)
                                                                                                      tcGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCAc                                                          >  1:370319/1‑134 (MQ=255)
                                                                                                         aTCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCggg                                        >  4:324845/1‑149 (MQ=255)
                                                                                                           cAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGc                                       >  3:73596/1‑148 (MQ=255)
                                                                                                            aGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCg                                      >  1:337664/1‑148 (MQ=255)
                                                                                                              gTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGaa                                    >  5:289796/1‑148 (MQ=255)
                                                                                                              gTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAAc                                   <  3:286004/149‑1 (MQ=255)
                                                                                                               tCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACa                                  <  5:220093/149‑1 (MQ=255)
                                                                                                                 tATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAtt                                >  7:367949/1‑149 (MQ=255)
                                                                                                                     aGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAttgttg                            <  4:14405/149‑1 (MQ=255)
                                                                                                                     aGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAttgtt                             <  3:80967/148‑1 (MQ=255)
                                                                                                                     aGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAttgtt                             <  1:413417/148‑1 (MQ=255)
                                                                                                                      ggTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAttg                               >  1:120504/1‑145 (MQ=255)
                                                                                                                      ggTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAttg                               <  2:120504/145‑1 (MQ=255)
                                                                                                                        tCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGc                         <  4:73596/149‑1 (MQ=255)
                                                                                                                          aGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCg                                                                                                                                       >  1:329247/1‑37 (MQ=255)
                                                                                                                          aGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCg                                                                                                                                       <  2:329247/37‑1 (MQ=255)
                                                                                                                          aGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGa                       <  7:219764/149‑1 (MQ=255)
                                                                                                                           gCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGAc                      <  1:90749/149‑1 (MQ=255)
                                                                                                                             gATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGAc                      >  1:58718/1‑147 (MQ=255)
                                                                                                                             gATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGAc                      <  2:58718/147‑1 (MQ=255)
                                                                                                                             gATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAAATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGACa                     >  6:385007/1‑148 (MQ=255)
                                                                                                                                 ggATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGACAAAccc                >  7:233891/1‑149 (MQ=255)
                                                                                                                                      aCCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAttgtt                             <  6:400934/131‑1 (MQ=255)
                                                                                                                                      aCCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACAttgtt                             >  5:400934/1‑131 (MQ=255)
                                                                                                                                       ccGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGACAAACCCGTTGcc          <  3:324845/149‑1 (MQ=255)
                                                                                                                                       ccGACTAACCCAACGGCAGTTTCGTAGCCTCTCGGGAAACGTGCGAGCTACCGCAGATCCAACTCGTAAAAAAGATCCTGGCGCGCGACCTGGATCCCCACGCGCTACAGCGAGCCGGGCGAACATTGTTGAGCGACAAACCCGTTGcc          >  5:399788/1‑149 (MQ=255)
                                                                                                                                              aCGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGCCAAACCCGTTGCCGCGTTc    >  2:299902/1‑148 (MQ=255)
                                                                                                                                               cGCTACGGCTGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGTTACCGAGGCGCTACAGCGTGGCGGGCGATCATTGTTGAGCGACAATCCGGTTGCCGGGTTCaa  <  6:399788/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTCAGTTGCGCGATCGCGGCGTCAGCCTGATTTTTGTCACTCACTTTCTCGATCAGGTCTATCAGGTCAGCGATCGGATCACCGTCTTACGCAACGGCAGTTTCGTAGGCTGTCGGGAAACGTGCGAGCTACCGCAGATCGAACTGGTAAAAATGATGCTGGGGCGCGAGCTGGATACCCACGCGCTACAGCGTGCCGGGCGAACATTGTTGAGCGACAAACCCGTTGCCGCGTTCAA  >  NC_000913/4451569‑4451860

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: