Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,845,218 T→G 8.4% P43P (CCT→CCG)  coaD → phosphopantetheine adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,845,2180TG 8.4% 176.3 / 4.1 57P43P (CCT→CCG) coaDphosphopantetheine adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base G (2/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCAGCCAGGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGG  >  NC_002947/5845070‑5845363
                                                                                                                                                    |                                                                                                                                                 
aCTCCAGCCAGGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCt                                                                                                                                                   >  3:303677/1‑149 (MQ=255)
aCTCCAGCCAGGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCt                                                                                                                                                   >  2:24456/1‑149 (MQ=255)
  tcaagccagCATGACGCGATGAAACGACTGTTGTACCCGGGAACTTTCTACCCCATTACCAAAGGCAATGGCGATCTGGTCGAGAGCGCTTACCTGCTATTCGACCACTTGATCATCGCCTTGGCGGCCAGCCCCAAGAAAAACCCatt                                                                                                                                                 >  4:108804/4‑149 (MQ=255)
  tCCAGCCAGGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCttt                                                                                                                                                 >  2:277392/1‑149 (MQ=255)
  tCCAGCCAGGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCtt                                                                                                                                                  >  6:24146/1‑148 (MQ=255)
    cagccagGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCtt                                                                                                                                                  <  5:24146/146‑1 (MQ=255)
    cagccagGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGt                                                                                                                                               >  4:131679/1‑149 (MQ=255)
       ccagGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTcc                                                                                                                                            >  1:290761/1‑149 (MQ=255)
       ccagGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTcc                                                                                                                                            >  1:46617/1‑149 (MQ=255)
            aTGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTg                                                                                                                                       <  4:174122/149‑1 (MQ=255)
              gACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCGTTGTTccaggggg                                                                                                                                      >  5:302581/1‑142 (MQ=255)
                cgcgATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGaa                                                                                                                                    >  5:18589/1‑148 (MQ=255)
                 gcgATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGaac                                                                                                                                   >  3:144746/1‑148 (MQ=255)
                       aaCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGGGGAACAAAggg                                                                                                                             <  6:302581/148‑1 (MQ=255)
                        aCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGtt                                                                                                                           <  6:80701/149‑1 (MQ=255)
                        aCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGtt                                                                                                                           >  4:1928/1‑149 (MQ=255)
                           gAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAg                                                                                                                        <  3:131679/149‑1 (MQ=255)
                              tgttgtACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGTTCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGGGGAACAAAGGGtttttctt                                                                                                                     <  6:225271/149‑7 (MQ=38)
                                 tgtACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGc                                                                                                                   <  3:108804/148‑1 (MQ=255)
                                   tACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCg                                                                                                                >  5:36559/1‑149 (MQ=255)
                                     cccGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGt                                                                                                               <  2:46617/148‑1 (MQ=255)
                                           tACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTg                                                                                                                                                <  6:266817/109‑1 (MQ=255)
                                           tACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTg                                                                                                                                                >  5:266817/1‑109 (MQ=255)
                                             cTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCAcc                                                                                                      <  1:277392/149‑1 (MQ=255)
                                             cTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCAcc                                                                                                      <  7:40859/149‑1 (MQ=255)
                                              tttCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAg                                                                                                            <  2:281609/142‑1 (MQ=255)
                                              tttCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAg                                                                                                            >  1:281609/1‑142 (MQ=255)
                                              tttCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCAcc                                                                                                      >  4:199911/1‑148 (MQ=255)
                                               ttCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGtt                                                                                                                           >  7:75087/1‑126 (MQ=255)
                                               ttCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGtt                                                                                                                           <  8:75087/126‑1 (MQ=255)
                                                 cGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCACGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAg                                                                                                   <  6:210724/148‑1 (MQ=255)
                                                     cccATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAGCGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCt                                                                                              <  1:299186/149‑1 (MQ=255)
                                                     cccATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTAATCGACCACGTGATCCGCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGACCTCCGGGTTGAGCTCGCCCGTGAGGTCAGCAAGCACCt                                                                                              >  7:214709/1‑149 (MQ=255)
                                                      ccATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAg                                                                                                            >  4:16797/1‑134 (MQ=255)
                                                      ccATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAg                                                                                                            <  3:16797/134‑1 (MQ=255)
                                                       cATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGc                                                                                            <  6:18589/149‑1 (MQ=255)
                                                        aTTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGcc                                                                                           >  1:244150/1‑149 (MQ=255)
                                                           aCCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCa                                                                                         <  5:60673/148‑1 (MQ=255)
                                                             cAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATg                                                                                      >  5:108172/1‑149 (MQ=255)
                                                                  gCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGaa                                                                                 <  1:24456/149‑1 (MQ=255)
                                                                   ccATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGGGGAACAAAGGGtt                                                                                                                           >  4:154396/1‑106 (MQ=38)
                                                                   ccATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGGGGAACAAAGGGtt                                                                                                                           <  3:154396/106‑1 (MQ=38)
                                                                        gCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATc                                                                           <  3:1928/149‑1 (MQ=255)
                                                                        gCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATc                                                                           <  2:244150/149‑1 (MQ=255)
                                                                                gTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGc                                                                        <  3:199911/144‑1 (MQ=255)
                                                                                gTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGc                                                                        <  2:290761/144‑1 (MQ=255)
                                                                                                     ttCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGc                                                                                                                   <  3:114419/80‑1 (MQ=255)
                                                                                                     ttCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGc                                                                                                                   >  4:114419/1‑80 (MQ=255)
                                                                                                     ttCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTc                                                <  4:144746/147‑1 (MQ=255)
                                                                                                     ttCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCg                                               >  8:48786/1‑148 (MQ=255)
                                                                                                         aCCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAg                                          <  6:36559/149‑1 (MQ=255)
                                                                                                                      cgcgGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCAAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGcc                              >  4:243215/1‑148 (MQ=255)
                                                                                                                               ggcCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTcc                    <  6:108172/149‑1 (MQ=255)
                                                                                                                                  cAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGc                 >  4:46335/1‑149 (MQ=255)
                                                                                                                                  cAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGc                 >  8:3394/1‑149 (MQ=255)
                                                                                                                                  cAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGc                 <  7:48786/149‑1 (MQ=255)
                                                                                                                                     ccccAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTg              >  4:95817/1‑149 (MQ=255)
                                                                                                                                      cccAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTgg             >  6:241468/1‑149 (MQ=255)
                                                                                                                                       ccAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGc            >  1:210281/1‑149 (MQ=255)
                                                                                                                                       ccAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTCAGCTCGCCCGTCAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGc            >  3:98719/1‑149 (MQ=255)
                                                                                                                                        cAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGcc           >  4:239078/1‑149 (MQ=255)
                                                                                                                                              aaaCCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGACAAGGAACAGGGCGCCGACGTCTTCCTGCGTGGCCtgcggg     >  3:238928/1‑147 (MQ=255)
                                                                                                                                                 cccTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCgg  >  7:197407/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACTCCAGCCAGGATGACGCGATGAACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGG  >  NC_002947/5845070‑5845363

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