Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,168,077 T→C 5.1% F381L (TTT→CTT)  PP_1921 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,168,0770TC 5.1% 233.1 / 2.1 78F381L (TTT→CTT) PP_1921hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base C (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

ACTTTGATGCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGG  >  NC_002947/2167931‑2168222
                                                                                                                                                  |                                                                                                                                                 
aCTTTGATGCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTAtttt                                                                                                                                                 >  4:348534/1‑149 (MQ=255)
  tttGATGCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTg                                                                                                                                               >  6:73718/1‑149 (MQ=255)
  tttGATGCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTAGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAATCGAATAGCATATATTTTTg                                                                                                                                               >  4:338549/1‑149 (MQ=255)
     gATGCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTggg                                                                                                                                             <  6:278387/148‑1 (MQ=255)
      aTGCAGGCGTGTCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTggggg                                                                                                                                           <  6:141547/149‑1 (MQ=255)
       tGCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGc                                                                                                                                          >  8:289006/1‑149 (MQ=255)
        gCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATCTTTGGGGGCt                                                                                                                                         >  5:444836/1‑149 (MQ=255)
          aGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCttt                                                                                                                                       <  4:314343/149‑1 (MQ=255)
           ggCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATCTTTGGGGGCTTTa                                                                                                                                      >  8:222340/1‑149 (MQ=255)
             cGTGGCAGACTTAGCCACATACCTTGGCTTGAAATACGCTTTGCCGGGTGGCATCGCGGCACTCGTACGATACTCTAGGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGTACAAAACGAATAGCATCTATTTTTGGGGGCTTTAtg                                                                                                                                    <  6:238989/149‑1 (MQ=255)
             cGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTAtg                                                                                                                                    <  3:418815/149‑1 (MQ=255)
              gTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTAtgc                                                                                                                                   >  3:57889/1‑149 (MQ=255)
                 ccAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTAtgctg                                                                                                                                 <  2:416516/148‑1 (MQ=255)
                  cAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCg                                                                                                                               >  6:142907/1‑149 (MQ=255)
                     aCTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGCCTAATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGa                                                                                                                            >  7:94213/1‑149 (MQ=255)
                           gCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTAcgc                                                                                                                      <  5:142907/149‑1 (MQ=255)
                            cacaGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTAcgcg                                                                                                                     <  7:289006/149‑1 (MQ=255)
                             acaGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCgg                                                                                                                    >  3:226544/1‑149 (MQ=255)
                             acaGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCgg                                                                                                                    >  6:52651/1‑149 (MQ=255)
                             acaGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCgg                                                                                                                    >  2:410826/1‑149 (MQ=255)
                                gACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCAt                                                                                                                 >  8:109195/1‑149 (MQ=255)
                                   cTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATa                                                                                                                >  4:27605/1‑147 (MQ=255)
                                      ggCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATagag                                                                                                             <  3:348534/147‑1 (MQ=255)
                                        cTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTAc                                                                                                         >  3:14977/1‑149 (MQ=255)
                                        cTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTAc                                                                                                         >  3:309308/1‑149 (MQ=255)
                                         ttGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGt                                                                                                                              >  5:183667/1‑127 (MQ=255)
                                         ttGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGt                                                                                                                              <  6:183667/127‑1 (MQ=255)
                                         ttGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTAcc                                                                                                        >  2:330590/1‑149 (MQ=255)
                                           gAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATCTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCg                                                                                                       <  6:444836/148‑1 (MQ=255)
                                            aaaTACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGa                                                                                                      >  2:187843/1‑148 (MQ=255)
                                             aaTACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGaa                                                                                                     >  5:218832/1‑148 (MQ=255)
                                              aTACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGaaaa                                                                                                   >  6:265584/1‑149 (MQ=255)
                                                  gCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAgcgc                                                                                               <  7:109195/149‑1 (MQ=255)
                                                  gCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAgcgc                                                                                               >  1:399451/1‑149 (MQ=255)
                                                      tggcGGCTGGCATCGCGGCACTAGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTACAGAGGTACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGaa                                                                                            <  2:141868/145‑1 (MQ=255)
                                                      tGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGaa                                                                                            <  5:52651/148‑1 (MQ=255)
                                                           gCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGt                                                                                      <  1:410826/149‑1 (MQ=255)
                                                             tGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTaa                                                                                    >  6:354227/1‑149 (MQ=255)
                                                               gCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAat                                                                                   <  1:187843/148‑1 (MQ=255)
                                                               gCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATCTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAat                                                                                   <  7:222340/148‑1 (MQ=255)
                                                                cATCGCGGCACTTGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCATTATGATGAGTGACTACGCGGCATAGAGATACCGAAAAGCGCGAATTAAGTaa                                                                                    >  6:30108/1‑146 (MQ=255)
                                                                 aTCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAACATc                                                                                >  3:293640/1‑149 (MQ=255)
                                                                    gcgGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTa                                                                             <  8:318172/149‑1 (MQ=255)
                                                                    gcgGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTa                                                                             <  7:33167/149‑1 (MQ=255)
                                                                             gTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAACGCGCGAATTAAGTAATATCTTATTTGGGGcc                                                                    >  3:390385/1‑149 (MQ=255)
                                                                              tACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCg                                                                   <  4:40272/149‑1 (MQ=255)
                                                                                 gATACTCTACGGCTCATGGCGTATCAATTGAGTATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGTGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTAt                                                                <  4:226544/149‑1 (MQ=255)
                                                                                 gATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTAt                                                                <  4:14977/149‑1 (MQ=255)
                                                                                      tctACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCt                                                           >  8:394925/1‑149 (MQ=255)
                                                                                      tctACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCt                                                           >  4:335744/1‑149 (MQ=255)
                                                                                       ctACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCt                                                           >  3:221993/1‑148 (MQ=255)
                                                                                          cGGCTCATGTCCTATCAATTTAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTTACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTg                                                        <  8:270921/148‑1 (MQ=255)
                                                                                            gCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGaaa                                                     >  4:182086/1‑149 (MQ=255)
                                                                                                aTGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGt                                                   >  2:42087/1‑147 (MQ=255)
                                                                                                      tATCAGTTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCg                                           >  7:394060/1‑149 (MQ=255)
                                                                                                        tCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGt                                         <  4:309308/149‑1 (MQ=255)
                                                                                                          aaTTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGc                                       <  2:399451/149‑1 (MQ=255)
                                                                                                           aTTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGc                                       >  3:24158/1‑148 (MQ=255)
                                                                                                              gAGGATCTCCAGAGGGACAAAACTAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGATCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAg                                    <  3:338549/148‑1 (MQ=255)
                                                                                                                ggATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTa                                                                                         <  8:375192/93‑1 (MQ=255)
                                                                                                                ggATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTa                                                                                         >  7:375192/1‑93 (MQ=255)
                                                                                                                ggATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCag                                 <  1:132130/149‑1 (MQ=255)
                                                                                                                ggATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCa                                  <  5:73718/148‑1 (MQ=255)
                                                                                                                 taTCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCag                                 <  4:352599/147‑1 (MQ=255)
                                                                                                                 gATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCTTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCaga                                <  5:417668/149‑1 (MQ=255)
                                                                                                                    ctcCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGc                              >  2:226748/1‑148 (MQ=255)
                                                                                                                       cAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAg                           <  4:24158/148‑1 (MQ=255)
                                                                                                                          agGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCAcc                       <  1:42087/149‑1 (MQ=255)
                                                                                                                            ggACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGc                     >  2:230402/1‑149 (MQ=255)
                                                                                                                              aCAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGATAGCAGAGATAGCACCGCtt                   >  7:382027/1‑149 (MQ=255)
                                                                                                                                 aaccGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTg                <  5:265584/146‑1 (MQ=255)
                                                                                                                                 aaaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTg                >  6:96502/1‑149 (MQ=255)
                                                                                                                                 aaaCGAATAGAATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCCACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTCCTTTGAAAGTAATATGAGGTGCTAGCAGAGCTAGCACCGCTATTg                >  8:158325/1‑149 (MQ=255)
                                                                                                                                      aaTAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCAtt           <  5:354227/149‑1 (MQ=255)
                                                                                                                                        tAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATtgt         <  5:96502/149‑1 (MQ=255)
                                                                                                                                            aTCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGaaa                                                     >  2:403175/1‑101 (MQ=255)
                                                                                                                                            aTCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGaaa                                                     <  1:403175/101‑1 (MQ=255)
                                                                                                                                            aTCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCtttt                 >  8:218163/1‑137 (MQ=255)
                                                                                                                                            aTCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCtttt                 <  7:218163/137‑1 (MQ=255)
                                                                                                                                                 tttttGGGGGCTTTATGCTGCGTTACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCgg  <  4:390385/147‑1 (MQ=255)
                                                                                                                                                 tttttGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCgg  <  7:394925/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACTTTGATGCAGGCGTGCCAGACTTAGGCACAGACCTTGGCTTGAAATACGCTTTGCCGGCTGGCATCGCGGCACTCGTACGATACTCTACGGCTCATGGCGTATCAATTGAGGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGG  >  NC_002947/2167931‑2168222

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑