Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,819 2 bp→CT intergenic (+56/+696) 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,8190TC94.6% 165.2 / ‑6.0 59intergenic (+56/+697)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  major base C (28/25);  minor base G (0/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029474,740,8200GT98.3% 190.0 / NA 59intergenic (+57/+696)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (29/29);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGG  >  NC_002947/4740676‑4740961
                                                                                                                                                ||                                                                                                                                              
ctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTATTTATTTGTCCTGTGcc                                                                                                                                           >  5:155014/1‑149 (MQ=255)
   cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAACTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGc                                                                                                                                       <  2:74614/149‑1 (MQ=255)
        cGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                                   >  4:148171/1‑148 (MQ=255)
           aCAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCAGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCATGTGCCGGCCTCTTcg                                                                                                                                >  1:233110/1‑149 (MQ=255)
               gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                                           >  7:145291/1‑149 (MQ=255)
               gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                                           >  3:218687/1‑149 (MQ=255)
                aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                                          <  2:6562/149‑1 (MQ=255)
                    gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGAGACATAGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAATGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGcc                                                                                                                      <  6:185617/149‑1 (MQ=255)
                       gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                                   <  2:4993/149‑1 (MQ=255)
                        gacgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                                   >  3:66453/3‑148 (MQ=255)
                          cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCCAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                <  7:159321/149‑1 (MQ=255)
                           gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                               <  5:104820/149‑1 (MQ=255)
                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                              <  2:267443/149‑1 (MQ=255)
                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                              <  2:296531/149‑1 (MQ=255)
                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                              <  2:8023/149‑1 (MQ=255)
                             gcCAGCTGTACACCGGCGAGAAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                             <  6:219235/149‑1 (MQ=38)
                                aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                          <  8:145291/149‑1 (MQ=255)
                                  cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                                        <  6:71071/149‑1 (MQ=255)
                                   tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                                        >  2:215942/1‑148 (MQ=255)
                                      acacCGGCGTGCAGAAGCGCGACATCGTGGCCCTGAAACACCCCAAGTAAGCTGCAGTAGCCCAACGCAACATGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                                    <  2:264615/149‑1 (MQ=255)
                                               gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                              >  6:44063/1‑131 (MQ=38)
                                               gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                              <  5:44063/131‑1 (MQ=38)
                                                aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcaccac                                                                                            >  1:254801/1‑148 (MQ=255)
                                                   aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCGGGAGAGGCCGGCACCGACACcacatct                                                                                        >  8:311102/1‑146 (MQ=255)
                                                      aGCGCCACATCGTTGCCCTGAACGACCGCACTTACGCTCCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAAcc                                                                                      <  4:218687/148‑1 (MQ=255)
                                                       gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAcc                                                                                      >  7:314920/1‑147 (MQ=255)
                                                           gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                                                <  2:254801/149‑1 (MQ=255)
                                                            aCATCGTTCCCCTGAAAGACCGCATGTAAGCTGCAGTAGCCGAACTCAAGAAGGATGCCCTCGGGCAGCTTTTTTTATTTGTCCTGTCCCGGCCTCATGGGGGCCTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGttt                                                                               <  3:137931/149‑1 (MQ=255)
                                                             cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtttt                                                                              >  8:200509/1‑149 (MQ=255)
                                                               tccttGCCATTAAAGACCGCAATTAAGCTGCATTAGCCGAACGAAAGAAGGCTGCCCTCGGCCAGCCTTTTTTATTTGTCCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                            <  4:72481/146‑1 (MQ=255)
                                                               tCGTTGCCCTGAATCACCGCATGTATGCTGCATTAGCCTATCTCAAGAAGGGTGCCCTCGGGCAGCTTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                            <  6:155014/149‑1 (MQ=255)
                                                                  ttGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGGGAAGAGGCCGGCACAGACACCACAAAACCTGAGTTTTTGTcc                                                                         >  7:235667/1‑149 (MQ=255)
                                                                         gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTCTCCGCCCGCt                                                                  >  3:108676/1‑149 (MQ=255)
                                                                             gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                              >  6:63902/1‑149 (MQ=255)
                                                                                         gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                  <  1:215942/149‑1 (MQ=255)
                                                                                          ctgctgTACCCTAACGCATGACGCCGGCCTTCGGGCAGAGTTTTTTATTTGTCGTGTGCGGGCATCTTCGCGGCCTTGCCCGGTATTAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                 <  8:56323/149‑1 (MQ=255)
                                                                                               gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                            <  7:259944/149‑1 (MQ=255)
                                                                                               gTAGCCGAACGCAACTAGGCGGCCCTCGGGCAGCCTTTTTTTTTTGTCCTGGGCCGGCCTCTGCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                            <  8:314920/149‑1 (MQ=255)
                                                                                                   ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGa                                          >  4:251222/1‑147 (MQ=255)
                                                                                                        cGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcg                                    >  6:337611/1‑148 (MQ=255)
                                                                                                           aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTAGCGGCCTTGGCCGCGAAGCGGCCGGCCCAGCCACCACAAAACCTCAGTTTTTCTCCGCCCCCGCACGCAGCCCTTTGAGGGGATTGAACGGCGCCTc                                >  6:207398/1‑149 (MQ=255)
                                                                                                            agacgccTGCCCTCGGCAAGCCTTTTTAATTTGACCTGGGCCGCCCTCTGCGCGCCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc                               <  7:235503/143‑1 (MQ=255)
                                                                                                             gaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATACGCCCGCTCACGCAGACCTTTAAGGGTATTGAACGGCGCATcca                              >  8:75301/1‑149 (MQ=255)
                                                                                                              aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCACCACAAAACCTCAGTTTTTCTACGCCCGCACACGCAGCCCTTTCAGGGTATTGAACGGCCCATcccc                             >  6:294400/1‑147 (MQ=255)
                                                                                                                ggCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCCGTTTTTGTGCGCCCGCTCACGCCGCCCGTTCAGGGTATTGAACGGCGCAAcgaccg                           >  8:6821/1‑148 (MQ=255)
                                                                                                                  cTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTATCCGCCCGCTCACGCAGCCCTTTAAGGGTATTGAAAGGCGCATccac                             >  4:281308/1‑145 (MQ=255)
                                                                                                                  cTGCCCGCGGGCAGCCTTTTTTATTTGTCCTGGGCCGCCCTCTTCGCGCCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                             <  3:281308/145‑1 (MQ=255)
                                                                                                                    gCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGACCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTGAGGGTATTGAACGGCGCATccacagctc                        >  1:238461/1‑144 (MQ=255)
                                                                                                                    gCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                              <  2:238461/142‑1 (MQ=255)
                                                                                                                     cccTCGGGCAGCTTTTTTTATTTGGCGTGGGCCGCCCGCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                      <  3:68704/149‑1 (MQ=255)
                                                                                                                       cTCGGGCAGCTTTTTTTATTTGTCCTGGGCCGGCCTCTTCGGGGCGTTGCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAAGGTATTGAACGGCGCATCCACCACGAACtt                    <  4:168550/149‑1 (MQ=255)
                                                                                                                              cagcgttttttATTTGGCGTGGGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc              <  8:235667/144‑1 (MQ=255)
                                                                                                                                   ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc          >  2:234375/1‑148 (MQ=255)
                                                                                                                                   ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc          >  8:259745/1‑148 (MQ=255)
                                                                                                                                   ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc          >  2:41423/1‑148 (MQ=255)
                                                                                                                                   ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCCGCACAGAAACCACAAAACCACACTTTTTATACGCCCGCTCACGCAGCCCTTTAAAGGTATTTAACGGAGCATACACCACGAACTTGTTGGccacc          >  6:180275/1‑148 (MQ=255)
                                                                                                                                   ttttttATTTGTCCTGTGCCGGCCACTTAGCGGCCATTCCCGCGAAGAGGCCGGCACAGCCACCACAAAACCGCAGTTTTTATCCCCCCCCGCACCCAGCCCATTAAGGGTATTTAACAGCGGATCCCACAAGAAATTGTTTGccaac          >  4:194417/1‑146 (MQ=255)
                                                                                                                                       ttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTGAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGc     >  7:89180/1‑149 (MQ=255)
                                                                                                                                          tttgtCTTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg  <  1:113439/149‑1 (MQ=255)
                                                                                                                                                ||                                                                                                                                              
CTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGG  >  NC_002947/4740676‑4740961

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