Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,203,150 G→A 6.9% intergenic (+373/‑100) PP_2809 → / → PP_2810 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,203,1500GA 6.9% 212.5 / 3.2 72intergenic (+373/‑100)PP_2809/PP_2810membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base A (2/3);  total (34/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 = 9.99e-01

ACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCC  >  NC_002947/3203010‑3203295
                                                                                                                                            |                                                                                                                                                  
acaTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCTCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGa                                                                                                                                            <  3:157125/149‑1 (MQ=255)
acaTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGa                                                                                                                                            >  4:104023/1‑149 (MQ=255)
acaTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGAAATGATCaactga                                                                                                                                            >  5:48995/1‑144 (MQ=255)
acaTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTAAGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGa                                                                                                                                            >  6:55135/1‑149 (MQ=255)
 aaTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTATAGCGTAATCCGCCTATGACGCATTCTGGCATGTGTCTTGCGTGTTGATCACATGACATGGTCAGTTGa                                                                                                                                            >  2:182118/2‑148 (MQ=255)
     tttCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTTATTTCTCTAGCGTAATCCGCCTATGACTAATTCTGTCATGTGGCTTGAGTGTTGATCAAATGAAATTGTAAGTTGATCATg                                                                                                                                       >  7:278359/1‑149 (MQ=255)
        ccTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCg                                                                                                                                    >  5:250561/1‑149 (MQ=255)
         cTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGt                                                                                                                                               <  5:85968/137‑1 (MQ=255)
         cTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGt                                                                                                                                               >  6:85968/1‑137 (MQ=255)
            aCTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCt                                                                                                                                  <  6:19540/147‑1 (MQ=255)
            aCTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCt                                                                                                                                  >  5:19540/1‑147 (MQ=255)
            aCTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcatg                                                                                                                                  <  6:48995/147‑4 (MQ=255)
            aCTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTAAGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCt                                                                                                                                  <  5:55135/147‑1 (MQ=255)
              ttCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGAt                                                                                                                              >  7:44926/1‑149 (MQ=255)
                 ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                                   >  3:256188/1‑139 (MQ=255)
                 ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                                   <  4:256188/141‑3 (MQ=255)
                  cAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAgg                                                                                                                          <  8:302896/149‑1 (MQ=255)
                       atGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATgta                                                                                                                    <  5:110524/149‑1 (MQ=255)
                                  ccAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGa                                                                                                         <  3:97719/149‑1 (MQ=255)
                                      gCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCtt                                                                                                     <  6:250561/149‑1 (MQ=255)
                                          gCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGAATGTAGTAGGCGCCGATCTTTCCa                                                                                                 <  7:51035/149‑1 (MQ=255)
                                                gTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATgt                                                                                                                     >  5:288609/1‑123 (MQ=255)
                                                gTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATgt                                                                                                                     <  6:288609/123‑1 (MQ=255)
                                                gTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGaa                                                                                           <  8:112782/149‑1 (MQ=255)
                                                 tCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTTAGTAGGCGCCGATCTTTCCACCCGaa                                                                                           <  1:4629/149‑1 (MQ=255)
                                                   gATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTAGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAggg                                                                                        >  1:274756/1‑149 (MQ=255)
                                                     taaataaaTTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTc                                                                                      <  8:44926/149‑1 (MQ=255)
                                                        ataaatTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGag                                                                                   >  7:254335/1‑149 (MQ=255)
                                                            atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGagcagc                                                                               >  3:349419/1‑149 (MQ=255)
                                                            atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGagcagc                                                                               <  6:77435/149‑1 (MQ=255)
                                                                 aCTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGa                                                                                                                               <  4:260406/97‑1 (MQ=255)
                                                                 aCTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGa                                                                                                                               >  3:260406/1‑97 (MQ=255)
                                                                 aCTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTc                                                                          <  8:254335/149‑1 (MQ=255)
                                                                   tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGaa                                                                                           <  5:224244/130‑1 (MQ=255)
                                                                   tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGaa                                                                                           >  6:224244/1‑130 (MQ=255)
                                                                       tGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATgtagta                                                                                                                 >  6:4816/1‑104 (MQ=255)
                                                                       tGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATgtagta                                                                                                                 <  5:4816/104‑1 (MQ=255)
                                                                        gATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTca                                                                                                                                    <  6:288339/85‑2 (MQ=38)
                                                                          tttCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCtt                                                                 >  5:348548/1‑149 (MQ=255)
                                                                           ttCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAg                                                                                          >  4:362286/1‑123 (MQ=255)
                                                                           ttCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAg                                                                                          <  3:362286/123‑1 (MQ=255)
                                                                            tctcTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTc                                                                <  3:122286/148‑1 (MQ=255)
                                                                            tctcTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTc                                                                <  2:274756/148‑1 (MQ=255)
                                                                                tAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATgta                                                                                                                    <  8:58365/92‑1 (MQ=255)
                                                                                tAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATgta                                                                                                                    >  7:58365/1‑92 (MQ=255)
                                                                                tAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGtt                                                             <  6:87450/147‑1 (MQ=255)
                                                                                tAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTg                                                           >  8:71320/1‑149 (MQ=255)
                                                                                tAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTg                                                           >  6:271977/1‑149 (MQ=255)
                                                                                 aGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGa                                                          >  8:245543/1‑149 (MQ=255)
                                                                                  gCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGa                                                          <  8:299868/148‑1 (MQ=255)
                                                                                        tCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGAtcct                                                   >  3:148143/1‑149 (MQ=255)
                                                                                          cGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGAtcctc                                                  >  3:131168/1‑148 (MQ=255)
                                                                                           gCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAg                                                         >  1:193183/1‑140 (MQ=255)
                                                                                           gCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAg                                                         <  2:193183/140‑1 (MQ=255)
                                                                                                    gCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACaa                                       <  7:245543/149‑1 (MQ=255)
                                                                                                       ttCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGa                                                          <  2:34040/127‑1 (MQ=255)
                                                                                                             cATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGttt                                                            <  1:260363/119‑1 (MQ=255)
                                                                                                             cATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGttt                                                            >  2:260363/1‑119 (MQ=255)
                                                                                                             cATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCa                              <  3:104023/149‑1 (MQ=255)
                                                                                                                 tgGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGt                          <  6:347342/149‑1 (MQ=255)
                                                                                                                 tgGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGt                          <  6:341918/149‑1 (MQ=255)
                                                                                                                   gCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCg                        >  5:2230/1‑149 (MQ=255)
                                                                                                                   gCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCg                        >  5:131465/1‑149 (MQ=255)
                                                                                                                     ttGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCgg                       >  2:325433/1‑148 (MQ=255)
                                                                                                                       gCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGa                      <  7:10063/147‑1 (MQ=255)
                                                                                                                                tCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTAt           >  7:210280/1‑149 (MQ=255)
                                                                                                                                   catgacatgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGACGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGAAAGGTCTTAAAGTGTctgatcag        >  8:320413/1‑142 (MQ=255)
                                                                                                                                   catgacatgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCCCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTACAAGCCCTTCGATTGAGGATCCTCCAAt                                              >  3:45013/1‑111 (MQ=255)
                                                                                                                                       acatgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATccgccc     <  1:325433/148‑1 (MQ=255)
                                                                                                                                       aaattgTCAGTTGATAATTTATATTTATCAGTATGTAGTAGGATCATATTTTTCTACACGAAGTGTCGAGCAGCAGGTCCAAGACCTTCGTTTGAGGATCCTACAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATccgccc     <  8:278359/143‑1 (MQ=255)
                                                                                                                                         atgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTc               >  7:283822/1‑136 (MQ=255)
                                                                                                                                         atgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTc               <  8:283822/136‑1 (MQ=255)
                                                                                                                                         atgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCcgcccgcc  <  4:131168/149‑1 (MQ=255)
                                                                                                                                           ggTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGa                                                          <  6:335267/91‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
ACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCC  >  NC_002947/3203010‑3203295

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: