Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,252,650 A→C 7.7% D184A (GAC→GCC)  leuD → 3‑isopropylmalate dehydratase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,252,6500AC 7.7% 280.0 / 11.6 91D184A (GAC→GCC) leuD3‑isopropylmalate dehydratase small subunit
Reads supporting (aligned to +/- strand):  ref base A (44/40);  new base C (3/4);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GTCGACGAGCTGTTCAAGCAGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGT  >  NC_002947/2252502‑2252796
                                                                                                                                                     |                                                                                                                                                   
gTCGACGAGCTGTTCAAGCAGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGa                                                                                                                                                     <  1:40469/149‑1 (MQ=255)
   gacgaGCTGTTCAAGCAGGTCGAAGCCAACCCGGGCTACCAGCTGACAATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGa                                                                                                                                                  >  5:191292/1‑149 (MQ=255)
   gacgaGCTGTTCAAGCAGGTCGAAGCCAACACGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGa                                                                                                                                                  >  1:243665/1‑149 (MQ=255)
           gTTCAAGCAGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCAGGCCGATATCgtc                                                                                                                                           >  3:302304/1‑146 (MQ=255)
               aaGCAGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCAGGCCGATATCGTCCAGAc                                                                                                                                     <  4:302304/149‑1 (MQ=255)
                aGCAGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGACATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCAGGCCGATATCGTCCAGAcc                                                                                                                                    <  5:219788/149‑1 (MQ=255)
                   aGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCAGGCCGATATCGTCCAGAccgtt                                                                                                                                 <  2:380010/149‑4 (MQ=255)
                   aGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCAGGCCGATATCGTCCAGAccgtt                                                                                                                                 >  1:380010/1‑146 (MQ=255)
                    ggTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGc                                                                                                                                <  4:215328/149‑1 (MQ=255)
                     gTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCa                                                                                                                               <  5:48074/149‑1 (MQ=255)
                     ctCGAAGCCAAACCGGGGTACCACCTGACCATAGACCTGCAGGCGCAGCCGCTGACCCCGCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCCCAAGCACTGCCAGCTCACCGGTCTGGACGATATCGGCCGGACGTTg                                                                                                                                 >  4:201386/2‑147 (MQ=255)
                      tggAATCCAACCCGGGGTACCAGCTGTACATCGCACTTCAGGCGGCCGGGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAg                                                                                                                              <  6:296985/147‑1 (MQ=255)
                      tCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAg                                                                                                                              >  4:150263/1‑149 (MQ=255)
                        gAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGa                                                                                                                            >  4:186639/1‑149 (MQ=255)
                             cAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGc                                                                                                                        >  2:15193/1‑148 (MQ=255)
                              aaCCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGa                                                                                                                      >  4:243088/1‑149 (MQ=255)
                              aaCCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGa                                                                                                                      >  3:388382/1‑149 (MQ=255)
                                 ccGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGAcgc                                                                                                                   <  2:243665/149‑1 (MQ=255)
                                   gggCTACCAGCTGACCATCGAGCTGCTGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGACATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGGCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACTCGa                                                                                                                 <  3:201386/149‑1 (MQ=255)
                                   gggCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGa                                                                                                                 <  2:57317/149‑1 (MQ=255)
                                           aGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGcc                                                                                                         >  5:23055/1‑149 (MQ=255)
                                             cTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCtt                                                                                                       <  3:186639/149‑1 (MQ=255)
                                               gACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCa                                                                                                                               >  4:234800/1‑123 (MQ=255)
                                               gACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCa                                                                                                                               <  3:234800/123‑1 (MQ=255)
                                                 ccATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGaa                                                                                                   <  8:138441/149‑1 (MQ=255)
                                                   aTCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGAAGCGATCAAGGCCTTCGAAg                                                                                                  >  5:106210/1‑148 (MQ=255)
                                                      gACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCaa                                                                                              <  3:301034/149‑1 (MQ=255)
                                                        ccTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGc                                                                                            <  1:15193/149‑1 (MQ=255)
                                                         cTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGAcc                                                                                                                                    <  8:33128/108‑1 (MQ=255)
                                                         cTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGAcc                                                                                                                                    >  7:33128/1‑108 (MQ=255)
                                                         cTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGTAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCa                                                                                           >  6:397794/1‑149 (MQ=255)
                                                         cTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCAGGCCGATATCGTCCAGAcc                                                                                                                                    <  1:80393/108‑1 (MQ=255)
                                                         cTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCAGGCCGATATCGTCCAGAcc                                                                                                                                    >  2:80393/1‑108 (MQ=255)
                                                                cgcaggcgGTGACCCGCCCCGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCAGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCgcgc                                                                                     >  5:228528/1‑148 (MQ=255)
                                                                 gcaggcgGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATTCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCgcgc                                                                                     >  1:291780/1‑148 (MQ=255)
                                                                 gcaggcgGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCgcca                                                                                   <  4:388382/149‑1 (MQ=255)
                                                                     gcgGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCgccagc                                                                                 >  4:393613/1‑147 (MQ=255)
                                                                     gcgGTGACACGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCgccagc                                                                                 >  2:243362/1‑147 (MQ=255)
                                                                            cccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAgg                                                                                                           <  4:13406/114‑1 (MQ=255)
                                                                            cccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAgg                                                                                                           >  3:13406/1‑114 (MQ=255)
                                                                            cccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGcc                                                                                                >  5:85959/1‑125 (MQ=255)
                                                                            cccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGcc                                                                                                <  6:85959/125‑1 (MQ=255)
                                                                            cccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTgg                                                                        <  2:124920/149‑1 (MQ=255)
                                                                             ccgcccgGATGGCAATGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCTATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGc                                                                       <  6:228528/149‑1 (MQ=255)
                                                                             ccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCgccagcc                                                                                <  3:23313/140‑1 (MQ=255)
                                                                             ccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCgccagcc                                                                                >  4:23313/1‑140 (MQ=255)
                                                                             ccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGc                                                                       >  2:348179/1‑149 (MQ=255)
                                                                             ccgcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGc                                                                       >  6:19795/1‑149 (MQ=255)
                                                                               gcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTgg                                                                        <  4:77675/146‑1 (MQ=255)
                                                                               gcccgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTgg                                                                        >  3:77675/1‑146 (MQ=255)
                                                                                  cgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACcgc                                                                                       >  8:347801/1‑128 (MQ=255)
                                                                                  cgGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACcgc                                                                                       <  7:347801/128‑1 (MQ=255)
                                                                                     aTGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGt                                                               <  3:150263/149‑1 (MQ=255)
                                                                                       ggCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGa                                                             >  3:214647/1‑149 (MQ=255)
                                                                                         cAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCgcc                                                                                    <  2:20343/124‑1 (MQ=255)
                                                                                         cAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCgcc                                                                                    >  1:20343/1‑124 (MQ=255)
                                                                                         cAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCtt                                                                          <  8:240486/134‑1 (MQ=255)
                                                                                         cAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCtt                                                                          >  7:240486/1‑134 (MQ=255)
                                                                                                ctgcACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGaa                                                     <  1:348179/148‑1 (MQ=255)
                                                                                                ctgcACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGaa                                                     <  1:243362/148‑1 (MQ=255)
                                                                                                  gcACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAgg                                                                                                           <  6:245696/92‑1 (MQ=255)
                                                                                                  gcACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAgg                                                                                                           >  5:245696/1‑92 (MQ=255)
                                                                                                   cACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGTAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTc                                                 <  5:397794/149‑1 (MQ=255)
                                                                                                      ttCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCt                                                                           <  1:5949/120‑1 (MQ=255)
                                                                                                      ttCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCt                                                                           >  2:5949/1‑120 (MQ=255)
                                                                                                        cGAGTTCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTa                                              <  2:291780/147‑1 (MQ=255)
                                                                                                        cGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTa                                              <  6:191292/147‑1 (MQ=255)
                                                                                                           gaTCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             >  5:131821/1‑145 (MQ=255)
                                                                                                           gaTCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             <  6:131821/145‑1 (MQ=255)
                                                                                                               cGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAg                                      >  6:388329/1‑149 (MQ=255)
                                                                                                                     gTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGAACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCgg                                 >  1:90820/1‑149 (MQ=255)
                                                                                                                     gTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCggg                                >  7:388133/1‑149 (MQ=255)
                                                                                                                        ccGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGcct                             >  6:152448/1‑149 (MQ=255)
                                                                                                                                aCTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTc                     <  2:43061/149‑1 (MQ=255)
                                                                                                                                       gcTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAg               >  2:64156/1‑148 (MQ=255)
                                                                                                                                        cTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGa              >  5:7253/1‑148 (MQ=255)
                                                                                                                                         tCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGa              >  1:160861/1‑147 (MQ=255)
                                                                                                                                          cAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAg           >  5:376736/1‑149 (MQ=255)
                                                                                                                                           aaCGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGa                                         >  6:257177/1‑118 (MQ=255)
                                                                                                                                           aaCGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGa                                         <  5:257177/118‑1 (MQ=255)
                                                                                                                                           aaCGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGaaa            <  4:214647/147‑1 (MQ=255)
                                                                                                                                            aCGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGtt         <  6:106210/149‑1 (MQ=255)
                                                                                                                                               gTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTAcg      >  1:51989/1‑149 (MQ=255)
                                                                                                                                                tCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTAcgc     <  2:90820/149‑1 (MQ=255)
                                                                                                                                                   ggACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGt  >  6:332230/1‑149 (MQ=255)
                                                                                                                                                    gACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             <  7:359215/105‑1 (MQ=255)
                                                                                                                                                    gACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             >  8:359215/1‑105 (MQ=255)
                                                                                                                                                    gACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             <  8:343144/105‑1 (MQ=255)
                                                                                                                                                    gACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             <  6:357205/105‑1 (MQ=255)
                                                                                                                                                    gACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             >  7:343144/1‑105 (MQ=255)
                                                                                                                                                    gACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAg                                             >  5:357205/1‑105 (MQ=255)
                                                                                                                                                    gACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCgg   <  8:388133/147‑1 (MQ=255)
                                                                                                                                                     aCGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTc                                                 >  2:385537/1‑100 (MQ=255)
                                                                                                                                                     aCGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTc                                                 <  1:385537/100‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTCGACGAGCTGTTCAAGCAGGTCGAAGCCAACCCGGGCTACCAGCTGACCATCGACCTGCAGGCGCAGGCGGTGACCCGCCCGGATGGCAAGGTGCTGCACTTCGAGATCGATGCGTTCCGCAAGCACTGCCTGCTCAACGGTCTGGACGATATCGGCCTGACCTTGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGT  >  NC_002947/2252502‑2252796

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: