Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,860,706 C→T 8.3% M77I (ATG→ATA) ‡ purN ← phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,860,7060CT 8.3% 169.6 / 3.4 60M77I (ATG→ATA) ‡purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (3/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGC  >  NC_002947/1860560‑1860853
                                                                                                                                                   |                                                                                                                                                    
tGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCAt                                                                                                                                                    <  7:238743/149‑1 (MQ=255)
   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                                                                                                                                            <  5:45909/95‑1 (MQ=255)
                                                             gATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                                            >  6:45909/1‑95 (MQ=255)
                                                              aTAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                       <  8:23111/147‑1 (MQ=255)
                                                                 gTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                      <  4:94580/129‑1 (MQ=255)
                                                                 gTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                      >  3:94580/1‑129 (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)
                                                                                                                                                   |                                                                                                                                                    
TGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGC  >  NC_002947/1860560‑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: ‑