Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 407,765 A→T 7.1% A80A (GCT→GCA aceF ← AceF‑S‑acetyldihydrolipoate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947407,7650AT 7.1% 245.2 / 7.3 85A80A (GCT→GCAaceFAceF‑S‑acetyldihydrolipoate
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base T (3/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

ACGAACCGATGTCCGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGTCAGCAGGCTCTG  >  NC_002947/407623‑407911
                                                                                                                                                |                                                                                                                                                  
acgaACCGATGTCCGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAgcggcg                                                                                                                                              >  3:215384/1‑149 (MQ=255)
     ccGATGTCCGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTc                                                                                                                                         <  1:310301/149‑1 (MQ=255)
       gctGTCCGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCGGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGCGCCGCTGCTGCGGcccctgag                                                                                                                                        <  8:138124/147‑5 (MQ=255)
       gATGTCCGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGCGCCGCTGCTGCGGcccctgag                                                                                                                                        <  3:256114/149‑5 (MQ=255)
           tCCGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTc                                                                                                                                   >  4:264924/1‑149 (MQ=255)
            ccGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTc                                                                                                                                   >  6:5213/1‑148 (MQ=255)
              ggCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCtga                                                                                                                                        >  7:146769/1‑139 (MQ=255)
              ggCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCtga                                                                                                                                        <  8:146769/141‑3 (MQ=255)
              ggCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAg                                                                                                                                >  5:157689/1‑149 (MQ=255)
                   gTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCa                                                                                                                                 <  2:253243/143‑1 (MQ=255)
                   gTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCATCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCa                                                                                                                                 >  1:253243/1‑143 (MQ=255)
                      gATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACcagca                                                                                                                        <  6:189302/149‑1 (MQ=255)
                      gATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACcagc                                                                                                                         <  7:9020/148‑1 (MQ=255)
                        tGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGt                                                                                                                      <  7:13470/149‑1 (MQ=255)
                           ccTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTtcg                                                                                                                   >  3:306985/1‑149 (MQ=255)
                             tGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGCGCCGCTGCTGCGGcccctgag                                                                                                                                        <  7:294116/127‑5 (MQ=255)
                             tGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGCGCCGCTGCTGCGGcccctgag                                                                                                                                        >  8:294116/1‑123 (MQ=255)
                              tggagCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTtcgtcg                                                                                                                >  4:311697/5‑149 (MQ=255)
                               cACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTtcgtc                                                                                                                 <  6:157689/147‑1 (MQ=255)
                               cACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTtcgtc                                                                                                                 <  6:289519/147‑1 (MQ=255)
                                  gCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCAACCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCATCGGCTTCCAGGACCAGCAGTTCGGCGCct                                                                                                             >  1:92940/1‑149 (MQ=255)
                                   ccggCCGCAGCCCCGGCTGCTGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCcttc                                                                                                           >  5:84569/3‑149 (MQ=255)
                                   ccggCCGCAGCCCCGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCcttc                                                                                                           >  7:213175/3‑149 (MQ=255)
                                     ggCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCcttctt                                                                                                         >  5:37862/1‑149 (MQ=255)
                                     ggCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCcttctt                                                                                                         >  7:22929/1‑149 (MQ=255)
                                      gccgcagcCCCGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCttt                                                                                                        <  8:213175/146‑1 (MQ=255)
                                       cagcagCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCttt                                                                                                        >  5:135333/1‑148 (MQ=255)
                                       cagcagCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTc                                                                                                       <  2:289700/149‑1 (MQ=255)
                                          catccGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACAAGCAGTTCGTCGCCTTCTTTCAg                                                                                                     <  7:203482/145‑1 (MQ=255)
                                           agcCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGc                                                                                                   >  6:33922/1‑149 (MQ=255)
                                            gcCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCg                                                                                                  <  8:296047/149‑1 (MQ=255)
                                               ggggCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTc                                                                                               <  5:5213/149‑1 (MQ=255)
                                               ccggCTGCTGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTc                                                                                               >  7:164292/3‑149 (MQ=255)
                                                gggCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCg                                                                                              <  3:264924/149‑1 (MQ=255)
                                                gggCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTCCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCg                                                                                              <  4:249574/149‑1 (MQ=255)
                                                         ccGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTg                                                                                     <  3:59725/149‑1 (MQ=255)
                                                         ccGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTg                                                                                     <  8:22929/149‑1 (MQ=255)
                                                          cGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGa                                                                                    >  6:12935/1‑149 (MQ=255)
                                                           ggAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGAc                                                                                   >  6:291352/1‑149 (MQ=255)
                                                           ggAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGAc                                                                                   >  4:7925/1‑149 (MQ=255)
                                                           ggAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGAc                                                                                   >  4:62618/1‑149 (MQ=255)
                                                            gAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAg                                                                                                     <  6:329896/130‑1 (MQ=255)
                                                            gAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAg                                                                                                     >  5:329896/1‑130 (MQ=255)
                                                                 gcagGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAg                                                                             >  7:72329/1‑149 (MQ=255)
                                                                    ggAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTc                                                                          <  4:306985/149‑1 (MQ=255)
                                                                     gAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTcc                                                                         >  2:70588/1‑149 (MQ=255)
                                                                     gAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTcc                                                                         <  3:7925/149‑1 (MQ=255)
                                                                         ggCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGa                                                                     <  3:311697/149‑1 (MQ=255)
                                                                          gCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCACTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACATTCAGTTCCTTGa                                                                     >  6:56105/1‑148 (MQ=255)
                                                                           cTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCtt                                                                       <  6:56464/145‑1 (MQ=255)
                                                                           cTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCtt                                                                       >  5:56464/1‑145 (MQ=255)
                                                                           cTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGAt                                                                    >  8:34029/1‑148 (MQ=255)
                                                                             tCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGTAGCAGCGGCGCCCTCGGCTTCCa                                                                                                                                 <  4:21790/85‑1 (MQ=255)
                                                                             tCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGTAGCAGCGGCGCCCTCGGCTTCCa                                                                                                                                 >  3:21790/1‑85 (MQ=255)
                                                                              cAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACg                                                                >  3:207762/1‑149 (MQ=255)
                                                                                    ctggcgctggcgCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGcc                                                          <  7:34029/149‑1 (MQ=255)
                                                                                     tggcgctggcgcCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCa                                                                                                      <  4:238750/104‑1 (MQ=255)
                                                                                     tggcgctggcgcCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCa                                                                                                      >  3:238750/1‑104 (MQ=255)
                                                                                       gcgctggcgcCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCtt                                                        <  1:70588/148‑1 (MQ=255)
                                                                                       gcgctggcgcCGCAGCCGGAGCAGCGGCAGCAGCCGACGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCtt                                                        <  5:321425/148‑1 (MQ=255)
                                                                                        cgctggcgcCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTgg                                                      <  8:164292/149‑1 (MQ=255)
                                                                                             gcgcCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGcc                                                 <  5:12935/149‑1 (MQ=255)
                                                                                             gcgcCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGcc                                                 <  6:84569/149‑1 (MQ=255)
                                                                                                   gcAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTggg                                                     >  3:196081/1‑140 (MQ=255)
                                                                                                   gcAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTggg                                                     <  4:196081/140‑1 (MQ=255)
                                                                                                     aGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCTAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAAttt                                          >  5:87256/1‑149 (MQ=255)
                                                                                                         ggAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCACCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGGCTTCTTTCAGGCGGTCGACCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCAt                                      >  1:173668/1‑149 (MQ=255)
                                                                                                          gAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATg                                     <  4:166730/149‑1 (MQ=255)
                                                                                                            gcagcGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGa                                   <  8:72329/149‑1 (MQ=255)
                                                                                                             cagcGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCGAGCCTTGGGCGCCGGAATTTCCATGGAg                                  <  1:185688/149‑1 (MQ=255)
                                                                                                                cGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGcc                               >  6:76468/1‑149 (MQ=255)
                                                                                                                 ggCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCt                              >  1:128857/1‑149 (MQ=255)
                                                                                                                    agcagcCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCGGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGt                           >  5:127612/1‑149 (MQ=255)
                                                                                                                        gcCGGCGCCTCAGCGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCggt                       <  6:161502/149‑2 (MQ=255)
                                                                                                                             cgcCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTc                    >  3:279523/1‑147 (MQ=255)
                                                                                                                               ccTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAgg                <  3:149916/149‑1 (MQ=255)
                                                                                                                               ccTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAgg                <  4:215384/149‑1 (MQ=255)
                                                                                                                                cTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTc                          <  4:86735/138‑1 (MQ=255)
                                                                                                                                cTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTc                          >  3:86735/1‑138 (MQ=255)
                                                                                                                                 tCAGTGGCCGCAGCAGCGGCGCCATCGGCTTCCATGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTACTTGATCACGCAAGCCTTGGGCGCCGGAATTTCCATGGAGTACTTGTCGGACTCAATGGt              >  6:248402/1‑149 (MQ=255)
                                                                                                                                 cgagggcgCCGCTGCTGCGGCCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGACGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGt              >  5:28759/8‑149 (MQ=39)
                                                                                                                                   aGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGt              >  4:302660/1‑147 (MQ=255)
                                                                                                                                   aGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGt              >  7:224940/1‑147 (MQ=255)
                                                                                                                                    gggcgCCGCTGCTGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGt              >  8:191088/5‑147 (MQ=255)
                                                                                                                                         ccGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGTCAGCAgg       <  2:92940/149‑1 (MQ=255)
                                                                                                                                           gcagcaGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGTCAGCAGGct     <  6:37862/149‑1 (MQ=255)
                                                                                                                                           gcagcaGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGTCAGCAGGc      >  6:31995/1‑148 (MQ=255)
                                                                                                                                              gcagcGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGTCAGCAGGCTCTg  <  1:130461/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ACGAACCGATGTCCGGCACGTGGATGTCCTGCACGCTGGCAGCAGCCGGGGCTGCGGCCGGAGCAGCAGGAGCGGCTTCAGCGGCTGGCGCTGGCGCCGCAGCCGGAGCAGCGGCAGCAGCCGGCGCCTCAGGGGCCGCAGCAGCGGCGCCCTCGGCTTCCAGGACCAGCAGTTCGTCGCCTTCTTTCAGGCGGTCGCCCAGCTTGACCTTCAGTTCCTTGATCACGCCAGCCTTGGGCGCCGGAATTTCCATGGAGGCCTTGTCGGACTCCAGGGTCAGCAGGCTCTG  >  NC_002947/407623‑407911

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: