Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,470,708 A→G 5.2% V161A (GTG→GCG)  PP_3961 ← transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,470,7080AG 5.2% 230.9 / 2.0 77V161A (GTG→GCG) PP_3961transporter
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base G (2/2);  total (35/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 = 7.26e-01

CCGACATGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGA  >  NC_002947/4470562‑4470852
                                                                                                                                                   |                                                                                                                                                   
ccGACATGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCAcc                                                                                                                                                   <  1:228220/149‑1 (MQ=255)
 cGACATGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCAc                                                                                                                                                    <  4:289762/147‑1 (MQ=255)
      tGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATc                                                                                                                                             <  1:231972/149‑1 (MQ=255)
        cTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcgg                                                                                                                                           <  1:295088/149‑1 (MQ=255)
           tCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcggcg                                                                                                                                         <  6:356064/148‑1 (MQ=255)
           tCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcggc                                                                                                                                          <  4:333024/147‑1 (MQ=255)
             ttGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGCTGTACCGCCACGCTCAGCCCCAGCGCACTCGGCGGGc                                                                                                                                      >  8:120531/1‑149 (MQ=255)
             ttGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGc                                                                                                                                      >  1:182880/1‑149 (MQ=255)
                cccGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCGCCGATTTCGGTggggct                                                                                                                                    >  2:160794/1‑145 (MQ=255)
                cccGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAg                                                                                                                                  <  7:352996/149‑1 (MQ=255)
                 ccGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGTGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAg                                                                                                                                  >  5:284062/1‑148 (MQ=255)
                 ccGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCGCCGATTTCGGTggggctga                                                                                                                                  <  1:160794/149‑6 (MQ=255)
                 ccGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcg                                                                                                                                            >  1:282419/1‑139 (MQ=255)
                 ccGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcg                                                                                                                                            <  2:282419/139‑1 (MQ=255)
                    gTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAAt                                                                                                                             <  5:344121/149‑1 (MQ=255)
                    gTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAAt                                                                                                                             <  4:361745/149‑1 (MQ=255)
                     tGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAAt                                                                                                                             >  2:11942/1‑148 (MQ=255)
                       gcggcATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGaa                                                                                                                          <  7:165060/149‑1 (MQ=255)
                          gcATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAgg                                                                                                                        >  3:93999/1‑148 (MQ=255)
                          gcATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAgg                                                                                                                        <  4:93999/148‑1 (MQ=255)
                          gcATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGCAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGt                                                                                                                       <  4:74745/149‑1 (MQ=255)
                                gCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTc                                                                                                                <  2:182880/149‑1 (MQ=255)
                                  ctctACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCg                                                                                                               >  3:270615/1‑148 (MQ=255)
                                     tACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTc                                                                                                             <  2:22402/147‑1 (MQ=255)
                                        tGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGaa                                                                                                        >  4:198263/1‑149 (MQ=255)
                                          gcggcTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGaa                                                                                                        <  7:64563/147‑1 (MQ=255)
                                              cTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGttt                                                                                                  >  3:17793/1‑149 (MQ=255)
                                                 ggTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTc                                                                                                 <  3:198263/147‑1 (MQ=255)
                                                  gTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAg                                                                                              >  3:254272/1‑149 (MQ=255)
                                                   tgcgtaCCTGCACCGCAATGCCTAACACGCCTATGGCATGTAACACGCAAGCTGCCAATGAACCGGCCTGCGTATGTACCGCCACAGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAg                                                                                              <  1:111432/144‑1 (MQ=255)
                                                   tGCCGACCTGCACCGCCAAGGCTAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGc                                                                                             <  4:270615/149‑1 (MQ=255)
                                                   tGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGc                                                                                             <  4:86241/149‑1 (MQ=255)
                                                    gCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGc                                                                                             >  2:71382/1‑148 (MQ=255)
                                                       gACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGt                                                                                         <  4:272133/149‑1 (MQ=255)
                                                          cTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAgg                                                                                      >  3:90150/1‑149 (MQ=255)
                                                           tGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGc                                                                                     >  7:122564/1‑149 (MQ=255)
                                                            gCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCt                                                                                    >  8:258757/1‑149 (MQ=255)
                                                             gaCCGCCATGGCGAACACGCCGATGGGCATGAACACGGCATGAGCCAATGTACCGGCCTGCGGAAGTACCTGCACGCTCAGCCCCACCGAAATCGGCGGGGCAGCACATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCt                                                                                    <  1:71289/148‑1 (MQ=255)
                                                             cACCGCCAAGGCGAACACGCCGATGGGCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTg                                                                                   >  2:159938/1‑149 (MQ=255)
                                                              aCCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGt                                                                                  >  1:110030/1‑149 (MQ=255)
                                                               ccGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTc                                                                                 >  6:250976/1‑149 (MQ=255)
                                                                cGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCt                                                                                >  7:196282/1‑149 (MQ=255)
                                                                  ccAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCATGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGCGTAGGCTGTCtt                                                                               <  1:91679/149‑1 (MQ=255)
                                                                        cGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAAAGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTgc                                                                        <  1:97523/149‑1 (MQ=255)
                                                                         gAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTgcg                                                                       <  3:188212/149‑1 (MQ=255)
                                                                                 cGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCCGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCcg                                                               <  6:383986/149‑1 (MQ=255)
                                                                                      gCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcggcg                                                                                                                                         <  6:63831/73‑1 (MQ=255)
                                                                                      gCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcggcg                                                                                                                                         >  5:63831/1‑73 (MQ=255)
                                                                                      gCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGATTTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGt                                                           <  1:102491/149‑1 (MQ=255)
                                                                                         aGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAAt                                                       >  6:270525/1‑149 (MQ=255)
                                                                                             acacGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATcgcg                                                   <  5:250976/149‑1 (MQ=255)
                                                                                              cacGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCgg                                                  <  7:120531/149‑1 (MQ=255)
                                                                                               acGGCAGCAGCCAATGTACCGGCCTGCGGATGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGtt                                                                          >  1:29744/1‑124 (MQ=255)
                                                                                               acGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGtt                                                                          <  2:29744/124‑1 (MQ=255)
                                                                                               acGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAGTCGCGGt                                                 >  3:363285/1‑149 (MQ=255)
                                                                                                  gcagcaGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTc                                                >  4:26546/1‑147 (MQ=255)
                                                                                                        gcCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgcgcgc                                        >  3:257769/1‑149 (MQ=255)
                                                                                                        gcCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgcgcgc                                        >  6:301486/1‑149 (MQ=255)
                                                                                                         ccAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgcgcgcg                                       <  1:11942/149‑1 (MQ=255)
                                                                                                         ccAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgcgcgcg                                       >  3:229518/1‑149 (MQ=255)
                                                                                                           aaTGTACCGGCCTCCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGCTAAATCGCGGTCCAgcgcgcgcg                                     >  7:233264/1‑149 (MQ=255)
                                                                                                            aTGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgcgcgcgcg                                     >  7:195018/1‑148 (MQ=255)
                                                                                                            aTGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGGTTCCAGCTAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCTGTAAATCGCGGTCCAgcgcgcgcgc                                    <  8:196282/149‑1 (MQ=255)
                                                                                                               tACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCAcc                                 <  1:159938/149‑1 (MQ=255)
                                                                                                                  cGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATc                              <  4:17793/149‑1 (MQ=255)
                                                                                                                  cGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATc                              <  1:71382/149‑1 (MQ=255)
                                                                                                                      cTGCGGAAGTACCGCCACGCTCAGCCGCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGcc                           >  1:200884/1‑148 (MQ=255)
                                                                                                                        gCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGc                                                                                             >  7:163481/1‑80 (MQ=255)
                                                                                                                        gCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGc                                                                                             <  8:163481/80‑1 (MQ=255)
                                                                                                                        gCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCGTCGCCGa                         >  2:258341/1‑148 (MQ=255)
                                                                                                                               tACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTgg                 <  5:270525/149‑1 (MQ=255)
                                                                                                                                  cGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCa              >  1:326715/1‑149 (MQ=255)
                                                                                                                                   gCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCAt             >  7:178717/1‑149 (MQ=255)
                                                                                                                                     cACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGa           <  4:254272/149‑1 (MQ=255)
                                                                                                                                              ggcccGCCGATTTCGGTGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGc    >  2:158954/3‑148 (MQ=255)
                                                                                                                                               gcccGCCGATTTCGGTGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGa  <  1:158954/148‑1 (MQ=255)
                                                                                                                                               ccccACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGa  <  4:257769/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGACATGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGA  >  NC_002947/4470562‑4470852

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: