Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,059,505 A→T 8.6% N443K (AAT→AAA PP_4458 ← opine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,059,5050AT 8.6% 235.4 / 10.4 81N443K (AAT→AAAPP_4458opine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base T (3/4);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCG  >  NC_002947/5059360‑5059644
                                                                                                                                                 |                                                                                                                                            
ccTGTTTCCCACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTc                                                                                                                                           <  2:360341/149‑1 (MQ=255)
ccTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTc                                                                                                                                           <  4:130030/149‑1 (MQ=255)
ccTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTc                                                                                                                                           <  3:248187/149‑1 (MQ=255)
  tGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCGTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCg                                                                                                                                         >  7:170160/1‑148 (MQ=255)
    tttCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGa                                                                                                                                        <  7:176967/147‑1 (MQ=255)
      tCATACCAAACAATGGGATGTCGCATACGTCGTCGTGAATCAGACCCTGCGCCATTTTCCTCAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGACCAATGTCCGAAATTCGACAt                                                                                                                                     <  2:106245/148‑1 (MQ=255)
        aGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATGATCAGCTACTGAGCAATTTCCGAAATTCGACATTg                                                                                                                                   <  6:350397/148‑1 (MQ=255)
         gACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGct                                                                                                                                 <  6:173044/149‑2 (MQ=255)
          aCCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctc                                                                                                                                <  7:37735/149‑3 (MQ=255)
          aCCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctc                                                                                                                                <  5:316119/149‑3 (MQ=255)
           ccAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGc                                                                                                                               >  6:295406/1‑149 (MQ=255)
            caaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCg                                                                                                                              >  8:112552/1‑149 (MQ=255)
            caaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCg                                                                                                                              >  3:104091/1‑149 (MQ=255)
             aaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCgg                                                                                                                             <  2:264492/149‑1 (MQ=255)
              aacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCggg                                                                                                                            >  5:96881/1‑149 (MQ=255)
                   aGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTgcgc                                                                                                                       <  6:164956/149‑1 (MQ=255)
                     ggATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGcc                                                                                                                      >  3:175844/1‑148 (MQ=255)
                      gATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGAt                                                                                                                   <  1:139592/149‑1 (MQ=255)
                        tGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATg                                                                                                                  >  4:254819/1‑149 (MQ=255)
                          gcgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATgg                                                                                                                 >  6:105146/1‑148 (MQ=255)
                          gcgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATgg                                                                                                                 >  6:308022/1‑148 (MQ=255)
                                 cgtcgtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGc                                                                                                          <  4:198548/148‑1 (MQ=255)
                                      tcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATc                                                                                                    >  8:118862/1‑149 (MQ=255)
                                        gtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCg                                                                                                   >  4:418922/1‑148 (MQ=255)
                                         tGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTa                                                                                                 <  2:262565/149‑1 (MQ=255)
                                           aaTCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTaa                                                                                                <  3:46031/148‑1 (MQ=255)
                                           aaTCTGAACCTGCGCCTTTTTCCACAAATCGACTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCTGATGGCATCATCAGCTACTGATCAAAGTCCGAAATTCGACATTCAAGCGGGTGCTCCGATGACAGAAGCCGTATCGtca                                                                                                >  7:168889/1‑146 (MQ=255)
                                            aTCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTaaaa                                                                                              <  3:254819/149‑1 (MQ=255)
                                            aTCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTaaa                                                                                               >  7:362347/1‑148 (MQ=255)
                                                gAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                          >  6:254210/1‑149 (MQ=255)
                                                 aaCCTGCGCCTTTTTCCACAAATCGAGATGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                          <  5:118000/148‑1 (MQ=255)
                                                  aCCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGt                                                                                        <  3:248744/149‑1 (MQ=255)
                                                   ccTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTc                                                                                       <  4:370045/149‑1 (MQ=255)
                                                   ccTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTc                                                                                       <  4:376155/149‑1 (MQ=255)
                                                     tGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAg                                                                                     >  7:351790/1‑149 (MQ=255)
                                                       cgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                    >  2:247674/1‑148 (MQ=255)
                                                         ccTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCag                                                                                 <  6:96881/149‑1 (MQ=255)
                                                         ccTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCag                                                                                 <  7:112552/149‑1 (MQ=255)
                                                             tttCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCTACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCagta                                                                               <  5:308022/147‑1 (MQ=255)
                                                             tttCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCagta                                                                               <  5:295406/147‑1 (MQ=255)
                                                             tttCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCagta                                                                               <  5:105146/147‑1 (MQ=255)
                                                                ccACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCgg                                                                          >  6:119903/1‑149 (MQ=255)
                                                                 cacaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgc                                                                         >  3:296194/1‑149 (MQ=255)
                                                                  acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTTGgcg                                                                        >  3:348696/1‑149 (MQ=255)
                                                                  acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGAAGAAGCCGAATCGTAAAACTCTGTCACCCAGTAGTCGgcg                                                                        >  6:111734/1‑149 (MQ=255)
                                                                     aaTCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGGTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAACCTCTGTCAGCCAGGAGTCGGCGccc                                                                     >  4:111680/1‑149 (MQ=255)
                                                                     aaTCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGccc                                                                     >  7:36252/1‑149 (MQ=255)
                                                                       tCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCgg                                                                   <  8:362347/149‑1 (MQ=255)
                                                                           gTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCAGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAATCCGAATCGTAAAACACTGTCAGCCAGTAGTCGGCGCCGGGatat                                                               >  3:101014/1‑149 (MQ=255)
                                                                            ttGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAACCTCTGTCTGCCAGTAGTCGGCGCCCGGATATc                                                              >  7:146171/1‑149 (MQ=255)
                                                                               aGTCTTCTGATCATCCTCGCCCATGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACCCTGTCAGCCATTAGTCGGCGCCCGGATATCGag                                                           <  4:101014/149‑1 (MQ=255)
                                                                               aGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGACCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCTAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                           >  5:92519/1‑149 (MQ=255)
                                                                                gTCTGCGGATAAGCCGCGAACCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAAATCGACATTGCAGCGCGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCATCCAGTAGTCGGCGCCCGGATATCGag                                                           <  8:170160/148‑1 (MQ=255)
                                                                                    gCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAgcagca                                                      >  4:418925/1‑149 (MQ=255)
                                                                                     cGGATCAGCCTAGACCCTTGCTTTACGTATGGAATCATCAGCTAGTGAGCAATTTCCGAAATTCGACATTGCAGCGGGTGCGCCTATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAgcagca                                                      <  8:168889/148‑1 (MQ=255)
                                                                                      ggctCAGCCTCGACCCTGGCGTTCCGGATGGCATCAGCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTTGGCGCCCGGATATCGAGCAGCAcc                                                    <  4:348696/146‑1 (MQ=255)
                                                                                      ggATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCAcc                                                    <  5:119903/149‑1 (MQ=255)
                                                                                            gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa                                              <  4:104091/149‑1 (MQ=255)
                                                                                             ccTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCAcc                                                    >  7:49051/1‑142 (MQ=255)
                                                                                             ccTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCAcc                                                    <  8:49051/142‑1 (MQ=255)
                                                                                                   cccTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGTGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGTATATCGAGCAGCACCGTAAAACACGATc                                       <  1:247674/149‑1 (MQ=255)
                                                                                                       ggCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCgcgc                                   >  3:286583/1‑149 (MQ=255)
                                                                                                         cTTTCCGGATGGCATCATCAGCTACAGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGctc                                 >  2:110298/1‑149 (MQ=255)
                                                                                                          tttCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGctct                                >  7:225294/1‑149 (MQ=255)
                                                                                                            tCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTgg                              >  2:346391/1‑149 (MQ=255)
                                                                                                                  tGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGa                        <  3:418925/149‑1 (MQ=255)
                                                                                                                     catcatCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACgg                     >  1:68171/1‑149 (MQ=255)
                                                                                                                       tcatcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCt                   >  7:35910/1‑149 (MQ=255)
                                                                                                                           caGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCgcgc                                   >  7:197205/1‑129 (MQ=255)
                                                                                                                           caGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCgcgc                                   <  8:197205/129‑1 (MQ=255)
                                                                                                                            aGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGTTCTTGACGGCTTTTg               <  8:36252/148‑1 (MQ=255)
                                                                                                                            aGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTgg              >  3:349105/1‑149 (MQ=255)
                                                                                                                             gCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCCCCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGc             >  7:136223/1‑149 (MQ=255)
                                                                                                                               tACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCt           >  5:213081/1‑149 (MQ=255)
                                                                                                                                  gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTgg         >  7:4574/4‑148 (MQ=255)
                                                                                                                                  gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTgg         >  2:212414/4‑148 (MQ=255)
                                                                                                                                   gAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt        >  7:110594/1‑148 (MQ=255)
                                                                                                                                    tgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt        >  2:98073/2‑147 (MQ=255)
                                                                                                                                    aGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt        <  3:418922/147‑1 (MQ=255)
                                                                                                                                    aGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTaa      >  1:427009/1‑149 (MQ=255)
                                                                                                                                         atgtcGAATTTCGGACATTGCAGCGGGTGCTCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTcg  <  8:4574/145‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCG  >  NC_002947/5059360‑5059644

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: