Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,665,001 C→G 5.2% R45R (CGC→CGG nuoI → NADH‑quinone oxidoreductase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,665,0010CG 5.2% 251.7 / 3.7 77R45R (CGC→CGGnuoINADH‑quinone oxidoreductase subunit I
Reads supporting (aligned to +/- strand):  ref base C (39/34);  new base G (2/2);  total (41/36)
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.37e-01

AGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGTTCTTCC  >  NC_002947/4664855‑4665149
                                                                                                                                                  |                                                                                                                                                    
aGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGAACCGGCAACCAGCTGCGCAGCCTGGCAATGGTGTTCTACCACGGGTTCCGCAAGCGCGACGCCCTGCAATACCCCGAAGAACCAGTGTACCTGCCGCCtcgcaa                                                                                                                                                    >  7:132478/1‑147 (MQ=255)
       gACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGc                                                                                                                                             >  7:240469/1‑149 (MQ=255)
                  aaGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTccc                                                                                                                                  >  7:169839/1‑148 (MQ=255)
                  aaGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCt                                                                                                                                  >  3:113126/1‑149 (MQ=255)
                   aGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCACACGGGTTCCGCAAGGGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTc                                                                                                                                 <  8:240469/149‑1 (MQ=255)
                       tatCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCTATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTc                                                                                                                                    >  8:256299/1‑142 (MQ=255)
                       tatCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCTATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTc                                                                                                                                    <  7:256299/142‑1 (MQ=255)
                       tatCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCAccc                                                                                                                             <  4:127856/149‑1 (MQ=255)
                         tCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCGGCCGCGGTAGCGGGGCGGCag                                                                                                                                        <  8:42404/136‑2 (MQ=255)
                         tCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCGGCCGCGGTAGCGCGGCGGCag                                                                                                                                        >  7:42404/1‑135 (MQ=255)
                         tCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCCCACCcgc                                                                                                                           <  8:169839/149‑1 (MQ=255)
                              gACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGAccc                                                                                                                      <  3:320538/149‑1 (MQ=255)
                                        aGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCg                                                                                                                                      <  1:321457/123‑1 (MQ=255)
                                        aGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCg                                                                                                                                      >  2:321457/1‑123 (MQ=255)
                                          ggcacgggcagcCAGCTGCGCATCCTGGCAATGGTGTTGTACCACGGGTTCCGCAATCGCGACGCCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCgagga                                                                                                          <  8:132478/149‑1 (MQ=255)
                                          ggcaccggcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCgag                                                                                                            >  2:196609/1‑147 (MQ=255)
                                           gcaccggcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCgaggag                                                                                                         >  4:299849/1‑149 (MQ=255)
                                           gcaccggcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCgaggag                                                                                                         >  7:193119/1‑149 (MQ=255)
                                           gcaccggcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCgaggag                                                                                                         >  6:25337/1‑149 (MQ=255)
                                                 gcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGc                                                                                                   <  3:299849/149‑1 (MQ=255)
                                                   accCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGc                                                                                                                                             <  3:41217/105‑1 (MQ=255)
                                                   accCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGc                                                                                                                                             >  4:41217/1‑105 (MQ=255)
                                                     ccAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCg                                                                                                              <  2:173849/134‑1 (MQ=255)
                                                     ccAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCg                                                                                                              >  1:173849/1‑134 (MQ=255)
                                                     ccAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAg                                                                                               <  5:25337/149‑1 (MQ=255)
                                                     ccAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAg                                                                                               >  6:153288/1‑149 (MQ=255)
                                                     ccAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAg                                                                                               >  6:158461/1‑149 (MQ=255)
                                                      cAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGc                                                                                              <  4:113126/149‑1 (MQ=255)
                                                        gCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGt                                                                                            >  5:276348/1‑149 (MQ=255)
                                                           gcgcAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGc                                                                                          >  5:338710/1‑148 (MQ=255)
                                                                   tGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGAcc                                                                                                                       <  2:169969/111‑1 (MQ=255)
                                                                   tGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGAcc                                                                                                                       >  1:169969/1‑111 (MQ=255)
                                                                       aaTGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCgg                                                                             >  5:74623/1‑149 (MQ=255)
                                                                        aTGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGt                                                                            <  3:149018/149‑1 (MQ=255)
                                                                         tGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTg                                                                           >  7:168469/1‑149 (MQ=255)
                                                                           gtgtTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTgg                                                                          <  8:193119/148‑1 (MQ=255)
                                                                                ctcCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGt                                                                                            >  1:45239/1‑125 (MQ=255)
                                                                                ctcCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGt                                                                                            <  2:45239/125‑1 (MQ=255)
                                                                                ctcCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGcc                                                                    >  3:173938/1‑149 (MQ=255)
                                                                                 tcCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGccc                                                                   >  8:23991/1‑149 (MQ=255)
                                                                                   ccACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCgg                                                                 >  1:53493/1‑149 (MQ=255)
                                                                                   ccACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCgg                                                                 >  7:309585/1‑149 (MQ=255)
                                                                                   ccACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCgg                                                                 >  7:85351/1‑149 (MQ=255)
                                                                                     aCGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTg                                                               >  7:196791/1‑149 (MQ=255)
                                                                                            ccGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGc                                                         <  7:23991/148‑1 (MQ=255)
                                                                                                aaGCGCGACACCCTGCAATACCCCGGAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTc                                                    >  1:115266/1‑149 (MQ=255)
                                                                                                aaGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGAGGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTc                                                    >  8:86867/1‑149 (MQ=255)
                                                                                                 aGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCg                                                   <  4:173938/149‑1 (MQ=255)
                                                                                                    gcgACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTg                                                <  5:158461/149‑1 (MQ=255)
                                                                                                    gcgACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTg                                                <  5:153288/149‑1 (MQ=255)
                                                                                                    gcgACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTg                                                >  2:133827/1‑149 (MQ=255)
                                                                                                        cacCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGa                                            <  6:276348/149‑1 (MQ=255)
                                                                                                           ccTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGCTGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAgg                                         >  5:212930/1‑149 (MQ=255)
                                                                                                           ccTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGGGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAgg                                         >  3:261263/1‑149 (MQ=255)
                                                                                                               cAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGTTGCAGAAGGCCga                                     <  2:53493/149‑1 (MQ=255)
                                                                                                               cAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGcc                                       <  1:196609/147‑1 (MQ=255)
                                                                                                                 aTACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCgaga                                   <  6:338710/149‑1 (MQ=255)
                                                                                                                   aCCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCgaga                                   >  4:328112/1‑147 (MQ=255)
                                                                                                                        gaagaaCCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGa                            >  6:269526/1‑149 (MQ=255)
                                                                                                                         aagaaCCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGc                                                            <  8:198044/116‑1 (MQ=255)
                                                                                                                         aagaaCCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGc                                                            >  7:198044/1‑116 (MQ=255)
                                                                                                                          agaaCCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCTAGACAGAGTACg                          <  8:309585/149‑1 (MQ=255)
                                                                                                                              cccGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCgaga                                   <  6:254025/136‑1 (MQ=255)
                                                                                                                              cccGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCgaga                                   >  1:136198/1‑136 (MQ=255)
                                                                                                                              cccGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCgaga                                   >  5:254025/1‑136 (MQ=255)
                                                                                                                              cccGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCg                      >  2:233148/1‑149 (MQ=255)
                                                                                                                               ccGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGc                     <  8:85351/149‑1 (MQ=255)
                                                                                                                                 gtgtACCTGCCGCCGCGCTACCGCGGCCGCTTCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTg                   >  1:28269/1‑149 (MQ=255)
                                                                                                                                      gatGCGGCCGCGGTAGCGCGGCGGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCga                                     <  5:15976/124‑1 (MQ=255)
                                                                                                                                      gatGCGGCCGCGGTAGCGCGGCGGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCga                                     >  6:15976/3‑126 (MQ=255)
                                                                                                                                      ccTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGTTGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTAcc              <  6:212930/149‑1 (MQ=255)
                                                                                                                                      ccTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGc                     >  2:23620/1‑142 (MQ=255)
                                                                                                                                      ccTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGc                     <  1:23620/142‑1 (MQ=255)
                                                                                                                                          ccgccgCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCg            <  4:261263/147‑1 (MQ=255)
                                                                                                                                          ccgccgCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGa          <  5:269526/149‑1 (MQ=255)
                                                                                                                                           cgccgcGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAg         <  1:233148/149‑1 (MQ=255)
                                                                                                                                             ccgcGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACg                          >  2:231039/1‑130 (MQ=255)
                                                                                                                                             ccgcGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACg                          <  1:231039/130‑1 (MQ=255)
                                                                                                                                                 gcTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGttcttc   <  1:133827/149‑1 (MQ=255)
                                                                                                                                                  cTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGTTCTTcc  <  6:74623/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGTTCTTCC  >  NC_002947/4664855‑4665149

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: