Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,270,150 C→T 45.2% A318V (GCA→GTA)  metZ → O‑succinylhomoserine sulfhydrylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,270,1500CT45.2% 27.2 / 112.5 93A318V (GCA→GTA) metZO‑succinylhomoserine sulfhydrylase
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base T (20/22);  total (41/52)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-02

CTGCGTATGCGTGCGCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGGTCGTCTGTCGCC  >  NC_002947/2270005‑2270297
                                                                                                                                                 |                                                                                                                                                    
cttcGTATGCGTGCGGACTGCGAAAGCGCCCAGGCGCTGGCCGTATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCAGGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGt                                                                                                                                                   <  5:210423/146‑1 (MQ=255)
   cGTATGCGTGCGCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAACGCCAGATGAGCGGTTTTGGTGCAGTg                                                                                                                                                  <  5:329649/147‑1 (MQ=255)
         cgtgcgCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCt                                                                                                                                           <  2:60004/148‑1 (MQ=255)
         cgtgcgCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCt                                                                                                                                           <  5:292426/148‑1 (MQ=255)
          gtgcgCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCa                                                                                                                                              >  4:34410/1‑144 (MQ=255)
          gtgcgCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCttt                                                                                                                                         >  6:305450/1‑149 (MQ=255)
          gtgcgCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCttt                                                                                                                                         <  5:305450/149‑1 (MQ=255)
          gtgcgCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGCGCAGTGGTCa                                                                                                                                              <  3:34410/144‑1 (MQ=255)
             cgcACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGaa                                                                                                                                      <  8:28159/149‑1 (MQ=255)
             cgcACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGaa                                                                                                                                      <  5:173409/149‑1 (MQ=255)
             cgcACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGaa                                                                                                                                      <  5:23612/149‑1 (MQ=255)
                 cTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCa                                                                                                                                  <  2:303851/149‑1 (MQ=255)
                           gcTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCa                                                                                                                        >  7:280448/1‑148 (MQ=255)
                                ggCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGGGGAGAAGGTGCATTACGCCGGCCGGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAagatg                                                                                                                  <  6:290827/149‑3 (MQ=255)
                                    cTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGgcgc                                                                                                              <  8:108159/149‑1 (MQ=255)
                                      ggCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCt                                                                                                            >  5:67199/1‑149 (MQ=255)
                                       gCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTg                                                                                                           >  3:280167/1‑149 (MQ=255)
                                           aaTGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGt                                                                                                       >  4:11573/1‑149 (MQ=255)
                                            aTGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGtt                                                                                                      <  2:52504/149‑1 (MQ=255)
                                               gCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTc                                                                                                    <  3:145855/148‑1 (MQ=255)
                                                cTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGGCAAAGAGGGCGCCTGGCGTTTCa                                                                                                   >  2:120219/1‑149 (MQ=255)
                                                cTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTcc                                                                                                   >  8:118029/1‑147 (MQ=255)
                                                   gAGCAGCTTGACGGGGTGGAGAAGGTTCATTACGCCGGCCTGCCCAGCCACCCTCAGCACGAACTGGGCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCgt                                                                                               <  2:14359/149‑2 (MQ=255)
                                                   gAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATc                                                                                                 <  5:163707/147‑1 (MQ=255)
                                                    agcagcAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGAc                                                                                              <  7:98545/149‑1 (MQ=255)
                                                        gcagGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTa                                                                                          >  1:223312/1‑149 (MQ=255)
                                                          agGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTAcc                                                                                        >  1:157978/1‑149 (MQ=255)
                                                           ggACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCaaa                                                                                                                      <  6:161969/118‑1 (MQ=255)
                                                           ggACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCaaa                                                                                                                      >  5:161969/1‑118 (MQ=255)
                                                           ggACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTAccc                                                                                       <  6:67199/149‑1 (MQ=255)
                                                            gACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACg                                                                                             <  3:303540/142‑1 (MQ=255)
                                                            gACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACg                                                                                             >  4:303540/1‑142 (MQ=255)
                                                             aCGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGa                                                                                     <  7:86761/149‑1 (MQ=255)
                                                             aCGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGa                                                                                     <  1:67527/149‑1 (MQ=255)
                                                             aCGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGa                                                                                     <  3:24969/149‑1 (MQ=255)
                                                                  gTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAagag                                                                                                                   >  4:5850/1‑114 (MQ=255)
                                                                  gTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAagag                                                                                                                   <  3:5850/114‑1 (MQ=255)
                                                                     gagaAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGAttt                                                                              <  3:11573/148‑1 (MQ=255)
                                                                     gagaAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGAttt                                                                              <  5:161256/148‑1 (MQ=255)
                                                                         aGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCAt                                                                          <  1:120219/149‑1 (MQ=255)
                                                                         aGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGCGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATc                                                                                                 <  7:213094/125‑1 (MQ=255)
                                                                         aGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGAGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATc                                                                                                 >  8:213094/1‑125 (MQ=255)
                                                                          ggTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCa                                                                        <  7:31896/149‑1 (MQ=255)
                                                                           gTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGt                                                                                                                                                >  8:256576/1‑77 (MQ=255)
                                                                           gTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTGTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCAc                                                                       >  5:58289/1‑149 (MQ=255)
                                                                            tGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCAc                                                                       >  1:328468/1‑148 (MQ=255)
                                                                                ttACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGt                                                                                   >  4:82315/1‑132 (MQ=255)
                                                                                ttACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGt                                                                                   <  3:82315/132‑1 (MQ=255)
                                                                                    gCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCt                                                              >  7:4562/1‑149 (MQ=255)
                                                                                      cGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTgg                                                            <  8:96876/149‑1 (MQ=255)
                                                                                      cGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGGGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTgg                                                            <  8:280448/149‑1 (MQ=255)
                                                                                            gCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGtt                                                                                                      >  6:46964/1‑101 (MQ=255)
                                                                                            gCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGtt                                                                                                      <  5:46964/101‑1 (MQ=255)
                                                                                             gccAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCCTCGCCGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAg                                                     >  7:247633/2‑149 (MQ=255)
                                                                                             cccAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGAc                                                       <  6:56381/147‑1 (MQ=255)
                                                                                             cccAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGAc                                                       <  2:157978/147‑1 (MQ=255)
                                                                                             cccAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGAc                                                       >  5:56381/1‑147 (MQ=255)
                                                                                             cccAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAg                                                     <  3:89528/149‑1 (MQ=255)
                                                                                             cccAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGTGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACg                                                                      <  5:80443/132‑1 (MQ=255)
                                                                                             cccAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGCTATTAGTACACTGGTCATCTTTTAAATAAAGTTGGGCAAAGAGAGAGCGTGGCGTTACAACGACGCTACACAACTGATTTCCATAACg                                                                      >  6:80443/1‑132 (MQ=255)
                                                                                               cAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGc                                                    >  3:86697/1‑148 (MQ=255)
                                                                                                      ccGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTg                                                                                                           >  2:307824/1‑86 (MQ=255)
                                                                                                      ccGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTg                                                                                                           <  1:307824/86‑1 (MQ=255)
                                                                                                       cGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAg                                                     >  5:69289/1‑139 (MQ=255)
                                                                                                       cGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAg                                                     <  6:69289/139‑1 (MQ=255)
                                                                                                            caCGAACTGGCCAAGCGCCAGATGAGCGGTGTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGc                                      <  6:58289/149‑1 (MQ=255)
                                                                                                                aaCTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCa                                                                           <  3:110741/108‑1 (MQ=255)
                                                                                                                aaCTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCa                                                                           >  4:110741/1‑108 (MQ=255)
                                                                                                                aaCTGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCAt                                  >  1:314833/1‑149 (MQ=255)
                                                                                                                aaCTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAac                                               >  5:239569/1‑136 (MQ=255)
                                                                                                                aaCTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAac                                               <  6:239569/136‑1 (MQ=255)
                                                                                                                   tGGCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCg                               >  3:82374/1‑149 (MQ=255)
                                                                                                                   tGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGt                                                          >  1:290092/1‑122 (MQ=255)
                                                                                                                   tGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGt                                                          <  2:290092/122‑1 (MQ=255)
                                                                                                                   tGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCg                               >  8:195304/1‑149 (MQ=255)
                                                                                                                     gCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGc                             >  5:119960/1‑149 (MQ=255)
                                                                                                                     gCCAAGCGCCAGATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGc                             <  4:86697/149‑1 (MQ=255)
                                                                                                                        aaGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGt                                                                                   >  8:171932/1‑92 (MQ=255)
                                                                                                                        aaGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGt                                                                                   <  7:171932/92‑1 (MQ=255)
                                                                                                                        aaGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGac                          >  1:236641/1‑149 (MQ=255)
                                                                                                                            gcCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACAAACCTGGGTGACAGCAAAACCACCATCGc                                      >  5:307211/1‑133 (MQ=255)
                                                                                                                            gcCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACAAACCTGGGTGACAGCAAAACCACCATCGc                                      <  6:307211/133‑1 (MQ=255)
                                                                                                                                gATGAGCGGTTTTGGTGTAGTGGTCAGCTTTGGAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTca                   >  6:216239/1‑148 (MQ=255)
                                                                                                                                  tGAGCGGTTTTGGTGTAGTGGGCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCacac                <  2:73305/149‑1 (MQ=255)
                                                                                                                                  tGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCacac                >  7:214843/1‑149 (MQ=255)
                                                                                                                                    aGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTCTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACgg              <  4:166144/149‑1 (MQ=255)
                                                                                                                                        gTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGgtcgt          >  3:248798/1‑149 (MQ=255)
                                                                                                                                         ttttGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGgtcgtc         <  2:236641/149‑1 (MQ=255)
                                                                                                                                           ttGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGGTCGTCt        <  7:174753/148‑1 (MQ=255)
                                                                                                                                              gtgtAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGGTCGTCTGTc     <  2:85988/148‑1 (MQ=255)
                                                                                                                                              gTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGGTCGTCTGTc     <  7:195304/148‑1 (MQ=255)
                                                                                                                                                gCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTg                                                             <  1:278587/90‑1 (MQ=255)
                                                                                                                                                gCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTg                                                             >  2:278587/1‑90 (MQ=255)
                                                                                                                                                gCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGGTCGTCTGTcgcc  >  6:244968/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGCGTATGCGTGCGCACTGCGAAAGCGCTCAGGCCCTGGCCGAATGGCTGGAGCAGCAGGACGGCGTGGAGAAGGTGCATTACGCCGGCCTGCCCAGCCACCCGCAGCACGAACTGGCCAAGCGCCAGATGAGCGGTTTTGGTGCAGTGGTCAGCTTTGAAGTCAAGGGGGGCAAAGAGGGCGCCTGGCGTTTCATCGACGCTACCCGAGTGATTTCCATCACGACCAACCTGGGTGACAGCAAAACCACCATCGCTCATCCGGCGACCACCTCACACGGTCGTCTGTCGCC  >  NC_002947/2270005‑2270297

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: