Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 232,405 T→C 5.5% P168P (CCT→CCC)  PP_0180 → cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947232,4050TC 5.5% 212.5 / 2.1 72P168P (CCT→CCC) PP_0180cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base T (32/36);  new base C (2/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

CGCGACTGGCGGTGGCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGCTGGACGACCGCCAACGCTGGGATGTGG  >  NC_002947/232257‑232547
                                                                                                                                                    |                                                                                                                                              
cgcgACTGGCGGTGGCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAAGTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCt                                                                                                                                                <  2:288492/149‑1 (MQ=255)
 gcgACTGGCGGTGGCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTg                                                                                                                                               >  1:180486/1‑149 (MQ=255)
   gACTGGCGGTGGCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCa                                                                                                                                             >  1:361768/1‑149 (MQ=255)
       ggCGGTGGCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCg                                                                                                                                         >  1:178253/1‑149 (MQ=255)
         cGGTGGCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGcc                                                                                                                                       <  2:180486/149‑1 (MQ=255)
             ggCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAAcc                                                                                                                                   >  3:298946/1‑149 (MQ=255)
                cTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTgc                                                                                                                                <  7:120625/149‑1 (MQ=255)
                   cGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGcgcg                                                                                                                             <  6:69904/149‑1 (MQ=255)
                      gaTCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATg                                                                                                                          >  4:110877/1‑149 (MQ=255)
                          aGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGtctc                                                                                                                      <  4:298946/149‑1 (MQ=255)
                          aGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGtctc                                                                                                                      >  7:12097/1‑149 (MQ=255)
                             cAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTcgcg                                                                                                                   >  2:137538/1‑149 (MQ=255)
                              aGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCgcgc                                                                                                                  >  8:181073/1‑149 (MQ=255)
                                    cGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGAc                                                                                                            >  4:329065/1‑149 (MQ=255)
                                    cGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGAc                                                                                                            >  2:293101/1‑149 (MQ=255)
                                       tGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAg                                                                                                         >  1:246898/1‑149 (MQ=255)
                                         aTTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGt                                                                                                        <  3:110877/148‑1 (MQ=255)
                                            aCCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTaa                                                                                                     <  5:983/148‑1 (MQ=255)
                                            aCCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTa                                                                                                      <  6:191027/147‑1 (MQ=255)
                                            aCCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTa                                                                                                      >  5:191027/1‑147 (MQ=255)
                                             ccccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCTCGGCGACACCGGCGCCGGCGATGGCCGGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGc                                                                                                   <  8:12097/149‑1 (MQ=255)
                                             ccccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTaa                                                                                                     >  3:348557/1‑147 (MQ=255)
                                             ccccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTaa                                                                                                     <  4:348557/147‑1 (MQ=255)
                                             ccccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGc                                                                                                   >  2:308320/1‑149 (MQ=255)
                                              cccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTg                                                                                                                                               <  6:73845/104‑1 (MQ=255)
                                              cccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTgcgc                                                                                                                              >  1:332991/1‑121 (MQ=255)
                                              cccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTgcgc                                                                                                                              <  2:332991/121‑1 (MQ=255)
                                              cccGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGcc                                                                                                  >  6:176123/1‑149 (MQ=255)
                                                cGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTg                                                                                                >  5:108232/1‑149 (MQ=255)
                                                  gACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGtt                                                                                              <  1:293101/149‑1 (MQ=255)
                                                       cgcccgcGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCTAcc                                                                                         <  3:265053/149‑1 (MQ=255)
                                                       cgcccgcGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGAcc                                                                                         >  6:143012/1‑149 (MQ=255)
                                                       cgcccgcGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCAAcc                                                                                         >  4:265053/1‑149 (MQ=255)
                                                       cgcccgcGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCCGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGAcc                                                                                         >  3:311947/1‑149 (MQ=255)
                                                            gcggcgCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCtat                                                                                    <  5:176123/149‑1 (MQ=255)
                                                            gcggcgCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCtat                                                                                    <  6:221452/149‑1 (MQ=255)
                                                                  ctgcgctgTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAcc                                                                              <  2:361768/149‑1 (MQ=255)
                                                                     cgctgTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccc                                                                             <  3:329065/147‑1 (MQ=255)
                                                                            ctcacaGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccctggccctgg                                                                    <  5:10343/147‑1 (MQ=255)
                                                                            cGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccctggccctgg                                                                    <  2:178253/149‑1 (MQ=255)
                                                                            cGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccctggccctgg                                                                    <  4:308584/149‑1 (MQ=255)
                                                                               acaGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCACCCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccctggccctgg                                                                    >  6:10343/1‑146 (MQ=255)
                                                                               acaGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGta                                                                 >  1:142574/1‑149 (MQ=255)
                                                                               acaGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGta                                                                 <  1:308320/149‑1 (MQ=255)
                                                                               acaGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGta                                                                 <  2:246898/149‑1 (MQ=255)
                                                                                 aGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGta                                                                 >  3:272453/1‑147 (MQ=255)
                                                                                   aaCTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccctgg                                                                          <  4:61514/136‑1 (MQ=255)
                                                                                   aaCTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccctgg                                                                          >  3:61514/1‑136 (MQ=255)
                                                                                    aCTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGAc                                                            <  4:57507/149‑1 (MQ=255)
                                                                                       gCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAac                                                                                 <  8:17400/125‑1 (MQ=255)
                                                                                       gCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAac                                                                                 >  7:17400/1‑125 (MQ=255)
                                                                                           gATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCAcc                                                      >  5:212456/1‑148 (MQ=255)
                                                                                                gCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGTGTGGGCCTGGAGCCCGCACCCGCCAACCTGCGCGAt                                                                                                                           <  5:192899/74‑1 (MQ=255)
                                                                                                gCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGTGTGGGCCTGGAGCCCGCACCCGCCAACCTGCGCGAt                                                                                                                           <  5:197933/74‑1 (MQ=255)
                                                                                                gCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGTGTGGGCCTGGAGCCCGCACCCGCCAACCTGCGCGAt                                                                                                                           >  6:197933/1‑74 (MQ=255)
                                                                                                gCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGTGTGGGCCTGGAGCCCGCACCCGCCAACCTGCGCGAt                                                                                                                           >  6:192899/1‑74 (MQ=255)
                                                                                                  cACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGcc                                              <  1:137538/149‑1 (MQ=255)
                                                                                                  cACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGcc                                              >  6:336244/1‑149 (MQ=255)
                                                                                                  cACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCCGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGcc                                              <  4:311947/149‑1 (MQ=255)
                                                                                                   aCGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCcg                                             <  4:272453/149‑1 (MQ=255)
                                                                                                    cGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTc                                                                                             <  3:268807/100‑1 (MQ=255)
                                                                                                    cGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTc                                                                                             >  4:268807/1‑100 (MQ=255)
                                                                                                    cGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCcgt                                            <  1:177131/149‑1 (MQ=255)
                                                                                                     ggCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCcgtc                                           <  6:108232/149‑1 (MQ=255)
                                                                                                           acCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGc                                     <  6:212456/149‑1 (MQ=255)
                                                                                                             cggcgccggcgATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCg                                   <  5:165200/149‑1 (MQ=255)
                                                                                                                   cggcgATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGc                             >  4:185636/1‑149 (MQ=255)
                                                                                                                     gcgATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACAccctggccctgg                                                                    <  4:308696/108‑1 (MQ=255)
                                                                                                                     gcgATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGCTg                           <  2:203860/149‑1 (MQ=255)
                                                                                                                       gATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGCTGga                         <  3:332481/149‑1 (MQ=255)
                                                                                                                          ggCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGCTGgacga                      >  4:278849/1‑149 (MQ=255)
                                                                                                                                            cTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGCTGGACGACCGCCAACGCTGGGAtgt    >  2:61100/1‑149 (MQ=255)
                                                                                                                                              ggAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGCTGGACGACCGCCAACGCTGGGATGTgg  >  3:118033/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CGCGACTGGCGGTGGCCTACGAGATCAGCCAGGCCCCGGTGATTACCCCGGACCCCGCCCGCGGCGCTGCGCTGTACGCACAGAACTGCTCGATCTGCCACGGCGACACCGGCGCCGGCGATGGCCCGGCTGGCGTGGGCCTGGAGCCTGCACCCGCCAACCTGCGCGATGTCTCGCGCCTGGACCAGTTAAGCCTGTTCGACCTCTATAACACCCTGGCCCTGGGTATAGACGGCACCGAAATGCCGTCGTTTGCCGACCAGCTGGACGACCGCCAACGCTGGGATGTGG  >  NC_002947/232257‑232547

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: ‑