Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 161,731 G→T 9.2% intergenic (‑226/+20) PP_0152 ← / ← PP_0153 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947161,7310GT 9.2% 178.2 / 3.0 71intergenic (‑226/+20)PP_0152/PP_0153hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (40/22);  new base T (2/6);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 5.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTGCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCACAGCAGCAG  >  NC_002947/161586‑161879
                                                                                                                                                 |                                                                                                                                                       
ctgTGCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGaa                                                                                                                                                       <  4:341876/148‑1 (MQ=255)
    gCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCgg                                                                                                                                               >  7:30956/1‑149 (MQ=255)
     ccAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCgg                                                                                                                                               <  3:306834/148‑1 (MQ=255)
      cAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAATTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGc                                                                                                                                              >  1:180850/1‑148 (MQ=255)
       aGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGc                                                                                                                                              >  3:214029/1‑147 (MQ=255)
       aGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGc                                                                                                                                              >  6:334720/1‑147 (MQ=255)
         aCAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGc                                                                                                                                          >  6:356231/1‑149 (MQ=255)
          cAACTGTAACGGCGAATTTTAATTTTTGTAAACCTGCCCGGGCCCCAACCTCACCCCCGCGACGGCGCGTCTCTCGCGCAAAAAAAACCCCGCGCCCGTCACCTCCACCGGGTTTGAGAAATTCCCGGTGCCGGCGAA‑‑‑GCCGGCAAGcc                                                                                                                                         <  2:188486/149‑1 (MQ=255)
           aaCTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCg                                                                                                                                        >  8:132817/1‑149 (MQ=255)
              tGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCgg                                                                                                                                       >  7:67281/1‑147 (MQ=255)
                   aGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCt                                                                                                                                  >  2:264674/1‑147 (MQ=255)
                          ttttaattttTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAgcgc                                                                                                                            >  5:322374/1‑146 (MQ=255)
                          ttttaattttTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAgcgc                                                                                                                            >  3:263215/1‑146 (MQ=255)
                             taattttTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCGGTCACCTGCAACGGGTTTGAGAAATACACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGa                                                                                                                        >  2:159011/1‑147 (MQ=255)
                              attttttGTAAAACTGTACGGGTCGAAATCTCAAGCACGCTACTGCGCGTCTCTCGGCCAATAAAAACCCCGCTGCAGTCACCTGCAAGGGTTTTTATAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGaaa                                                                                                                      <  4:230706/146‑1 (MQ=255)
                                    tGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCTGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTg                                                                                                                <  4:97077/148‑1 (MQ=255)
                                     gTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCt                                                                                                              >  2:135512/1‑149 (MQ=255)
                                      tAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCTTC‑‑‑GCCGGCaccgt                                                                                                                                         >  5:21785/1‑117 (MQ=255)
                                      tAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCTTC‑‑‑GCCGGCaccgt                                                                                                                                         <  6:21785/121‑5 (MQ=255)
                                       aaacCTGAACTGGTCTCAAGCTTAAGCCGGCGACGGCGCGTCTCTCGGGAAATAAAAAACCCGCTGAAGTCACCTGAATCGGGTTTTATAAATTCCCGGTGCCGGCTAC‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTc                                                                                                             <  8:69674/148‑1 (MQ=255)
                                       aaaaCTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTc                                                                                                             >  7:69674/1‑148 (MQ=255)
                                        aaaCTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAg                                                                                                           <  8:67281/149‑1 (MQ=255)
                                                   gTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAAGGTGTAGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCg                                                                                                    >  1:124329/1‑148 (MQ=255)
                                                   gTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGa                                                                                                                                     >  1:103345/1‑112 (MQ=255)
                                                    tCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACg                                                                                               >  1:361343/1‑149 (MQ=255)
                                                      gaatccTGAACCCCGCGTCTGGGCGTCTGTCGGGAAATAAAACCCCCTGTGCATTCACTTTCACGGTGTTTTATAATTTCACGTGGCGGGGCAC‑‑‑GCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCt                                                                                             <  2:191593/144‑1 (MQ=255)
                                                       cAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTg                                                                                            >  5:377964/1‑149 (MQ=255)
                                                              aaGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGGAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCg                                                                                     >  1:281075/1‑149 (MQ=255)
                                                              aaGCACCCCACTCCCCGTCTATCGGGAAATAAAAAACCCGCTCACGTCCACTTCAAGGGTTTTTATAAATAAACGGTCCCGGCTAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTc                                                                                      <  4:263215/148‑1 (MQ=255)
                                                                  aCGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGa                                                                                 >  5:16435/1‑149 (MQ=255)
                                                                   cgcgACTGCGCGTCTCTCGGGCAATAAAAACCCCGCTGCAGTCACCTGCAAGGGGTTTTATAAATTCACGGGGCCGGCGAC‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGAt                                                                                <  1:264674/149‑1 (MQ=255)
                                                                    gctaCTGCGCTTCTCTCGGGAAAAAAAAAACCCGCTCCAGACACCTGCAAGGGTTTTTATAAATTCACGGTGCCGGCGAC‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATa                                                                               <  7:132817/146‑1 (MQ=255)
                                                                    gcgACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATa                                                                               >  3:107861/1‑149 (MQ=255)
                                                                            gcgTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTc                                                                       >  3:337315/1‑149 (MQ=255)
                                                                                 tctcGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGc                                                                    >  5:69516/1‑147 (MQ=255)
                                                                                     gggCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCtt                                                                                                                 >  3:270041/1‑98 (MQ=255)
                                                                                     gggCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCtt                                                                                                                 <  4:270041/98‑1 (MQ=255)
                                                                                      ggCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATAACGGCGTCTGCCAGCAgg                                                             >  3:192484/1‑149 (MQ=255)
                                                                                        cAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCg                                                           >  4:4884/1‑149 (MQ=255)
                                                                                         aaaaaaacaCCCGCTGCCGCCCCCTCCTCCGGTTTTTATACATCCCCGGCGCGGGCGAT‑‑‑CCCGCCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                          <  4:192484/146‑1 (MQ=255)
                                                                                         aaaaaaaacCCCGCTGCAGTCACCTGCAACGGGTTTGATAAATTCACGGGGCCGGCGAC‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                          <  2:124329/146‑1 (MQ=255)
                                                                                         aaTAAAAACCCCCCGGCCGTCACCTCCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                          <  8:185038/149‑1 (MQ=255)
                                                                                         aaTAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTTATAAATTCACGTGGCGGGCTAA‑‑‑GCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                          <  8:30956/149‑1 (MQ=255)
                                                                                         aaTAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                          >  1:330201/1‑149 (MQ=255)
                                                                                           tAAAAAACCCGCTGCAGTCCCCTGCAACGGTTTTGAGAAATGCACGGGGCCGGCGAC‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGa                                                        <  5:356231/149‑1 (MQ=255)
                                                                                           aaaaacccccGCTGCCGTCACCTGCAACGGGTTTGAGAAATTCCCGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGa                                                        <  8:152422/142‑1 (MQ=255)
                                                                                           aaaaaaaCCCGCTGCAGTCACCTCCAAGGGTTTTTATAATTTCCGGGGGCGGGCTAC‑‑‑GCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGa                                                        <  6:16435/148‑1 (MQ=255)
                                                                                           aaaaaaaCCCGCTCCAGCCACCTCAAAGGGGTTTTATAAATTCCCGGTCCGGGCTAC‑‑‑GCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGa                                                        <  5:328438/148‑1 (MQ=255)
                                                                                             aaaaaCCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGa                                                        >  7:152422/1‑147 (MQ=255)
                                                                                              caaaCCCGTTGCAGGTGACTGCAGCGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCAGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTc                                                     >  2:232258/2‑149 (MQ=255)
                                                                                                     gCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTg                                                                                                                >  1:354099/1‑83 (MQ=255)
                                                                                                         cAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGc                                            >  1:333036/1‑147 (MQ=255)
                                                                                                           gTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGt                                        >  1:187946/1‑149 (MQ=255)
                                                                                                               ccTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTgac                                    >  6:369104/1‑149 (MQ=255)
                                                                                                                cTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTgac                                    >  2:23923/1‑148 (MQ=255)
                                                                                                                   cAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTgacgac                                 <  2:281075/148‑1 (MQ=255)
                                                                                                                   cAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGacgacg                                >  6:345847/1‑149 (MQ=255)
                                                                                                                        ggTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTgg                           <  4:107861/149‑1 (MQ=255)
                                                                                                                          tttGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGt                          >  3:217566/1‑148 (MQ=255)
                                                                                                                              ataatttCACGGTGCCGGCGAC‑‑‑GCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATa                      <  2:187946/143‑1 (MQ=255)
                                                                                                                              agaAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCCGGCGGGAGTCCAGGCGAGCCAGTTGACGCCGGCTGGTCAtcc                     >  6:117484/1‑147 (MQ=255)
                                                                                                                              agaAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATAc                     >  5:279139/1‑149 (MQ=255)
                                                                                                                              agaAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATAc                     >  2:34783/1‑149 (MQ=255)
                                                                                                                                 aaTTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACgg                   <  6:69516/148‑1 (MQ=255)
                                                                                                                                      aCGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTg              <  6:377964/148‑1 (MQ=255)
                                                                                                                                         ttGCCGGCTTC‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCca           >  1:138895/2‑148 (MQ=255)
                                                                                                                                          tGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCcaca         <  5:334720/149‑1 (MQ=255)
                                                                                                                                            ccGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCAcagc       <  6:279139/149‑1 (MQ=255)
                                                                                                                                              ggCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGgc                  >  1:23800/1‑136 (MQ=255)
                                                                                                                                              ggCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGgc                  <  2:23800/136‑1 (MQ=255)
                                                                                                                                               gCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCAcagcag     >  7:168882/1‑148 (MQ=255)
                                                                                                                                                 gAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCAcagcagcag  <  1:23923/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CTGTGCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCCGCTGCAGTCACCTGCAACGGGTTTGAGAAATTCACGGTGCCGGCGAA‑‑‑GCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCACAGCAGCAG  >  NC_002947/161586‑161879

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: