Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,657,318 C→G 6.8% intergenic (+33/‑47) nuoB → / → nuoC NADH‑quinone oxidoreductase subunit B/NADH‑quinone oxidoreductase subunit C/D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,657,3180CG 6.8% 232.4 / 6.5 73intergenic (+33/‑47)nuoB/nuoCNADH‑quinone oxidoreductase subunit B/NADH‑quinone oxidoreductase subunit C/D
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base G (2/3);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

CCTGGGTTGTTGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGATGCATTCGTCGC  >  NC_002947/4657173‑4657465
                                                                                                                                                    |                                                                                                                                                    
ccTGGGTTGTTGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCcgg                                                                                                                                                    <  8:265109/148‑1 (MQ=255)
ccTGGGTTGTTGGTGATCAAGGTATTTACAGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCcggc                                                                                                                                                   >  4:38505/1‑149 (MQ=255)
    ggTTGTTGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCGGCTCCTTCGGGAGCGggcggca                                                                                                                                               <  8:86653/149‑2 (MQ=255)
    ggTTGTTGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCGGCTCCTTCGGGAGCGggcggca                                                                                                                                               <  6:86327/149‑2 (MQ=255)
    ggTTGTTGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCGGCTCCTTCGGGAGCGggcggc                                                                                                                                                <  3:215637/148‑1 (MQ=255)
     gttgttGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCt                                                                                                                                              >  4:177660/1‑149 (MQ=255)
     gttgttGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCt                                                                                                                                              <  2:275975/149‑1 (MQ=255)
           ggTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCtt                                                                                                                                         >  7:274707/1‑148 (MQ=255)
              taTCAAGGGCTCTACCCTGCCGAGATGCAAGCCCAGAAAGACCGGCGTCGGGAGCTGCGCAGTGGAGTTACCAACCTGCGCAGGCCCGACGAAGTCTGATCCACC‑‑GACCTGCCGCCAGCTCCCGCTGGAGCCGGCGGCCTGGCTTCAtt                                                                                                                                     <  8:189190/148‑1 (MQ=255)
               aTCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCAtt                                                                                                                                     <  5:35804/148‑1 (MQ=255)
                 cAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCt                                                                                                                                   >  4:136082/1‑148 (MQ=255)
                  aaGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCGGCTCCTTCGGGAgcgggcgg                                                                                                                                                 >  1:269361/1‑133 (MQ=255)
                  aaGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCGGCTCCTTCGGGAgcgggcgg                                                                                                                                                 <  2:269361/133‑1 (MQ=255)
                   aGGTATTTACCGTGCCGAGATGCCAGCCCAGAAATACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCGTaa                                                                                                                                >  7:157158/1‑149 (MQ=255)
                   aGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTaa                                                                                                                                >  7:162123/1‑149 (MQ=255)
                   aGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTaa                                                                                                                                <  3:268341/149‑1 (MQ=255)
                       aTTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGtt                                                                                                                            >  4:321990/1‑149 (MQ=255)
                          tACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGa                                                                                                                          <  8:7609/148‑1 (MQ=255)
                           aCCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGAcc                                                                                                                        >  7:270525/1‑149 (MQ=255)
                               tGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCaaa                                                                                                                    <  8:270525/149‑1 (MQ=255)
                                   gagaTGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAAc                                                                                                                               >  7:73584/1‑134 (MQ=255)
                                   gagaTGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAAc                                                                                                                               <  8:73584/134‑1 (MQ=255)
                                                          gCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAAc                                                                                         >  1:95410/1‑149 (MQ=255)
                                                          gCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAAc                                                                                         >  1:90562/1‑149 (MQ=255)
                                                                tGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGCGAGACATGCCGCCCGCTCCCGAAGGAGCCGGCGGCATGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGCCCATGACAGCGGACAACGGTa                                                                                     >  1:287641/1‑149 (MQ=255)
                                                                  agcagcGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTAtttt                                                                                 <  3:38505/149‑1 (MQ=255)
                                                                    cagcGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCa                                                                               <  3:172170/149‑1 (MQ=255)
                                                                          aTTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGc                                                                         >  7:158398/1‑149 (MQ=255)
                                                                               aGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTa                                                                     >  7:329376/1‑148 (MQ=255)
                                                                                gTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTAc                                                                    >  5:164872/1‑148 (MQ=255)
                                                                                 tAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGc                                                                         >  6:219259/1‑142 (MQ=255)
                                                                                 tAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGc                                                                         <  5:219259/142‑1 (MQ=255)
                                                                                 tAACCAACCTGCGCAGACCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACaa                                                                  <  8:157158/149‑1 (MQ=255)
                                                                                  aaccaaccTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACaa                                                                  <  8:162123/148‑1 (MQ=255)
                                                                                   accaaccTGCGCAGCCCCGACGAAGTCGTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCATGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGc                                                                         >  1:205071/1‑141 (MQ=255)
                                                                                   accaaccTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGc                                                                         <  2:205071/140‑1 (MQ=255)
                                                                                   accaaccTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAgg                                                                >  2:163142/1‑149 (MQ=255)
                                                                                   accaaccTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAgg                                                                >  1:77702/1‑149 (MQ=255)
                                                                                    ccaaccTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGAcc                                                                                                              <  3:391/102‑1 (MQ=255)
                                                                                    ccaaccTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGAcc                                                                                                              >  4:391/1‑102 (MQ=255)
                                                                                    ccaaccTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGc                                                               <  2:77702/149‑1 (MQ=255)
                                                                                                 ccccGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGAt                                                  >  6:78785/1‑149 (MQ=255)
                                                                                                    cgacgaAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGAtgtg                                               <  3:136082/149‑1 (MQ=255)
                                                                                                    cgacgaAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGAtgtg                                               >  5:155421/1‑149 (MQ=255)
                                                                                                    cgacgaAGTCTGATCCAGC‑‑GACCTGACGCGCGGTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTATCGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACTAGGCCGACGACCacatcgtt                                               <  2:287641/149‑8 (MQ=255)
                                                                                                     gacgaAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTATTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTACGCCCTACAAGGCTGACGACCAGGATGTgg                                              >  6:53684/1‑149 (MQ=255)
                                                                                                       cgaAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGgtc                                            <  3:321990/149‑1 (MQ=255)
                                                                                                          aGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGgtcgtc                                         <  8:158398/149‑1 (MQ=255)
                                                                                                           gTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGtcgtcg                                        <  6:289575/149‑1 (MQ=255)
                                                                                                               tGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAAc                                     >  1:233065/1‑149 (MQ=255)
                                                                                                                gATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAAc                                     <  4:71517/148‑1 (MQ=255)
                                                                                                                gATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAAc                                     >  3:71517/1‑148 (MQ=255)
                                                                                                                 aTCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACt                                    >  3:214092/1‑148 (MQ=255)
                                                                                                                 aTCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACt                                    >  1:220297/1‑148 (MQ=255)
                                                                                                                  tCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGa                                  <  5:78785/149‑1 (MQ=255)
                                                                                                                    cAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGgtc                                            <  7:340812/137‑1 (MQ=255)
                                                                                                                    cAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGgtc                                            >  8:340812/1‑137 (MQ=255)
                                                                                                                             aggccgccGGCTCCTTCGGGAGCGGGCGGCCTGGCTTCATTCTTACCGTTGCCCCAACGCGAGCGCGACCATGACAGCGGACAGCGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGAGGTGGTCGTCGAACTGAACAACCGtt                         >  1:19518/3‑149 (MQ=255)
                                                                                                                             aggccgccGGCTCCTTCGGGAGCGGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGtt                         >  6:65881/3‑149 (MQ=255)
                                                                                                                               gccgccCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGtttt                       >  7:22172/1‑149 (MQ=255)
                                                                                                                                ccgcccgcTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGtttt                       <  8:329376/148‑1 (MQ=255)
                                                                                                                                ccgcccgcTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTg                      >  1:354801/1‑149 (MQ=255)
                                                                                                                                atcccatcTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGTCCCAAACCGTCCGACACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTTAACAACCGtttt                       <  2:330556/141‑1 (MQ=255)
                                                                                                                                    ccgcTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGAcc                                                                                                                        >  2:327246/1‑47 (MQ=255)
                                                                                                                                    ccgcTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGAcc                                                                                                                        <  1:327246/47‑1 (MQ=255)
                                                                                                                                        tCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGAt             >  3:47742/1‑149 (MQ=255)
                                                                                                                                          ccGAAGGAGCCGGCGTCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGtttt                       <  2:330555/138‑1 (MQ=255)
                                                                                                                                          ccGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGAt             <  3:177660/147‑1 (MQ=255)
                                                                                                                                          ccGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGAt             <  8:274707/147‑1 (MQ=255)
                                                                                                                                           cGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGATGCa          >  5:116751/1‑149 (MQ=255)
                                                                                                                                           cGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGATGCa          >  1:50055/1‑149 (MQ=255)
                                                                                                                                            gAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGcgacgaatgca         >  5:122331/1‑144 (MQ=255)
                                                                                                                                            gAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGATGCAt         >  7:167266/1‑149 (MQ=255)
                                                                                                                                            gAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGATGCAt         >  6:340511/1‑149 (MQ=255)
                                                                                                                                             aaGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGgcgc                  <  3:184645/140‑1 (MQ=255)
                                                                                                                                             aaGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGgcgc                  >  4:184645/1‑140 (MQ=255)
                                                                                                                                                gAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGATGCATtcgt     >  4:232332/1‑149 (MQ=255)
                                                                                                                                                    cggcggCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGACGAATGCAtcggcgcc  <  6:122331/149‑8 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTGGGTTGTTGGTGATCAAGGTATTTACCGTGCCGAGATGCCAGCCCAGAAAGACCTGCGTCGTGAGCAGCGCATTGCAGTAACCAACCTGCGCAGCCCCGACGAAGTCTGATCCAGC‑‑GACCTGCCGCCCGCTCCCGAAGGAGCCGGCGGCCTGGCTTCATTCTTAACGTTGACCCAAAGCGACCGAGACCATGACAGCGGACAACGCTATTTTCATTCCGCCCTACAAGGCTGACGACCAGGATGTGGTCGTCGAACTGAACAACCGTTTTGGCGCCGATGCATTCGTCGC  >  NC_002947/4657173‑4657465

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: