Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,804,866 A→G 5.2% C68R (TGC→CGC)  PP_4233 ← oxidoreductase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,804,8660AG 5.2% 232.9 / 2.1 77C68R (TGC→CGC) PP_4233oxidoreductase small subunit
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/2);  total (38/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 = 9.50e-01

GTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCG  >  NC_002947/4804719‑4805014
                                                                                                                                                   |                                                                                                                                                    
gTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAg                                                                                                                                                     <  6:42872/149‑1 (MQ=255)
        tAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGTCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCAc                                                                                                                                             <  6:219809/149‑1 (MQ=255)
              cACTGGAAGGCGAACAGGTCGACAAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACttcttc                                                                                                                                       >  8:33141/1‑149 (MQ=255)
                      ggCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGt                                                                                                                                                   >  5:153296/1‑129 (MQ=255)
                      ggCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGt                                                                                                                                                   <  6:153296/129‑1 (MQ=255)
                           aCAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCTTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCg                                                                                                                          >  3:46343/1‑149 (MQ=255)
                              ggTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCgggg                                                                                                                       <  5:398739/149‑1 (MQ=255)
                               gTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGt                                                                                                                      >  4:56559/1‑149 (MQ=255)
                                tCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGtt                                                                                                                     >  7:304055/1‑149 (MQ=255)
                                tCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGtt                                                                                                                     >  3:160496/1‑149 (MQ=255)
                                   acgaACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTc                                                                                                                  <  8:304055/149‑1 (MQ=255)
                                     gaACTTCTGCTAGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTcga                                                                                                                >  1:75308/1‑149 (MQ=255)
                                       aCTTCTGCTGGATCGGCTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTcgacg                                                                                                              <  3:55932/149‑1 (MQ=255)
                                           ctgctgGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGAt                                                                                                            <  8:311477/147‑1 (MQ=255)
                                            tgctgGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTg                                                                                                                                >  4:288758/1‑126 (MQ=255)
                                            tgctgGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTg                                                                                                                                <  3:288758/126‑1 (MQ=255)
                                                 gATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACg                                                                                                     >  1:410874/1‑148 (MQ=255)
                                                   tggggTTGAGGGTCAGGTTGCCGGCCTCGTGTTTCTTGGCATGGCCTTCGATTGTGAGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATTCAGGTg                                                                                                                                                  <  4:342328/99‑1 (MQ=255)
                                                   tCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTg                                                                                                                                                  >  3:342328/1‑101 (MQ=255)
                                                    cGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGAt                                                                                                            >  3:175329/1‑138 (MQ=255)
                                                    cGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGAt                                                                                                            <  4:175329/138‑1 (MQ=255)
                                                      ggTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCg                                                                                               >  4:104844/1‑149 (MQ=255)
                                                       gTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACttcttc                                                                                                                                       <  8:264428/108‑1 (MQ=255)
                                                       gTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACttcttc                                                                                                                                       >  7:264428/1‑108 (MQ=255)
                                                        ttGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTg                                                                                             <  4:46343/149‑1 (MQ=255)
                                                        ttGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGACGTCGGGGTTGTCGACGATCACCACGCAGGCGTg                                                                                             >  5:35794/1‑149 (MQ=255)
                                                            gCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTgg                                                                                            >  5:28073/1‑146 (MQ=255)
                                                               tCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATg                                                                                      >  5:304736/1‑149 (MQ=255)
                                                                      gCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTgg                                                                                <  1:233255/148‑1 (MQ=255)
                                                                      gCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTgg                                                                                <  5:254712/148‑1 (MQ=255)
                                                                      gCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTgg                                                                                <  7:47648/148‑1 (MQ=255)
                                                                       ccGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCGGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGc                                                                               >  8:363268/1‑148 (MQ=255)
                                                                       ccGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCGGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGc                                                                               <  7:363268/148‑1 (MQ=255)
                                                                       ccGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGcc                                                                              <  7:202321/149‑1 (MQ=255)
                                                                       ccGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGc                                                                               <  5:104671/148‑1 (MQ=255)
                                                                       ccGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGc                                                                               <  7:33141/148‑1 (MQ=255)
                                                                       ccGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGc                                                                               <  3:56559/148‑1 (MQ=255)
                                                                               gTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAAcc                                                                        <  3:368143/147‑1 (MQ=255)
                                                                                  tttcttGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAg                                                                     <  8:231049/147‑1 (MQ=255)
                                                                                  tttcttGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAg                                                                     >  7:231049/1‑147 (MQ=255)
                                                                                      ctGGCATGGCCTTCTATGTCGCGCACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGa                                                               <  6:289018/148‑1 (MQ=255)
                                                                                        ggCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGa                                                               <  3:104844/147‑1 (MQ=255)
                                                                                           aTGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCgg                                                          <  6:304736/149‑1 (MQ=255)
                                                                                            tGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGt                                                         <  7:53253/149‑1 (MQ=255)
                                                                                              gCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGt                                                         >  3:66578/1‑147 (MQ=255)
                                                                                                     aTGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTg                                                <  4:160496/149‑1 (MQ=255)
                                                                                                      tGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTgg                                               >  2:371923/1‑149 (MQ=255)
                                                                                                           cGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCGGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTc                                          >  3:278239/1‑149 (MQ=255)
                                                                                                           cGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTc                                          >  2:97714/1‑149 (MQ=255)
                                                                                                             aaCCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGt                                                                                                                                 >  2:336594/1‑60 (MQ=255)
                                                                                                             aaCCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGt                                                                                                                                 >  2:332069/1‑60 (MQ=255)
                                                                                                             aaCCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTg                                                                                             >  4:10886/1‑96 (MQ=255)
                                                                                                             aaCCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTg                                                                                             <  3:10886/96‑1 (MQ=255)
                                                                                                              aCCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGt                                        <  2:75308/148‑1 (MQ=255)
                                                                                                               ccTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTg                                       >  1:253745/1‑148 (MQ=255)
                                                                                                                cttcttGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTg                                       >  4:262704/1‑147 (MQ=255)
                                                                                                                   cttGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAgg                                   <  6:28073/148‑1 (MQ=255)
                                                                                                                   cttGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAg                                    <  6:194083/147‑1 (MQ=255)
                                                                                                                      gCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATg                                                                                      >  4:311099/1‑94 (MQ=255)
                                                                                                                      gCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATg                                                                                      <  3:311099/94‑1 (MQ=255)
                                                                                                                       ccTTCGAAGTAATGCGCACCGGTGATGCGGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGt                               <  4:278239/148‑1 (MQ=255)
                                                                                                                       ccTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTc                              >  1:382124/1‑149 (MQ=255)
                                                                                                                       ccTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTc                              >  1:388287/1‑149 (MQ=255)
                                                                                                                        cTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCAACACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTc                              >  3:32986/1‑148 (MQ=255)
                                                                                                                         ttCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGa                            >  3:309635/1‑149 (MQ=255)
                                                                                                                             aaGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCTGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGatcat                        >  2:108088/1‑149 (MQ=255)
                                                                                                                              aGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGatcatc                       <  4:66578/149‑1 (MQ=255)
                                                                                                                                 aaTGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGATCATCGcc                    >  8:141204/1‑149 (MQ=255)
                                                                                                                                   tGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGATCATCGCCAg                  <  2:215435/149‑1 (MQ=255)
                                                                                                                                      gcACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTgg                                                                                <  6:347840/84‑1 (MQ=255)
                                                                                                                                      gcACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTg                                                <  8:62667/116‑1 (MQ=255)
                                                                                                                                      gcACCGGTGATGCAGGTGCGCAATTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTg                                                >  7:62667/1‑116 (MQ=255)
                                                                                                                                        aCCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGt                                        <  2:378620/122‑1 (MQ=255)
                                                                                                                                        aCCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGATCATCGCCAGGTCcg             >  4:70359/1‑149 (MQ=255)
                                                                                                                                            gTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGATCATCGCCAGGTCCGCgg          >  2:75098/1‑148 (MQ=255)
                                                                                                                                                 gCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAgg                                   <  3:94790/118‑1 (MQ=255)
                                                                                                                                                 gCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAgg                                   >  4:94790/1‑118 (MQ=255)
                                                                                                                                                   aGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCg  <  3:70359/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTGCGCAGTAGCTGCACTGGAAGGCGAACAGGTCGACGAACTTCTGCTGGATCGGGTTGAGCGTCAGGTTGCCGGCCTCGTCTTTCTTGGCATGGCCTTCGATGGTGCGAACCTTCTTGCCTTCGAAGTAATGCGCACCGGTGATGCAGGTGCGCACTTCTTCGCTGGTGCCGTCGGGGTTGTCGACGATCACCACGCAGGCGTGGCAGATGCCCTGGCCGCAACCCAGGCGCGAGCCGGTGAGGTTCTGGTGTTCGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCG  >  NC_002947/4804719‑4805014

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: