Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,816 T→G intergenic (+53/+700) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8160TG98.3% 196.3 / ‑4.4 59intergenic (+53/+700)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (28/30);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

CAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740679‑4740947
                                                                                                                                         |                                                                                                                                      
cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCTTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑c                                                                                                                         >  3:321459/1‑149 (MQ=255)
      gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTc                                                                                                                   <  7:293204/149‑1 (MQ=255)
       tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTcg                                                                                                                  >  6:85527/1‑149 (MQ=255)
        aCAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTcgc                                                                                                                 >  5:83247/1‑149 (MQ=255)
         cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTcgcg                                                                                                                >  7:289156/1‑149 (MQ=255)
              tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCt                                                                                                            <  7:237163/148‑1 (MQ=255)
                 gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGcc                                                                                                        >  6:353562/1‑149 (MQ=255)
                    gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCcgc                                                                                                     <  5:353562/149‑1 (MQ=255)
                       cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTTCTGTGCCGG‑‑‑CCTCTTCGCGGGCAAGGCCGCGaa                                                                                                  <  7:329592/149‑1 (MQ=255)
                       cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGGCAAGGCCGCGaa                                                                                                  <  7:61241/149‑1 (MQ=255)
                       cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGaa                                                                                                  <  3:140930/149‑1 (MQ=255)
                       cccgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGaa                                                                                                  <  2:289205/147‑1 (MQ=255)
                         cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                <  4:134421/149‑1 (MQ=38)
                         cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                <  1:89529/149‑1 (MQ=255)
                          gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               <  8:119582/149‑1 (MQ=255)
                          gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTAAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               <  8:39835/149‑1 (MQ=255)
                           ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAgg                                                                                              <  5:164083/149‑1 (MQ=255)
                              ggtGTACACGGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGGTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCGTTGCCGGCGAAGAGGCCg                                                                                           <  2:138928/147‑1 (MQ=255)
                               cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTTCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                          >  6:193208/1‑149 (MQ=255)
                                tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGc                                                                                         >  6:272202/1‑149 (MQ=255)
                                 gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGGCAAGGCCGCGAagag                                                                                               <  7:175465/142‑1 (MQ=38)
                                 gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGGCAAGGCCGCGAagag                                                                                               >  8:175465/1‑142 (MQ=38)
                                   acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                        >  1:9100/1‑147 (MQ=255)
                                        ggCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGGACCTATTAGCGGCCTTGCCCGCGAAGAGGCCGGccaaga                                                                                    >  1:129132/1‑144 (MQ=255)
                                               aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCCGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                          >  8:91728/1‑149 (MQ=255)
                                                gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAAc                                                                         >  7:363516/1‑149 (MQ=255)
                                                gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                         >  7:271004/1‑149 (MQ=255)
                                                  aGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCt                                                                       <  5:85527/149‑1 (MQ=255)
                                                      cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCACTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGACACCACAAAACCACAgg                                                                   >  7:326078/1‑148 (MQ=255)
                                                       gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCCGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                                  <  7:91728/149‑1 (MQ=255)
                                                           tCGTTCCCCTGAAAGCCCGCAAGTACGCTGCAGTAGCCGAACGCACGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                              <  8:363516/149‑1 (MQ=255)
                                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTATTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCGCAGTTTTTGTccgccc                                                       >  8:57438/1‑149 (MQ=255)
                                                                         gACCGAATGTACGCTCCAGTAGCCTAACGCAAGACGGCTGCCCTCGGGCAGCTTTTTTTTTTTGTCCTGTGCCGG‑‑‑CCTCTGCGGGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                <  2:9100/149‑1 (MQ=255)
                                                                          aCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                >  4:305151/1‑148 (MQ=255)
                                                                            cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTAGCGGCCTTGCCCGCGCAGAGGCCGGCACAGacac                                                                                 >  3:348265/1‑113 (MQ=255)
                                                                              cAAGTACGCTGCCGTAGCCGAAGGCATGACGGCTGCCCTCGGGCTGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTGCGCGGCCTGGCCCGCCAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                           <  6:83247/149‑1 (MQ=255)
                                                                              cAAGTAAGCTGAAGTAGCCGAACGCATGAAGGCTGCCGTCGGGAAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCGGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                           <  6:311764/149‑1 (MQ=255)
                                                                                     gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGTGCAGCCTTTTTTATTTGTCCTTTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                    <  5:88568/149‑1 (MQ=255)
                                                                                      gtGCAGTCCCCGACCCCAAGCAGGCGGCCTTGGGGCAGCTTTTTTTTTTTGTCGTGGGCCGG‑‑‑CCTCTGCGCGGCCTGGCCCGCGATGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                   <  4:321459/148‑1 (MQ=255)
                                                                                      cggCAGTAGCCGACCGCATGACGGCTGCCCTCGGGAAGCCTTTTTTATTTGTCCTGTCCCGG‑‑‑CCGCTTCGCGGCCTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                   <  8:271004/147‑1 (MQ=255)
                                                                                        gCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCGTGGCCGCGAAGAGGGCGGCCCAGACACCACAAAACCTCAGTTTTTATCCGCCCGCGGACGCAGCCCGTTCAgggg                                 >  7:145212/1‑148 (MQ=255)
                                                                                           gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGCCACCACAAAACCTCAGTGTTTGTACGCCCGCTCACGCAGCCCGGTGAGGGTAtt                              >  3:95692/1‑149 (MQ=255)
                                                                                            tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCACTTAGCGGCCTTGGCCGCGAATAGTCCGGGCCAGCCACCACAAAAACGCCCTTTTTATCAGCCCCCTCCCGCCCCCCATTCAGGGGAtt                              >  4:198856/1‑148 (MQ=255)
                                                                                              gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGa                            <  3:305151/148‑1 (MQ=255)
                                                                                              gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCGTGCCCGGGAAGAGGGCGGGGCCGACCCCACAAAACCTCAGGTTTTCTCCGCCCGCTCCCGCAGCCCGTTGAGGGTATTgta                           >  7:137000/1‑147 (MQ=255)
                                                                                               cccaaCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                          <  7:57438/146‑1 (MQ=255)
                                                                                               ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                          >  7:283798/1‑149 (MQ=255)
                                                                                                 gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACCCCACAAAACCTCAGTTTTTCTCCGCCCGCTGACGCAGCCCGTTCAGGGGTTTGAACgg                        >  1:154106/1‑149 (MQ=255)
                                                                                                      aaaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGCGGCCGGCACCGACACCACAAAACCGCAGTTTTTGTCCGCCCGCGCACGCAGCCCTTTCAGGGGATTGAACGGCGCAt                   >  7:28431/5‑149 (MQ=255)
                                                                                                       aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGGCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                              <  7:272231/137‑1 (MQ=255)
                                                                                                        aaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                  >  2:82566/3‑148 (MQ=255)
                                                                                                             gCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                                       >  3:34679/1‑106 (MQ=38)
                                                                                                                 cccgagggCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc        <  1:204343/144‑1 (MQ=255)
                                                                                                                  ccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGGCCGGGAAGGGGCCGGCACAGACACCACAAAACCTGAGTTTTTCTGCGCCCGCTCACGCCGCCCGTTCAGGGTATTGAACGGCGCATCCACCCCGAACt       >  8:212879/1‑149 (MQ=255)
                                                                                                                      gggCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt   <  5:272202/149‑1 (MQ=255)
                                                                                                                      gggCAGCCTTTTTTATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt   <  8:283798/149‑1 (MQ=255)
                                                                                                                        gtcagcattttttATTTGGCCGGGGCCGC‑‑‑CCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg  <  4:40131/142‑1 (MQ=255)
                                                                                                                                  tttATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt   >  5:357219/1‑138 (MQ=255)
                                                                                                                                  tttATTTGTCCTGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt   <  6:357219/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGG‑‑‑CCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740679‑4740947

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