Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 487,675 G→T 6.2% N209K (AAC→AAA)  pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947487,6750GT 6.2% 205.9 / 3.7 64N209K (AAC→AAA) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base T (2/2);  total (35/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 = 9.03e-01

GTGCTCCAGATATTTGGGGGTAAACAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCAC  >  NC_002947/487531‑487815
                                                                                                                                                |                                                                                                                                            
gTGCTCCAGATATTTGGGGGTAAACAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGa                                                                                                                                          <  3:42613/149‑1 (MQ=255)
    tCCAGATATTTGGGGGTAAACAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGTTTGAGGcc                                                                                                                                      >  5:46903/1‑149 (MQ=255)
    tCCAGATATTTGGGGGTAAACAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGTTTGAGGcc                                                                                                                                      >  5:51840/1‑149 (MQ=255)
    tCCAGATATTTGGGGGTAAACAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGcc                                                                                                                                      >  8:105328/1‑149 (MQ=255)
              tGGGGGTAAACAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAg                                                                                                                             <  4:35080/148‑1 (MQ=255)
                        cAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGt                                                                                                                  >  7:67160/1‑149 (MQ=255)
                        cAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGt                                                                                                                  >  7:263765/1‑149 (MQ=255)
                        cAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGt                                                                                                                  >  7:258988/1‑149 (MQ=255)
                         aGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGtt                                                                                                                 <  2:33574/149‑1 (MQ=255)
                          gggTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTg                                                                                                                <  8:263765/149‑1 (MQ=255)
                          gggTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTg                                                                                                                <  8:258988/149‑1 (MQ=255)
                           ggTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTgg                                                                                                               <  8:17875/149‑1 (MQ=255)
                            gTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGc                                                                                                              <  4:57230/149‑1 (MQ=255)
                               tCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAAt                                                                                                           <  8:68586/149‑1 (MQ=255)
                               tCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAAt                                                                                                           <  6:136214/149‑1 (MQ=255)
                                 ggccggcAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATg                                                                                                          <  7:105328/148‑1 (MQ=255)
                                  gccggcAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGc                                                                                                         >  8:103402/1‑148 (MQ=255)
                                  gccggcAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGc                                                                                                         >  4:308114/1‑148 (MQ=255)
                                      gcAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAAc                                                                                                    <  7:103402/149‑1 (MQ=255)
                                          cGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGt                                                                                                >  6:323112/1‑149 (MQ=255)
                                                        ccATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCgg                                                                                                                    >  1:273312/1‑115 (MQ=255)
                                                        ccATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCgg                                                                                                                    <  2:273312/115‑1 (MQ=255)
                                                         cATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCa                                                                                  >  5:79260/1‑148 (MQ=255)
                                                          aTGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGATCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATg                                                                                >  8:17202/1‑149 (MQ=255)
                                                             ccTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGTTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCa                                                                              <  6:51840/148‑1 (MQ=255)
                                                             ccTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGTTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCa                                                                              <  6:46903/148‑1 (MQ=255)
                                                                       gcAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGt                                                                    >  7:251804/1‑148 (MQ=255)
                                                                             gggCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCTTGCGGGGTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGt                                                             <  6:79674/149‑1 (MQ=255)
                                                                              ggCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGtt                                                            >  7:62249/1‑149 (MQ=255)
                                                                               gCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTc                                                           <  5:323112/149‑1 (MQ=255)
                                                                                         gCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGt                                                                                                >  5:191605/1‑102 (MQ=255)
                                                                                         gCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGt                                                                                                <  6:191605/102‑1 (MQ=255)
                                                                                         gCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGctct                                                 >  6:5426/1‑149 (MQ=255)
                                                                                             gatgatGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGc                                               <  1:214375/147‑1 (MQ=255)
                                                                                             gatgatGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCa                                              <  5:5426/148‑1 (MQ=255)
                                                                                             gatgatGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAg                                             <  8:62249/149‑1 (MQ=255)
                                                                                                 atgTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCTCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGa                                          <  3:260739/148‑1 (MQ=255)
                                                                                                  tgTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGa                                          >  4:223323/1‑147 (MQ=255)
                                                                                                  tgTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATa                                        <  8:67160/149‑1 (MQ=255)
                                                                                                   gTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAg                                       >  7:122941/1‑149 (MQ=255)
                                                                                                    tCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAg                                       <  8:122941/148‑1 (MQ=255)
                                                                                                      gATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCga                                    >  1:282423/1‑148 (MQ=255)
                                                                                                           cttcGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGAc                                                        >  4:141404/1‑124 (MQ=255)
                                                                                                            ttcGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCCTGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCaa                              >  5:235748/1‑149 (MQ=255)
                                                                                                            ttcGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGCCGCGCTCTGCAGTGATAGCGTCGGCaa                              >  8:119547/1‑149 (MQ=255)
                                                                                                             tcGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCaa                              <  7:17202/148‑1 (MQ=255)
                                                                                                                  gCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCgg                                 <  3:202156/140‑1 (MQ=255)
                                                                                                                  gCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCgg                                 >  4:202156/1‑140 (MQ=255)
                                                                                                                     cAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACg                       >  4:301793/1‑147 (MQ=255)
                                                                                                                     cAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACg                       >  1:87330/1‑147 (MQ=255)
                                                                                                                       ggTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATACAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGt                   <  6:79260/149‑1 (MQ=255)
                                                                                                                        gTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTg                  <  4:48051/149‑1 (MQ=255)
                                                                                                                         tGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTgg                 <  2:87330/149‑1 (MQ=255)
                                                                                                                             cgccCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAgg             >  5:205102/1‑149 (MQ=255)
                                                                                                                             cgccCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAgg             >  4:12491/1‑149 (MQ=255)
                                                                                                                             cgccCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAAGGCCAAACTAGTCGCGCATGTCGGAATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGGAGTGGCAgg             >  7:206690/1‑149 (MQ=255)
                                                                                                                             cgccCTCGCCCGAATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAgg             >  3:39900/1‑149 (MQ=255)
                                                                                                                              gccCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGt            <  5:60990/149‑1 (MQ=255)
                                                                                                                                  tCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGt        >  3:210982/1‑149 (MQ=255)
                                                                                                                                    gCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTgg      >  3:199338/1‑149 (MQ=255)
                                                                                                                                       cGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACTGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCa   >  7:211848/1‑149 (MQ=255)
                                                                                                                                        gCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCCCGCAGTGGCAGGTGCGGGGTCAc  >  8:299140/1‑149 (MQ=255)
                                                                                                                                        gCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCAc  >  3:199472/1‑149 (MQ=255)
                                                                                                                                          aTGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCAc  >  3:317136/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTGCTCCAGATATTTGGGGGTAAACAGGGTATCGGCCGGCAGCGGCCCGCGCAGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGTCAC  >  NC_002947/487531‑487815

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