Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 133,400 C→G 6.2% L60V (CTG→GTG)  PP_0128 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947133,4000CG 6.2% 261.7 / 6.0 81L60V (CTG→GTG) PP_0128hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base G (3/2);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCACTTGGGG  >  NC_002947/133255‑133546
                                                                                                                                                  |                                                                                                                                                  
ccTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACctgc                                                                                                                                                 <  6:169403/149‑1 (MQ=255)
   gCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCtgctgc                                                                                                                                              >  2:73193/1‑149 (MQ=255)
    cACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGcc                                                                                                                                             >  7:91964/1‑149 (MQ=255)
     aCGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTAACGCCATTACCTGACCCGTAGCTGGCAGCACCTGCTGCCg                                                                                                                                            >  5:205096/1‑149 (MQ=255)
               ggTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  >  7:130863/1‑149 (MQ=255)
                gTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCg                                                                                                                                            <  1:335373/138‑1 (MQ=255)
                gTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCGGCTGCCg                                                                                                                                            >  2:335373/1‑138 (MQ=255)
                 tCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGt                                                                                                                                <  3:362752/149‑1 (MQ=255)
                 tCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGt                                                                                                                                <  2:191850/149‑1 (MQ=255)
                       aGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACCGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCgg                                                                                                                           <  1:73193/149‑1 (MQ=255)
                         cATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCAGGTGCTGcc                                                                                                                                             >  6:176252/1‑128 (MQ=255)
                         cATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCAGGTGCTGcc                                                                                                                                             <  5:176252/128‑1 (MQ=255)
                           tCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAAcc                                                                                                                      >  8:21340/1‑149 (MQ=255)
                           tCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAAcc                                                                                                                      >  8:26883/1‑149 (MQ=255)
                               gggCAGGCCCGGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGCGCGTGCCGGTAACCTGc                                                                                                                   >  8:340235/3‑148 (MQ=255)
                               cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  <  1:356396/133‑1 (MQ=255)
                               cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  >  2:350630/1‑133 (MQ=255)
                               cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  >  2:356396/1‑133 (MQ=255)
                               cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  <  1:350630/133‑1 (MQ=255)
                                  cAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCag                                                                                                                 >  6:7352/1‑147 (MQ=255)
                                                  ccGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGAGGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGACGGTAACCTGCAGAAGATCGGCAGGctgctg                                                                                               >  7:246461/1‑149 (MQ=255)
                                                       gATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGa                                                                                          >  6:94994/1‑149 (MQ=255)
                                                                   cGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCt                                                                                                                     <  5:119862/110‑1 (MQ=255)
                                                                   cGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtg                                                                                <  8:91964/147‑1 (MQ=255)
                                                                   cGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCCCTACCGCCATTACCTGACCCCTAGCTGGCAGCACCTTCTGCCGCATAACGGGCGTGCCGGTAAcca                                                                                                                     >  6:119862/1‑109 (MQ=255)
                                                                         ctCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGccc                                                                         <  6:205096/148‑1 (MQ=255)
                                                                         ctCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCt                                                                        >  3:141644/1‑149 (MQ=255)
                                                                             gTTTCAACATCCAGGTGGGCATCATCACCGAGCGCTACCGCCATTACCTTACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAg                                                                    <  7:340235/149‑1 (MQ=255)
                                                                             gTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTc                                                                                      >  3:108626/1‑131 (MQ=255)
                                                                             gTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTc                                                                                      <  4:108626/131‑1 (MQ=255)
                                                                              tttCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAg                                                                    >  6:244386/1‑148 (MQ=255)
                                                                               ttCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   >  7:227882/1‑148 (MQ=255)
                                                                                tCAACCTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   >  7:346077/1‑147 (MQ=255)
                                                                                 cAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAg                                                                >  6:135102/1‑149 (MQ=255)
                                                                                  aaCATCCCGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACTCGTCGCTGTCAGCACCTGCTGCCGCATAACGTGCGTGCCGGTAACCTGGAGATGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGc                                                               <  6:200529/149‑1 (MQ=255)
                                                                                  aaCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTg                                                                       <  8:157160/141‑1 (MQ=255)
                                                                                  aaCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTg                                                                       >  7:157160/1‑141 (MQ=255)
                                                                                  aaCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGc                                                               <  8:130863/149‑1 (MQ=255)
                                                                                   aCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   >  5:266648/1‑144 (MQ=255)
                                                                                   aCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   <  6:266648/144‑1 (MQ=255)
                                                                                     aTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGa                                                            <  7:206473/149‑1 (MQ=255)
                                                                                     aTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGa                                                            <  5:94994/149‑1 (MQ=255)
                                                                                     aTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGa                                                            >  5:69132/1‑149 (MQ=255)
                                                                                                 aTCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCt                                                >  8:336433/1‑149 (MQ=255)
                                                                                                        ccGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCg                                         >  4:356368/1‑149 (MQ=255)
                                                                                                        ccGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCg                                         <  5:244386/149‑1 (MQ=255)
                                                                                                        ccGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCg                                         <  4:141644/149‑1 (MQ=255)
                                                                                                         cGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGATGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCgg                                        >  7:46157/1‑149 (MQ=255)
                                                                                                          gAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCggg                                       <  3:356368/149‑1 (MQ=255)
                                                                                                              gcTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTAc                                   <  7:26883/149‑1 (MQ=255)
                                                                                                              gcTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTAc                                   <  7:21340/149‑1 (MQ=255)
                                                                                                                 aCCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            >  6:131027/1‑105 (MQ=255)
                                                                                                                 aCCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            <  5:131027/105‑1 (MQ=255)
                                                                                                                   cGCCTTAACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCAGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCaa                              >  2:112531/1‑149 (MQ=255)
                                                                                                                   cGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCa                                                    >  8:159496/1‑127 (MQ=255)
                                                                                                                   cGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCa                                                    <  7:159496/127‑1 (MQ=255)
                                                                                                                    gCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAAc                             <  5:7352/149‑1 (MQ=255)
                                                                                                                     ccATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGAc                                                                                   >  1:319644/1‑94 (MQ=255)
                                                                                                                     ccATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGAc                                                                                   <  2:319644/94‑1 (MQ=255)
                                                                                                                     ccATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAAc                             >  5:8794/1‑148 (MQ=255)
                                                                                                                      cATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCa                           <  8:346077/149‑1 (MQ=255)
                                                                                                                       aTTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAg                          <  8:246461/149‑1 (MQ=255)
                                                                                                                       aTTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAg                          <  5:135102/149‑1 (MQ=255)
                                                                                                                       aTTACCTGACGCGTAGCTGGCAGCACCTGCTGCCACATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAg                          >  5:180990/1‑149 (MQ=255)
                                                                                                                        ttACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAgg                         >  3:289061/1‑149 (MQ=255)
                                                                                                                          aCCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGGCCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGt                        <  8:39988/148‑1 (MQ=255)
                                                                                                                            cTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACgg                                                                                                                                    >  4:183335/1‑38 (MQ=255)
                                                                                                                            cTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACgg                                                                                                                                    <  3:183335/38‑1 (MQ=255)
                                                                                                                            cTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCg                      >  3:333586/1‑148 (MQ=255)
                                                                                                                             tGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGCAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGaa                    >  6:79508/1‑149 (MQ=255)
                                                                                                                                  gcgTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            <  1:254476/89‑1 (MQ=255)
                                                                                                                                  gcgTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            >  2:254476/1‑89 (MQ=255)
                                                                                                                                     tAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTgg             <  1:362684/149‑1 (MQ=255)
                                                                                                                                     tAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTgg             <  8:46157/149‑1 (MQ=255)
                                                                                                                                       gCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTggtt           >  8:107406/1‑147 (MQ=255)
                                                                                                                                        gcggCAGCAGGTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGc                                                                                                        <  2:147938/53‑1 (MQ=255)
                                                                                                                                        gcggCAGCAGGTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGc                                                                                                        >  1:147938/3‑55 (MQ=38)
                                                                                                                                        cTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTggttg          <  7:107406/149‑4 (MQ=255)
                                                                                                                                         cggCAGCAGGTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCa         >  2:38400/2‑149 (MQ=255)
                                                                                                                                             agcaCCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGcc           <  8:108129/143‑1 (MQ=255)
                                                                                                                                             agcaCCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGcc           >  7:108129/1‑143 (MQ=255)
                                                                                                                                                aCCTGCTGCCACATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCACTTgggg  <  6:180990/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCACTTGGGG  >  NC_002947/133255‑133546

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: