Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,599,307 G→T 9.5% R168R (CGG→AGG)  PP_4071 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,599,3070GT 9.5% 237.8 / 13.8 84R168R (CGG→AGG) PP_4071membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base T (3/5);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

GACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAGG  >  NC_002947/4599159‑4599455
                                                                                                                                                     |                                                                                                                                                    
gACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                                      <  3:422620/149‑1 (MQ=255)
     tACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGcccc                                                                                                                                                 >  2:358462/1‑149 (MQ=255)
       ccGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGt                                                                                                                                               >  6:404695/1‑149 (MQ=255)
        cGAAGACTCAACGTTAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                                              >  8:83633/1‑149 (MQ=255)
        cGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                                              >  1:411301/1‑149 (MQ=255)
         gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGCTAAAGGCACCGGGCGGGGCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgg                                                                                                                                             <  4:282044/149‑2 (MQ=37)
         gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGc                                                                                                                                             <  4:391401/149‑1 (MQ=255)
          aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGcc                                                                                                                                            >  6:374378/1‑149 (MQ=255)
          aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTggg                                                                                                                                            <  4:77446/149‑3 (MQ=38)
           aGACTCAACGCCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgggg                                                                                                                                           <  7:35316/149‑4 (MQ=37)
             aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                         >  8:8978/1‑149 (MQ=255)
             aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                         >  8:2671/1‑149 (MQ=255)
             aCTCAACGTCAGGGGAGGCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                                                         <  3:151650/149‑1 (MQ=255)
              cTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcg                                                                                                                                        <  2:35839/149‑1 (MQ=255)
               tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                       >  7:20697/1‑149 (MQ=255)
               tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                       >  5:201125/1‑149 (MQ=255)
               tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                                 >  7:191986/1‑139 (MQ=38)
               tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                                                 <  8:191986/139‑1 (MQ=38)
                cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                                      <  7:2671/149‑1 (MQ=255)
                cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                                      <  7:8978/149‑1 (MQ=255)
                 aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                                       <  6:378962/147‑1 (MQ=255)
                 aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                                     >  7:105538/1‑149 (MQ=255)
                 aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                                     >  5:52232/1‑149 (MQ=255)
                 aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                                     >  5:332642/1‑149 (MQ=255)
                           gAGTCACTGCCCAGCGAATCTGTGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCg                                                                                                                           >  6:12035/1‑149 (MQ=255)
                           gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGc                                                                                                                            <  6:81427/148‑1 (MQ=255)
                              tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                                        >  8:51009/1‑149 (MQ=255)
                              tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                                                        <  5:54930/149‑1 (MQ=255)
                                  tGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                                    >  7:214131/1‑149 (MQ=255)
                                    cccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                                    >  7:377068/1‑147 (MQ=255)
                                      cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                                                 >  7:132207/1‑148 (MQ=255)
                                      cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACAAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                                                 >  7:339325/1‑148 (MQ=255)
                                       aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                                >  5:122338/1‑148 (MQ=255)
                                        gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                                              <  6:122338/149‑1 (MQ=255)
                                            aTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACttt                                                                                                          <  2:411301/149‑1 (MQ=255)
                                                 cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTg                                                                                                         <  5:441286/145‑1 (MQ=255)
                                                 cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATc                                                                                                      <  6:201125/148‑1 (MQ=255)
                                                    ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  <  5:12035/149‑1 (MQ=255)
                                                    ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  >  2:69170/1‑149 (MQ=255)
                                                    ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  <  6:52232/149‑1 (MQ=255)
                                                     gCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                                                  >  8:31085/1‑148 (MQ=255)
                                                      ccGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCa                                                                                                <  6:332642/149‑1 (MQ=255)
                                                       cGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAg                                                                                               <  7:51009/149‑1 (MQ=255)
                                                         ttCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTc                                                                                             <  7:83633/149‑1 (MQ=255)
                                                             tgcagATAAACGCACCGGGGGGGTCGATGATCCCCGCATCATTGGTGATACACCCCTCCCTCGCGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCaa                                                                                           <  1:358462/144‑1 (MQ=255)
                                                                  aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                                       <  6:234841/130‑1 (MQ=255)
                                                                  aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                                       >  5:234841/1‑130 (MQ=255)
                                                                     aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                                      <  4:321757/80‑1 (MQ=255)
                                                                     aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                                      >  3:321757/1‑80 (MQ=255)
                                                                     aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgctgct                                                                                 <  3:212510/149‑1 (MQ=255)
                                                                      aGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGctc                                                                                <  1:231808/149‑1 (MQ=255)
                                                                            cgggcgggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCg                                                                          <  7:31085/149‑1 (MQ=255)
                                                                                  ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                                    >  4:425101/1‑149 (MQ=255)
                                                                                  ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                                    <  6:325981/149‑1 (MQ=255)
                                                                                  ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                                                    <  5:374378/149‑1 (MQ=255)
                                                                                        tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                                              <  2:6542/149‑1 (MQ=255)
                                                                                        tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                                              <  8:20697/149‑1 (MQ=255)
                                                                                         gatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcgc                                                             >  6:313327/1‑149 (MQ=255)
                                                                                           tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                                 <  2:401858/143‑1 (MQ=255)
                                                                                           tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                                 >  1:401858/1‑143 (MQ=255)
                                                                                            cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                                >  2:44751/1‑95 (MQ=255)
                                                                                            cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                                <  1:44751/95‑1 (MQ=255)
                                                                                             aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTATGCGTGCGCCTTCGCTGGGCGGGCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCa                                                          <  2:299253/148‑1 (MQ=255)
                                                                                             aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAg                                                         >  4:247087/1‑149 (MQ=255)
                                                                                                   tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATGAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                                                   <  8:48266/149‑1 (MQ=255)
                                                                                                   tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                                                   >  7:257281/1‑149 (MQ=255)
                                                                                                    catTGCTGATCCACCAC‑CCCAGGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGGCAACGAGAACCTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGCgg                                                  >  8:386002/1‑149 (MQ=255)
                                                                                                     atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAgg                                                  <  3:421283/148‑1 (MQ=255)
                                                                                                     atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGt                                                 <  8:257281/149‑1 (MQ=255)
                                                                                                         cTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCa                                             >  7:300901/1‑149 (MQ=255)
                                                                                                                 ccac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATc                                     <  5:270398/149‑1 (MQ=255)
                                                                                                                  cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                                    >  7:352522/1‑149 (MQ=255)
                                                                                                                  cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                                    >  4:430497/1‑149 (MQ=255)
                                                                                                                   ac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCgg                                   >  6:135943/1‑149 (MQ=255)
                                                                                                                    c‑cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATATAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGt                                  <  4:107898/149‑1 (MQ=255)
                                                                                                                      cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa                                 <  8:132207/149‑1 (MQ=255)
                                                                                                                       gcACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa                                 >  1:99909/2‑148 (MQ=38)
                                                                                                                           ggggccggCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcg                            >  4:214933/7‑149 (MQ=38)
                                                                                                                           ggggTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcg                            <  8:339325/149‑1 (MQ=255)
                                                                                                                                 ggCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGAc                      <  3:247087/149‑1 (MQ=255)
                                                                                                                                         ccTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg              >  5:286600/1‑149 (MQ=255)
                                                                                                                                         ccTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg              >  2:91759/1‑149 (MQ=255)
                                                                                                                                         caggcgcTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg              <  2:99909/146‑1 (MQ=255)
                                                                                                                                          cTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCgg             >  5:169797/1‑149 (MQ=255)
                                                                                                                                              gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTTCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt         <  8:214131/149‑1 (MQ=255)
                                                                                                                                              gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGGCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt         >  7:274543/1‑149 (MQ=255)
                                                                                                                                              gcTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCAtt         >  5:167623/1‑149 (MQ=255)
                                                                                                                                                    cggccCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAg   <  8:377068/149‑1 (MQ=255)
                                                                                                                                                     ggccCCGTGCCTGAGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAgg  <  8:105538/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGCATTGTCCAGG  >  NC_002947/4599159‑4599455

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: ‑