Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,860,710 A→G 8.4% L76P (CTA→CCA)  purN ← phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,860,7100AG 8.4% 160.7 / 4.8 59L76P (CTA→CCA) purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base G (3/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGC  >  NC_002947/1860563‑1860853
                                                                                                                                                     |                                                                                                                                               
gTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑aa                                                                                                                                                <  8:287800/149‑2 (MQ=255)
    gCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                                            >  6:56874/1‑149 (MQ=255)
     cAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                                            <  5:224240/148‑1 (MQ=255)
         cccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCCGCGCTAAT‑GGcgcgcatc                                                                                                                                        <  1:64689/149‑5 (MQ=255)
          ccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCg                                                                                                                                      <  8:212532/149‑1 (MQ=255)
           cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATc                                                                                                                                       <  8:136835/147‑1 (MQ=255)
           cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGa                                                                                                                                     >  7:136835/1‑149 (MQ=255)
              gTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGaa                                                                                                                                    <  8:83323/147‑1 (MQ=255)
                   tGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCa                                                                                                                             >  8:129761/1‑149 (MQ=255)
                      gcagcaGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACgg                                                                                                                          <  4:278008/149‑1 (MQ=255)
                          cagcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGt                                                                                                                      <  1:82154/149‑1 (MQ=255)
                             cGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCg                                                                                                                    >  1:175118/1‑148 (MQ=255)
                               aCGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAAc                                                                                                                 <  6:276291/149‑1 (MQ=255)
                                  ggTGCAGGTTCAGCAGGCGGCCCTGATAGTGGCCCACGAAGCCGCCACTGAGGAGGCGCATCAAACCGCCCAGCACCACCAGGTCCGGGGAAAAGCCGTCTATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                              <  2:175118/149‑1 (MQ=255)
                                   gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCttt                                                                                                                              <  5:64671/132‑2 (MQ=37)
                                   gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCttt                                                                                                                              >  6:64671/1‑131 (MQ=37)
                                       aTGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                                         >  6:339426/1‑149 (MQ=255)
                                                   cggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGa                                                                                                    <  5:8035/142‑1 (MQ=255)
                                                   cggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGa                                                                                                    >  6:8035/1‑142 (MQ=255)
                                                          gATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                                            >  6:45909/1‑95 (MQ=255)
                                                          gATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                                            <  5:45909/95‑1 (MQ=255)
                                                           aTAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                       <  8:23111/147‑1 (MQ=255)
                                                              gTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                      >  3:94580/1‑129 (MQ=255)
                                                              gTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                      <  4:94580/129‑1 (MQ=255)
                                                                       gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTg                                                                         >  4:258618/1‑149 (MQ=255)
                                                                       gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTg                                                                         >  4:297269/1‑149 (MQ=255)
                                                                            cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                              >  8:227977/1‑107 (MQ=255)
                                                                            cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                              <  7:227977/107‑1 (MQ=255)
                                                                            cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgcggcg                                                                    >  4:25371/1‑149 (MQ=255)
                                                                              ccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGcgcg                                                                  <  6:148733/149‑1 (MQ=255)
                                                                                 cTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTg                                                               >  7:63592/1‑149 (MQ=255)
                                                                                    aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                             <  7:129761/148‑1 (MQ=255)
                                                                                      gATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGc                                                          >  3:236325/1‑149 (MQ=255)
                                                                                        tGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgtgt                                                                                                          <  4:268229/99‑1 (MQ=255)
                                                                                        tGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgtgt                                                                                                          >  3:268229/1‑99 (MQ=255)
                                                                                           gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTGAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCg                                                        >  1:277678/1‑147 (MQ=255)
                                                                                           gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                     <  5:56874/149‑1 (MQ=255)
                                                                                           gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGCTGTGGTCGAGCACGACCCTGTCGATTCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                     >  3:248594/1‑149 (MQ=255)
                                                                                           gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATAGAAGGCTTCAAGGCCGTCGAACTGGGTGTGGTCGAGCACCACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGGCGGAg                                                     >  2:173536/1‑149 (MQ=255)
                                                                                               gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGAGCGGCATCGAACGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                 >  5:61848/1‑149 (MQ=255)
                                                                                                 aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCATGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                               >  8:37886/1‑149 (MQ=255)
                                                                                                 aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCATGCCGTAGGCATCg                                               >  4:310290/1‑149 (MQ=255)
                                                                                                 aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                               <  4:236325/149‑1 (MQ=255)
                                                                                                         aGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGtt                                       <  5:339426/149‑1 (MQ=255)
                                                                                                          gCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                 >  6:148792/1‑138 (MQ=255)
                                                                                                          gCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                 <  5:148792/138‑1 (MQ=255)
                                                                                                           caccacCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTg                                      >  7:156848/1‑148 (MQ=255)
                                                                                                             ccaccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCCTT‑AGCGCGGGATGGAAGGCTTGACGGGCGGCGAACGGGGGGGGGTCGAGGACGACGCTGTCGCTGCCGGCCGCTGCGGGGGGGTGGAGGCCGGAGGCATCGGCGCGGTTggc                                    >  2:170975/1‑147 (MQ=255)
                                                                                                              caccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGa                                    >  7:99926/1‑147 (MQ=255)
                                                                                                               accaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCAAT‑AGCGCGCCATCCAAGGCTAAACGGCCGGCCAACAGCGAGTGGACGAGCACGAAGCACTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAAAGGCCCGCTTGGAAAg                                 >  2:124088/1‑149 (MQ=255)
                                                                                                                    gTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcg                            <  1:173536/149‑1 (MQ=255)
                                                                                                                         gggCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCgg                       <  8:156848/149‑1 (MQ=255)
                                                                                                                         gggCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCgg                       <  8:99926/149‑1 (MQ=255)
                                                                                                                           gCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGAt                     <  3:258618/149‑1 (MQ=255)
                                                                                                                             aaaGCCGTCGATG‑CGCGCCATA‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATAc                   <  1:170975/149‑1 (MQ=255)
                                                                                                                               aGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGc                 <  4:248594/149‑1 (MQ=255)
                                                                                                                               aGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGc                 <  3:25371/149‑1 (MQ=255)
                                                                                                                                gCCTTCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCa                >  3:22125/1‑149 (MQ=38)
                                                                                                                                 ccttCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc               >  7:255816/4‑149 (MQ=38)
                                                                                                                                               ccATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGAGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGc  <  5:258438/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGC  >  NC_002947/1860563‑1860853

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: ‑