Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,777,525 C→G 7.6% intergenic (+259/+33) PP_5476 → / ← PP_1585 hypothetical protein/antidote protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,777,5250CG 7.6% 240.0 / 8.0 79intergenic (+259/+33)PP_5476/PP_1585hypothetical protein/antidote protein
Reads supporting (aligned to +/- strand):  ref base C (29/44);  new base G (2/4);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GTCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCA  >  NC_002947/1777382‑1777670
                                                                                                                                                |                                                                                                                                                 
gTCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTgg                                                                                                                                              >  7:209582/1‑149 (MQ=255)
 tCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGgc                                                                                                                                             <  5:55118/149‑1 (MQ=255)
      ttAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGcccc                                                                                                                                        >  5:63251/1‑149 (MQ=255)
        aaCTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGGGCGTGGCGccgcg                                                                                                                                       <  6:145467/148‑4 (MQ=255)
         aCTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGGGCGTGGCGccgcga                                                                                                                                      <  1:117745/148‑5 (MQ=255)
             ggCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTa                                                                                                                                  <  4:82312/148‑1 (MQ=255)
             ggCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTa                                                                                                                                  <  5:343257/148‑1 (MQ=255)
             ggCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCACGCGGCGCCACGCCCGTGGCGCCCCTGttt                                                                                                                                   <  1:32737/147‑1 (MQ=255)
              gCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTAtt                                                                                                                                <  4:315091/149‑1 (MQ=255)
              gCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTAt                                                                                                                                 <  1:81852/148‑1 (MQ=255)
                  cccAGGCTTTGATGCCGTAGAGAGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGTGGGGCACCTGTTTATTGAt                                                                                                                             >  8:227205/1‑148 (MQ=255)
                   ccAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGAt                                                                                                                             >  4:273021/1‑147 (MQ=255)
                    cAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAg                                                                                                                          >  1:25988/1‑149 (MQ=255)
                     aGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAg                                                                                                                          >  1:358690/1‑148 (MQ=255)
                      ggCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATc                                                                                                                            <  5:59640/145‑1 (MQ=255)
                      ggCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATc                                                                                                                            >  6:59640/1‑145 (MQ=255)
                        cTTTGAAGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAt                                                                                                                       <  7:87762/148‑1 (MQ=255)
                           tGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAtgt                                                                                                                     <  7:65709/147‑1 (MQ=255)
                               gCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAtgt                                                                                                                     >  7:107056/1‑143 (MQ=255)
                               gCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAtgt                                                                                                                     <  8:107056/143‑1 (MQ=255)
                                ccGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCaa                                                                                                               <  6:264478/148‑1 (MQ=255)
                                 cGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGc                                                                                                             <  8:209582/149‑1 (MQ=255)
                                     gagCGTGTAGGTCATTGGTTCCATTGGGGTTGCAACTAGCCCCAATCTCTCCGTTTTTCGTAGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGAc                                                                                                          <  8:154202/148‑1 (MQ=255)
                                        cGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGAAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGa                                                                                                      >  2:289429/1‑149 (MQ=255)
                                               gTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGga                                                                                               >  5:123375/1‑149 (MQ=255)
                                                  aTTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCaa                                                                                                               >  8:44989/1‑130 (MQ=255)
                                                  aTTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCaa                                                                                                               <  7:44989/130‑1 (MQ=255)
                                                   ttGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccga                                                                                           <  3:131888/149‑1 (MQ=255)
                                                     gtgtACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccgacc                                                                                         >  5:110782/1‑149 (MQ=255)
                                                        tACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCga                                                                                      >  3:239495/1‑149 (MQ=255)
                                                        tACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCga                                                                                      <  3:230912/149‑1 (MQ=255)
                                                          cATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCgatg                                                                                    <  3:273021/149‑1 (MQ=255)
                                                          cATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGAAACGATCCTCGCGGCGCCACGACCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCga                                                                                      <  1:289429/148‑1 (MQ=255)
                                                           aTTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCgatga                                                                                   <  4:239495/149‑1 (MQ=255)
                                                            ttGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCTCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCgatgat                                                                                  <  2:330200/149‑1 (MQ=255)
                                                              ggggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATAAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGAtct                                                                                >  3:60690/1‑149 (MQ=255)
                                                                     cttgcaaccccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGcc                                                                         <  5:261053/142‑1 (MQ=255)
                                                                     cAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGcc                                                                         <  6:110782/149‑1 (MQ=255)
                                                                      aaGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCg                                                                        <  6:123375/149‑1 (MQ=255)
                                                                         tAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTTTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCgtgt                                                                     <  6:11227/149‑1 (MQ=255)
                                                                         tAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAt                                                                                                                       <  6:294615/99‑1 (MQ=255)
                                                                         tAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAt                                                                                                                       >  5:294615/1‑99 (MQ=255)
                                                                            ccccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGc                                                                  >  1:23221/1‑149 (MQ=255)
                                                                             cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGAGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCt                                                                 >  2:275216/1‑149 (MQ=255)
                                                                              ccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGTGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGTCCAGACGCTGGACCGACCCGATGAtctc                                                                               <  1:157399/134‑1 (MQ=255)
                                                                              ccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCg                                                                                                            <  8:229094/105‑1 (MQ=255)
                                                                              ccAATCTCTCCGTTTTTCGATGGCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCg                                                                                                            >  7:229094/1‑105 (MQ=255)
                                                                               cAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGt                                                               >  2:118823/1‑149 (MQ=255)
                                                                               cAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGt                                                               >  6:179493/1‑149 (MQ=255)
                                                                                aaTCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGAGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGtt                                                              <  6:63251/149‑1 (MQ=255)
                                                                                    tctcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGc                                                          >  7:40317/1‑149 (MQ=255)
                                                                                     ctcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCt                                                         >  5:242744/1‑149 (MQ=255)
                                                                                      tcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCtg                                                        >  3:15883/1‑149 (MQ=255)
                                                                                      tcCGTTTTTCGATGTCAGTCGCGCTATTATGCTGGAACGATCCTCGCGGCGACCCGCCCGTGGCCGCCCTGTTTATTGAAAAGGATGTCTTCAATCGACCAGACGCTGGGCCGAACCGCTGATCTCATCGCCGTGTTTCTGTTAGGCtg                                                        >  6:4387/1‑149 (MQ=255)
                                                                                        cGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGAGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCtgt                                                       <  1:275216/148‑1 (MQ=255)
                                                                                         gTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGCCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTgc                                                     >  1:253842/1‑149 (MQ=255)
                                                                                          tttttCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTgcg                                                    <  6:242744/149‑1 (MQ=255)
                                                                                             ttCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGc                                                          <  1:152924/140‑1 (MQ=255)
                                                                                             ttCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACCCCCGTTGCGCCCCAGTTTTTTGATGAGGATGTGTTAAACCGACCCAACCCGGGACCCGCCCGATGATCTAAAAGCCGTGTTTCTGTTGGGCtaggaagg                                                  >  2:152924/1‑141 (MQ=255)
                                                                                                    tCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGaa                                           <  4:15883/148‑1 (MQ=255)
                                                                                                      aGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAgg                                         <  6:344152/148‑1 (MQ=255)
                                                                                                       gTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGcc                                       <  3:325707/149‑1 (MQ=255)
                                                                                                        tCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCg                                      >  4:329011/1‑149 (MQ=255)
                                                                                                              gATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCa                                <  8:40317/149‑1 (MQ=255)
                                                                                                                  aTGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGt                             <  2:358690/148‑1 (MQ=255)
                                                                                                                     cGGGAACGATCCTCGCGTCGGCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGGGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGa                          <  5:244940/148‑1 (MQ=255)
                                                                                                                      gggAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGc                        <  3:329011/149‑1 (MQ=255)
                                                                                                                         aaCGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGc                        >  7:91735/1‑146 (MQ=255)
                                                                                                                                 tCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAg              >  8:2722/1‑148 (MQ=255)
                                                                                                                                   gcggcgCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAgg             <  5:179493/147‑1 (MQ=255)
                                                                                                                                   gcggcgCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAg              <  7:2722/146‑1 (MQ=255)
                                                                                                                                    gggcgCCACGGGCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGggt          >  8:179113/2‑149 (MQ=255)
                                                                                                                                     ggcgCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGAACAGACGCTGGACCGACCCGATTATCTGATCGCCGTGTTGCTGTTCGGCTGTGCGCTGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGggtg         <  5:4387/149‑1 (MQ=255)
                                                                                                                                      gcgcCACGGGCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAggg            <  5:190805/145‑1 (MQ=255)
                                                                                                                                      gcgcCACGGGCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAggg            >  6:190805/1‑145 (MQ=255)
                                                                                                                                       cgcCACGGGCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGggtggt       <  7:179113/149‑1 (MQ=255)
                                                                                                                                       cgcCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGggtggt       <  2:23221/149‑1 (MQ=255)
                                                                                                                                            cgcgcgTGGCGCCCCTGTTTATTTTTCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCa  <  4:60690/149‑1 (MQ=255)
                                                                                                                                            cGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCa  <  1:118823/149‑1 (MQ=255)
                                                                                                                                            cGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCa  <  2:25988/149‑1 (MQ=255)
                                                                                                                                               cccGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgac                                                                                              <  6:35882/55‑1 (MQ=255)
                                                                                                                                               cccGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgac                                                                                              >  5:35882/1‑55 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCA  >  NC_002947/1777382‑1777670

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: