Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,949,791 T→C 8.2% N416D (AAT→GAT)  PP_3486 ← cytochrome c

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,949,7910TC 8.2% 210.2 / 5.5 73N416D (AAT→GAT) PP_3486cytochrome c
Reads supporting (aligned to +/- strand):  ref base T (35/32);  new base C (3/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGACCACTGCGCCGGCAAGCGCGCCA  >  NC_002947/3949643‑3949937
                                                                                                                                                    |                                                                                                                                                  
cggcgACCAGGGCCAAAGCCCAGTGCCGGTGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTAtt                                                                                                                                                    <  6:280859/149‑1 (MQ=255)
cggcgACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTAtt                                                                                                                                                    <  7:294012/149‑1 (MQ=255)
    gACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGAcc                                                                                                                                                >  4:103772/1‑149 (MQ=255)
       cAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTg                                                                                                                                              >  7:284053/1‑148 (MQ=255)
           gCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAAt                                                                                                                                          >  5:198971/1‑148 (MQ=255)
            ccAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGa                                                                                                                                        >  8:15262/1‑149 (MQ=255)
            ccAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGa                                                                                                                                        <  5:264754/149‑1 (MQ=255)
             cAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGAc                                                                                                                                       >  4:139774/1‑149 (MQ=255)
             cAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGAc                                                                                                                                       <  4:297489/149‑1 (MQ=255)
             cAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGAc                                                                                                                                       >  2:24930/1‑149 (MQ=255)
                aGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGcc                                                                                                                                    >  2:258160/1‑149 (MQ=255)
                    cAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCt                                                                                                                                               >  1:345387/1‑134 (MQ=255)
                    cAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCt                                                                                                                                               <  2:345387/134‑1 (MQ=255)
                    cAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGgc                                                                                                                                >  2:341514/1‑149 (MQ=255)
                      gTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGgcgc                                                                                                                              <  3:103772/149‑1 (MQ=255)
                          cGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGAt                                                                                                                          >  1:59258/1‑149 (MQ=255)
                           ggggTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATc                                                                                                                         <  8:268159/149‑1 (MQ=255)
                            gggTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATcc                                                                                                                        <  4:117686/149‑1 (MQ=255)
                              gTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAg                                                                                                                      >  7:37887/1‑149 (MQ=255)
                               tCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAgg                                                                                                                     <  2:69953/149‑1 (MQ=255)
                               tCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAgg                                                                                                                     <  8:37887/149‑1 (MQ=255)
                                   gCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATg                                                                                                                                         <  5:123639/125‑1 (MQ=255)
                                   gCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATg                                                                                                                                         >  6:123639/1‑125 (MQ=255)
                                     cGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTgg                                                                                                               <  1:24930/149‑1 (MQ=255)
                                            cacaTCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGCCGCCTGGCGCGGATCCAGGCCCTGGCgctgct                                                                                                        >  6:351304/1‑149 (MQ=255)
                                             acaTCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCgctgct                                                                                                        <  6:198971/148‑1 (MQ=255)
                                             acaTCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGctgctg                                                                                                       >  7:48038/1‑149 (MQ=255)
                                             acaTCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGctgctg                                                                                                       <  5:169419/149‑1 (MQ=255)
                                                   tGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGccc                                                                                                 >  8:322720/1‑149 (MQ=255)
                                                       ggCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGc                                                                                              <  5:380098/148‑1 (MQ=255)
                                                       ggCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGc                                                                                              >  6:380098/1‑148 (MQ=255)
                                                        gCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGa                                                                                            <  3:215451/149‑1 (MQ=255)
                                                        gCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGa                                                                                            <  2:59258/149‑1 (MQ=255)
                                                          cTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGc                                                                                              <  6:174893/145‑1 (MQ=255)
                                                          cTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGc                                                                                              >  5:174893/1‑145 (MQ=255)
                                                             cGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAgcgc                                                                                        <  1:258160/148‑1 (MQ=255)
                                                                  gcagcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCa                                                                                       >  4:56573/1‑144 (MQ=255)
                                                                  gcagcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAg                                                                                  <  8:284053/149‑1 (MQ=255)
                                                                  gcagcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACCCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCa                                                                                       <  3:56573/144‑1 (MQ=255)
                                                                    agcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGc                                                                                     >  7:81062/1‑144 (MQ=255)
                                                                    agcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGc                                                                                     <  8:81062/144‑1 (MQ=255)
                                                                    agcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGc                                                                                <  5:234708/149‑1 (MQ=255)
                                                                    agcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGc                                                                                >  6:234708/1‑149 (MQ=255)
                                                                     gcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGcc                                                                               >  4:151672/1‑149 (MQ=255)
                                                                     gcaATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGcc                                                                               >  6:156170/1‑149 (MQ=255)
                                                                       aaTCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCtt                                                                             >  1:15749/1‑149 (MQ=255)
                                                                            gTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCg                                                                        >  7:273208/1‑149 (MQ=255)
                                                                            gTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCg                                                                        >  7:332810/1‑149 (MQ=255)
                                                                                   gcgGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCTGATAGCGTATCGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGACAGGCCTTGCTCGGTCCAGt                                                                 >  5:259186/1‑149 (MQ=255)
                                                                                   gcgGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATCGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGt                                                                 >  5:259182/1‑149 (MQ=255)
                                                                                   gcgGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATCGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTAGGTCCAGt                                                                 >  7:308191/1‑149 (MQ=255)
                                                                                        tGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCt                                                              >  2:36961/1‑147 (MQ=255)
                                                                                               ccGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACg                                                      >  5:176498/1‑148 (MQ=255)
                                                                                                 gTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAg                                                   >  7:196694/1‑149 (MQ=255)
                                                                                                   ccAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCg                                                                                                >  2:243506/1‑102 (MQ=255)
                                                                                                   ccAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCg                                                                                                <  1:243506/102‑1 (MQ=255)
                                                                                                    cAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAggg                                                 <  1:341514/148‑1 (MQ=255)
                                                                                                      aTCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATCGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTAGGTCCAGTCCTCAGGCACGTAGGGGTa                                              <  8:308191/149‑1 (MQ=255)
                                                                                                      aTCTGCTTCAGCGCGTCGCTCAGCAGGCGCAACGCGCATAGCGTATCGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTa                                              <  6:259182/149‑1 (MQ=255)
                                                                                                      aTCTGCTTCAGCGCATCGCTAAGCAGGCGCAACGCGCATAGCGTATCGACATGCAATGACGCCTTGCGCGGATCCAGTCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTa                                              <  6:259186/149‑1 (MQ=255)
                                                                                                       tCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAg                                                           <  4:239611/135‑1 (MQ=255)
                                                                                                       tCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAg                                                           >  3:239611/1‑135 (MQ=255)
                                                                                                           cTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAg                                             <  7:311119/145‑1 (MQ=255)
                                                                                                              cAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCCGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGa                                       >  6:79901/1‑148 (MQ=255)
                                                                                                                gcgcTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAAc                                    <  3:151672/149‑1 (MQ=255)
                                                                                                                gcgcTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAAc                                    <  5:156170/149‑1 (MQ=255)
                                                                                                                   cTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGc                                 >  1:282354/1‑149 (MQ=255)
                                                                                                                    ttCGCTCAGCAGGCGCAACGCGGATAGAGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCATCAGGCACGTAGGGGTAGGCCAGAAACAGc                                 >  7:234132/1‑148 (MQ=255)
                                                                                                                           agcagGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAATCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGc                           <  7:322720/147‑1 (MQ=255)
                                                                                                                            gcagGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGAc                        >  2:366813/1‑149 (MQ=255)
                                                                                                                             cagGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGAcc                       <  2:15749/149‑1 (MQ=255)
                                                                                                                             cagGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGAc                        >  2:169914/1‑148 (MQ=255)
                                                                                                                                      cgcgCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGACCACTGCGCCg              <  8:332810/149‑1 (MQ=255)
                                                                                                                                            tAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGACCACTGCGCCGGCAAgc        <  1:366813/149‑1 (MQ=255)
                                                                                                                                              gCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGACACTGCGCCGGCAAgcgc      <  2:317164/148‑1 (MQ=255)
                                                                                                                                                   ttGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGACCACTGCGCCGGCAAGCGCGCca  <  1:169914/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTTCAGCGCTTCGCTCAGCAGGCGCAACGCGCATAGCGTATTGACCTGCAATGACGCCTGGCGCGGATCCAGGCCCTGGCGCTGCTGCAGCCCGGCGAGCGCAGCCAGGCCTTGCTCGGTCCAGTCCTCAGGCACGTAGGGGTAGGCCAGAAACAGCCGCTGCGACCACTGCGCCGGCAAGCGCGCCA  >  NC_002947/3949643‑3949937

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: