Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,853,166 C→T 6.1% A272T (GCA→ACA)  ydjE ← putative MFS sugar transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,853,1660CT 6.1% 203.2 / 2.6 65A272T (GCA→ACA) ydjEputative MFS sugar transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base C (34/27);  new base T (2/2);  total (36/29)
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.77e-01

AATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAATTTGC  >  NC_000913/1853020‑1853314
                                                                                                                                                  |                                                                                                                                                    
aataatTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCaa                                                                                                                                                    >  6:205222/1‑149 (MQ=255)
    attaatGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTa                                                                                                                                                 >  2:243663/1‑148 (MQ=255)
     ttaatGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTaa                                                                                                                                                >  6:171049/1‑148 (MQ=255)
        atGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACa                                                                                                                                            <  8:85848/149‑1 (MQ=255)
          gCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAg                                                                                                                                           >  1:395837/1‑148 (MQ=255)
           ccGCAATAAATATTCCTACCGGAGCGCCATTCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGt                                                                                                                                          >  5:388120/1‑148 (MQ=255)
            cGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTaa                                                                                                                                        <  2:395837/149‑1 (MQ=255)
              cAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATc                                                                                                                                      <  8:329438/149‑1 (MQ=255)
               aataaataTTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATcg                                                                                                                                     <  5:205222/149‑1 (MQ=255)
                   aataTTCCTACCGGAGCGCCATTCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGAc                                                                                                                                 <  5:200775/149‑1 (MQ=255)
                    atatTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGAc                                                                                                                                 >  4:152044/1‑148 (MQ=255)
                      atTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTaa                                                                                                                              >  4:314143/1‑149 (MQ=255)
                        tCCTACCGGAGCGCCTCTCATAATAACATCGGGCATTAATAATGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATTCGGTTAGGGTATAAACCGATATGTTCATTGCAATTTAAACAGTAATCGCGACTAAt                                                                                                                             <  2:216503/148‑1 (MQ=255)
                           tACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCaa                                                                                                                                                    <  2:49818/122‑1 (MQ=255)
                           tACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCaa                                                                                                                                                    >  1:49818/1‑122 (MQ=255)
                                  gcgcCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGGCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTaa                                                                                                                  <  6:287086/149‑1 (MQ=255)
                                     ccAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGc                                                                                                                                                      >  4:8975/1‑110 (MQ=255)
                                     ccAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGc                                                                                                                                                      <  3:8975/110‑1 (MQ=255)
                                       aaTCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACAtt                                                                                                              >  7:186968/1‑148 (MQ=255)
                                        aTCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTc                                                                                                            <  3:152044/149‑1 (MQ=255)
                                           ataataACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCTGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGAAATCGCGACTAATGTACGTCGTAACATTTCAcc                                                                                                         >  3:289726/1‑149 (MQ=255)
                                           ataataACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCAc                                                                                                          <  3:365301/148‑1 (MQ=255)
                                              ataaCAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACATATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCtt                                                                                                       >  8:379133/1‑148 (MQ=255)
                                                      ggTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGa                                                                                              <  3:25016/149‑1 (MQ=255)
                                                       gTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCAc                                                                                                          <  7:134431/136‑1 (MQ=255)
                                                       gTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCAc                                                                                                          >  8:134431/1‑136 (MQ=255)
                                                         cATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTa                                                                                                                          <  6:351209/118‑1 (MQ=255)
                                                         cATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTa                                                                                                                          >  5:351209/1‑118 (MQ=255)
                                                              atatTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGa                                                                                        >  3:237028/1‑147 (MQ=255)
                                                                atTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAAGATTTCACCTTTAAACAGGAGCCAGAAAg                                                                                     >  1:133029/1‑148 (MQ=255)
                                                                atTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAg                                                                                     >  7:305651/1‑148 (MQ=255)
                                                                atTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAg                                                                                     >  8:191748/1‑148 (MQ=255)
                                                                           gACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTa                                                                         >  7:264854/1‑149 (MQ=255)
                                                                               tCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACtt                                                                      >  2:345528/1‑148 (MQ=255)
                                                                                cAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTa                                                                    >  4:300235/1‑149 (MQ=255)
                                                                                 aaTGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTaa                                                                   >  7:14632/1‑149 (MQ=255)
                                                                                 aaTGCCGGAGTTAACAAATATGGGCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTa                                                                    <  4:6367/148‑1 (MQ=255)
                                                                                       ggAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGc                                                             >  3:183311/1‑149 (MQ=255)
                                                                                         aGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTg                                                           <  3:314143/149‑1 (MQ=255)
                                                                                                     tGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATAc                                               >  8:351668/1‑149 (MQ=255)
                                                                                                       gTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACg                                              >  6:46454/1‑148 (MQ=255)
                                                                                                       gTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCAACTTTAAAAAGGAGCCAGAAAGTACACTTGATTACTTTAACCTTGCTGTTGCTATtataag                                              >  6:40907/1‑146 (MQ=255)
                                                                                                          ggTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAg                                           >  6:180101/1‑148 (MQ=255)
                                                                                                          ggTATCAATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAg                                           >  8:234735/1‑148 (MQ=255)
                                                                                                                cATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCgg                                    <  4:289726/149‑1 (MQ=255)
                                                                                                                 aTACGGTGATGGTATAAAGCGAAATGTTCATTGTAATTAAAACAGTAATCGCGAc                                                                                                                                 >  5:37873/1‑55 (MQ=255)
                                                                                                                 aTACGGTGATGGTATAAAGCGAAATGTTCATTGTAATTAAAACAGTAATCGCGAc                                                                                                                                 <  6:37873/55‑1 (MQ=255)
                                                                                                                 aTACGGTGATGGTATAAAGCGAAATGTTCATTGTAATTAAAACAGTAATCGCGAc                                                                                                                                 <  6:32964/55‑1 (MQ=255)
                                                                                                                 aTACGGTGATGGTATAAAGCGAAATGTTCATTGTAATTAAAACAGTAATCGCGAc                                                                                                                                 >  5:32964/1‑55 (MQ=255)
                                                                                                                  tACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCgggg                                  <  5:171049/149‑1 (MQ=255)
                                                                                                                    cGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTAGCTCTGATACGAAGTAAGCggggg                                 <  1:119981/148‑1 (MQ=255)
                                                                                                                             tataAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATa                       >  8:59662/1‑149 (MQ=255)
                                                                                                                             tataAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCCACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTTCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATa                       <  8:186968/149‑1 (MQ=255)
                                                                                                                               taAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATAct                     >  5:307513/1‑149 (MQ=255)
                                                                                                                               taAAGCGAAACGTTCATTGCAATTAAAACCGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATAct                     <  8:200467/149‑1 (MQ=255)
                                                                                                                                  aGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACtctt                  <  6:388120/149‑1 (MQ=255)
                                                                                                                                   gCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACtcttc                 >  5:124932/1‑149 (MQ=255)
                                                                                                                                          gTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACtctt                  >  3:244684/1‑141 (MQ=255)
                                                                                                                                          gTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACtctt                  <  4:244684/141‑1 (MQ=255)
                                                                                                                                           ttCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGTGGTAAACGAATACTCTTCTCTCttt          <  6:307513/148‑1 (MQ=255)
                                                                                                                                              aTTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCCCGTTTAAAAGGGAGCAAGAAAGTACCGTTGTTTAGTTTAACATTGCTGTTGCTCTGATACGAATTAAGCGGGGGTAAAGGAATAATCTTCTCTCTTTAAAt      >  3:251907/1‑149 (MQ=255)
                                                                                                                                              aTTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCaa       <  1:243663/148‑1 (MQ=255)
                                                                                                                                              aTTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAAt      <  4:183311/149‑1 (MQ=255)
                                                                                                                                              aTTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAAt      <  8:264854/149‑1 (MQ=255)
                                                                                                                                               ttGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAAtt     <  8:305651/149‑1 (MQ=255)
                                                                                                                                                tGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAAttt    >  5:190232/1‑149 (MQ=255)
                                                                                                                                                 gCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAATTtg   >  3:190490/1‑149 (MQ=255)
                                                                                                                                                  cAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAATTtg   <  2:360131/148‑1 (MQ=255)
                                                                                                                                                  cAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGAGGGGGTAAACGAATACTCTTATCTCTTTCAAGTtgc  >  3:228758/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATACGGTGATGGTATAAAGCGAAATGTTCATTGCAATTAAAACAGTAATCGCGACTAATGTACGTCGTAACATTTCACCTTTAAACAGGAGCCAGAAAGTACCCTTGATTACTTTAACCTTGCTGTTGCTCTGATACGAAGTAAGCGGGGGTAAACGAATACTCTTCTCTCTTTCAATTTGC  >  NC_000913/1853020‑1853314

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: