Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,772,819 C→G 8.6% V48V (GTG→GTC) ‡ betT‑I ← choline transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,772,8190CG 8.6% 199.4 / 9.8 70V48V (GTG→GTC) ‡betT‑Icholine transporter
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base G (2/4);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

GCTGAAGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAGGCAGAACTCATTGCA  >  NC_002947/5772678‑5772966
                                                                                                                                             |                                                                                                                                                   
gCTGAAGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCa                                                                                                                                              >  6:337957/1‑149 (MQ=255)
 cTGAAGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAg                                                                                                                                             >  1:220681/1‑149 (MQ=255)
  tGAAGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCGTGGCCCAg                                                                                                                                             >  2:89221/1‑148 (MQ=255)
   gAAGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCa                                                                                                                                              <  7:40255/146‑1 (MQ=255)
   gAAGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCa                                                                                                                                              >  8:40255/1‑146 (MQ=255)
    aaGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCGTGGCCCAgaa                                                                                                                                           <  5:182998/148‑3 (MQ=255)
    aaGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCGTGGCCCAgaa                                                                                                                                           <  1:89221/148‑3 (MQ=255)
    aaGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGc                                                                                                                                            <  2:220681/147‑1 (MQ=255)
      gTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCGTGGCCCAgaagtg                                                                                                                                        <  5:291372/149‑6 (MQ=255)
       tCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTg                                                                                                                                       >  3:10434/1‑149 (MQ=255)
       tCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTg                                                                                                                                       >  7:180675/1‑149 (MQ=255)
         aGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGc                                                                                                                                     >  6:169672/1‑149 (MQ=255)
         aGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGc                                                                                                                                     <  5:289953/149‑1 (MQ=255)
          ggCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGcc                                                                                                                                    >  7:86803/1‑149 (MQ=255)
           gCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGGCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCcg                                                                                                                                   >  8:338154/1‑149 (MQ=255)
                  gtcCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCGTGGCCCAgaa                                                                                                                                           >  8:179860/1‑132 (MQ=255)
                  gtcCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCGTGGCCCAgaa                                                                                                                                           <  7:179860/134‑3 (MQ=255)
                  gtcCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGttt                                                                                                                            >  7:53558/1‑149 (MQ=255)
                  gtcCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGttt                                                                                                                            >  6:24782/1‑149 (MQ=255)
                     cTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCgg                                                                                                                         <  5:132059/149‑1 (MQ=255)
                        gCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTa                                                                                                                      >  7:100571/1‑149 (MQ=255)
                         ccACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTa                                                                                                                      >  5:317007/1‑148 (MQ=255)
                            cccAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCGCAGCACCTGGCCGGCGGTTTCGGGTACAGc                                                                                                                  <  5:337957/149‑1 (MQ=255)
                            cccAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGc                                                                                                                  >  3:274299/1‑149 (MQ=255)
                                    gAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGTGGTTTCGGGTACAGCGATCAGc                                                                                                           >  6:203799/1‑148 (MQ=255)
                                     aGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAg                                                                                                                   >  6:271921/1‑139 (MQ=255)
                                     aGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAg                                                                                                                   <  5:271921/139‑1 (MQ=255)
                                         tGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCTGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAgg                                                                                                       >  2:298635/1‑147 (MQ=255)
                                             gTAGTCGGAGAAAGCCAGGTAGACGAAGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCa                                                                                                 <  8:100571/149‑1 (MQ=255)
                                                gTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGaa                                                                                              <  4:313237/149‑1 (MQ=255)
                                                gTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGaa                                                                                              <  8:86803/149‑1 (MQ=255)
                                                            cAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACa                                                                                                                    <  4:38472/115‑1 (MQ=255)
                                                            cAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACa                                                                                                                    >  3:38472/1‑115 (MQ=255)
                                                            cAGGTAGACGACGAACAGCAGGTATCCGCAGATGACAAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCACGGCAGTCAGAACAAGCAtcatca                                                                                  >  8:297896/1‑149 (MQ=255)
                                                                     gacgaACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAgg                                                                         >  2:271649/1‑149 (MQ=255)
                                                                              cagGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACTCCCAGCACCTGGCCGGTGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAaca                                                                <  5:203799/149‑1 (MQ=255)
                                                                                      cGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCgg                                                        >  1:149519/1‑149 (MQ=255)
                                                                                      cGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCgg                                                        >  3:329498/1‑149 (MQ=255)
                                                                                            tGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGAcc                                                  <  2:7664/149‑1 (MQ=255)
                                                                                               ccAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGt                                               <  5:169672/149‑1 (MQ=255)
                                                                                                  gCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTAcc                                            >  5:296452/1‑149 (MQ=255)
                                                                                                         gtaCCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTcg                                      <  4:10434/148‑1 (MQ=255)
                                                                                                         gtaCCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTcg                                      <  4:274299/148‑1 (MQ=255)
                                                                                                         ctaCCAGCCAAACGTGCGCGTCAGCCACTTCTGGGCCACTCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTcg                                      <  8:180675/147‑1 (MQ=255)
                                                                                                           aCCAGCCACAGGTGCGCGGCAGCCACTTCTGGGCAACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTcggcg                                   <  8:327992/149‑1 (MQ=255)
                                                                                                            ccAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGatc                                                                                      >  2:279717/1‑97 (MQ=255)
                                                                                                            ccAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTcggcgg                                  <  6:38943/149‑1 (MQ=255)
                                                                                                            ccAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGAGCATCAGCGCGGAGGTGAAAAACCCCCCGCGGTTGACCCGTACCCTCgcggcg                                   >  7:327992/1‑148 (MQ=255)
                                                                                                             cAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCTGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACACGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCggg                                 <  5:282688/149‑1 (MQ=255)
                                                                                                             cAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACCAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCggg                                 <  1:298635/149‑1 (MQ=255)
                                                                                                             cAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAACACCACGCGGTTGACCCGTACCCTCTCGGCgggg                                >  3:35286/1‑149 (MQ=255)
                                                                                                              aGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCgggg                                <  6:317007/149‑1 (MQ=255)
                                                                                                               gCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCggggg                               <  5:24782/149‑1 (MQ=255)
                                                                                                                  aaaGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGTGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGCGGGCt                            >  8:284118/1‑149 (MQ=255)
                                                                                                                  aaaGGTGCGCGTCAGCAACTTCTGGGCCACGCCCAGCAGCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAAAATGGATCATCAGCGCGGAGGTGAAAAACCCCACGCGGTTGACCCGTCCCCGGTCGGCGGGGGGcg                            >  6:272517/1‑148 (MQ=255)
                                                                                                                   aaGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCtt                           >  2:32788/1‑149 (MQ=255)
                                                                                                                      gTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGCACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTgg                         >  5:250018/1‑148 (MQ=255)
                                                                                                                      gTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGt                        <  8:53558/149‑1 (MQ=255)
                                                                                                                         cgcgTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGt                    <  4:35286/149‑1 (MQ=255)
                                                                                                                          gcgTCAGCCACTTCTTGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGt                    <  6:296452/149‑1 (MQ=255)
                                                                                                                          gcgTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGt                    >  7:292740/1‑149 (MQ=255)
                                                                                                                                gCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAg                                                                                                >  6:191158/1‑67 (MQ=255)
                                                                                                                                gCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAGGCa              <  2:338597/149‑1 (MQ=255)
                                                                                                                                 ccACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAgg                >  5:155753/1‑146 (MQ=255)
                                                                                                                                  cACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACACTCTCGGCGGGGGGCTTGGTAAGTGAGGCAg             >  6:94249/1‑148 (MQ=255)
                                                                                                                                   aCTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCATTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAGGCAGaa           <  7:338154/149‑1 (MQ=255)
                                                                                                                                   aCTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAGGCAGaa           <  1:173639/149‑1 (MQ=255)
                                                                                                                                          ggCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAGGCAGAACTCAtt     <  7:297896/148‑1 (MQ=255)
                                                                                                                                           gCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTa                                              >  6:311203/1‑106 (MQ=255)
                                                                                                                                           gCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTa                                              <  5:311203/106‑1 (MQ=255)
                                                                                                                                            ccACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAGGCAGAACTCATTGCa  >  5:35450/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCTGAAGTCAGGCTGGTCGTCCTTGCCACCCAGCTTGAGGCTGCCGTAGTCGGAGAAAGCCAGGTAGACGACGAACAGCAGGTATCCGCAGATGACCAGCATGTAGTACCAGCCAAAGGTGCGCGTCAGCCACTTCTGGGCCACGCCCAGCACCTGGCCGGCGGTTTCGGGTACAGCGATCAGCAAGGCAGTCAGAACAAGGATCATCAGCGCGGAGGTGAAAAACACCACGCGGTTGACCCGTACCCTCTCGGCGGGGGGCTTGGTAAGTGAGGCAGAACTCATTGCA  >  NC_002947/5772678‑5772966

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: