Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,473,817 G→A 8.3% D326D (GAC→GAT proA ← gamma‑glutamyl phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,473,8170GA 8.3% 164.9 / 4.6 60D326D (GAC→GATproAgamma‑glutamyl phosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (25/30);  new base A (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-02

AAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTTCTGGCA  >  NC_002947/5473674‑5473965
                                                                                                                                               |                                                                                                                                                    
aaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGAc                                                                                                                                                 <  2:274103/149‑1 (MQ=255)
aaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGACGGCCTGAACCAGGCCATCGag                                                                                                                                                 <  5:170608/149‑2 (MQ=255)
          gggTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGGCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATc                                                                                                                                       >  1:263894/1‑149 (MQ=255)
               ttGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCGTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGGCGACCACGCGAATCGACa                                                                                                                                   <  1:282362/148‑1 (MQ=255)
               ttGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGt                                                                                                                                                     >  7:289858/1‑130 (MQ=255)
               ttGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGt                                                                                                                                                     <  8:289858/130‑1 (MQ=255)
                       gACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCAAAAGGATCGCAt                                                                                                                          >  1:271250/1‑149 (MQ=255)
                        aCCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATc                                                                                                                         >  1:147846/1‑149 (MQ=255)
                         ccGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATc                                                                                                                         <  2:263894/148‑1 (MQ=255)
                          cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATc                                                                                                                         >  1:121336/1‑147 (MQ=255)
                             cGCCTAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGt                                                                                                                    <  5:14293/149‑1 (MQ=255)
                             cGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGt                                                                                                                    <  6:19837/149‑1 (MQ=255)
                                    tCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCgt                                                                                                             <  2:121336/149‑1 (MQ=255)
                                      gACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCgtg                                                                                                            <  4:305457/148‑1 (MQ=255)
                                       aCCTCAGCCATGAACTGACGGGCTTCACCCTTGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                                      <  3:92624/137‑1 (MQ=255)
                                       aCCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                                      >  4:92624/1‑137 (MQ=255)
                                        ccTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGc                                                                                                         <  6:275509/149‑1 (MQ=255)
                                        ccTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGc                                                                                                         <  6:148470/149‑1 (MQ=255)
                                         ctctGCCATGAACTGACGGGCTTCACCGTGGGGTTCGCTGAGGTTCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCTGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTTTTCCGTGTGcc                                                                                                        <  2:312174/149‑1 (MQ=255)
                                                    aattaggggCTTCTCTCTGGTGTTCGCTGATGATCTAGTCTGTGTGGTTCGAGCCATTGTGGTTGATGTGCTCGATTGCGTGGTTCAGGCCGTCGACCACGGGAATCGACAGGATCGGATCGAGGTATTCCGTGTGCCAGTCGGCTTCg                                                                                             <  7:99161/143‑1 (MQ=255)
                                                    aCTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATATTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGAc                                                                                                                                    <  6:91617/110‑1 (MQ=255)
                                                    aCTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAATGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGAc                                                                                                                                    >  5:91617/1‑110 (MQ=255)
                                                       gACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTg                                                                                          >  5:275611/1‑149 (MQ=255)
                                                          gggCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGt                                                                                           >  5:50983/1‑145 (MQ=255)
                                                          gggCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGt                                                                                           <  6:50983/145‑1 (MQ=255)
                                                            gCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCgg                                                                                     <  6:275611/149‑1 (MQ=255)
                                                              ttCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGtt                                                                                   >  5:215469/1‑149 (MQ=255)
                                                                   ccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGgcg                                                                              >  7:50077/1‑149 (MQ=255)
                                                                     tGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCt                                                                            >  8:39034/1‑149 (MQ=255)
                                                                       gtgtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtga                                                                          >  6:252049/1‑149 (MQ=255)
                                                                       gtgtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtg                                                                           <  7:39034/148‑1 (MQ=255)
                                                                            cGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATa                                                                     >  4:288054/1‑149 (MQ=255)
                                                                             gCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAg                                                                    >  6:44652/1‑149 (MQ=255)
                                                                                   gatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTggg                                                              >  7:116327/1‑149 (MQ=255)
                                                                                          tCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgat                                                                         <  4:101176/131‑1 (MQ=255)
                                                                                          tCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgat                                                                         >  3:101176/1‑131 (MQ=255)
                                                                                          tCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTc                                                       <  5:252049/149‑1 (MQ=255)
                                                                                            ggTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATAGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCg                                                      >  1:53275/1‑148 (MQ=255)
                                                                                             gtgtGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGc                                                     >  4:129786/1‑148 (MQ=255)
                                                                                                   tgCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACcgcg                                              <  2:143284/149‑1 (MQ=255)
                                                                                                    gCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCAAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACcgcg                                              >  4:28063/1‑148 (MQ=255)
                                                                                                     cGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCgcgc                                             <  8:116327/148‑1 (MQ=255)
                                                                                                         ccATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCt                                         <  6:215469/148‑1 (MQ=255)
                                                                                                              gTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACg                                    <  2:147846/148‑1 (MQ=255)
                                                                                                                 attcgcgtgGTCGACGGCCTGAACCAGGCCATCGAGCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCt                                >  8:57728/7‑149 (MQ=38)
                                                                                                                  ttcgcgtgGTCGACGGCCTGAACCAGGCCATCGAGCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCtt                               >  5:36154/6‑149 (MQ=38)
                                                                                                                        tggtCGACGGCCTGAACCAGGCCATCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCg                                                      <  2:27991/117‑1 (MQ=255)
                                                                                                                        tggtCGACGGCCTGAACCAGGCCATCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCg                                                      >  1:27991/4‑120 (MQ=255)
                                                                                                                        tgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACTGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCa                          <  2:271250/148‑1 (MQ=255)
                                                                                                                        tgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCAc                         >  4:266277/1‑149 (MQ=255)
                                                                                                                          cTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACg                        >  4:244670/1‑148 (MQ=255)
                                                                                                                          cTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACg                        >  2:139063/1‑148 (MQ=255)
                                                                                                                             gATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgatga                                                                       <  7:222813/98‑1 (MQ=255)
                                                                                                                             gATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgatga                                                                       >  8:222813/1‑98 (MQ=255)
                                                                                                                             gATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAg                     <  5:44652/148‑1 (MQ=255)
                                                                                                                                  ccTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGcggcgg               >  7:54509/1‑149 (MQ=255)
                                                                                                                                           ggCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTTCt      <  3:266277/149‑1 (MQ=255)
                                                                                                                                           ggCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTTCt      <  8:50077/149‑1 (MQ=255)
                                                                                                                                             ccGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTTCTgg    <  3:28063/149‑1 (MQ=255)
                                                                                                                                             ccGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTTCTgg    <  1:139063/149‑1 (MQ=255)
                                                                                                                                               gTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTTCTGGCa  <  8:54509/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTTCTGGCA  >  NC_002947/5473674‑5473965

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: