Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,161,206 A→T 9.4% Q26L (CAG→CTG)  PP_2776 → homocysteine S‑methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,161,2060AT 9.4% 163.9 / 5.3 61Q26L (CAG→CTG) PP_2776homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base T (3/3);  total (30/31)
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

GGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCTTTGCGCGC  >  NC_002947/3161064‑3161354
                                                                                                                                                |                                                                                                                                                    
gggcgccgcCAATCAGGCGACCGTCGATCGTTGGGTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCGCGATGGTGGGATGGGACGAGAACTCACGCGGATCGGCGCACCCTTCCGGCAGCCg                                                                                                                                                  <  4:160240/144‑1 (MQ=255)
 gTCGCCGCCAATCAGGCGACCGTCGCTCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGaa                                                                                                                                               <  7:151272/149‑1 (MQ=255)
    gccgccAATCAGGCGTCCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCTGCCGGaagg                                                                                                                                             <  3:4302/148‑3 (MQ=255)
     ccgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGt                                                                                                                                           <  5:227509/149‑1 (MQ=255)
      cgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTc                                                                                                                                          <  8:297281/149‑1 (MQ=255)
       gccAATCAGGCGACCGGCGATCGTTGGTTTTTTTCCACTGGCACCCATGTGGCGTCCCCCAACCCCCCCCCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGt                                                                                                                                           <  1:308234/149‑1 (MQ=255)
          aaTCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCa                                                                                                                                      >  4:355268/1‑149 (MQ=255)
                gCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGgcg                                                                                                                                >  2:99910/1‑149 (MQ=255)
                         gATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGaa                                                                                                                       <  6:334802/149‑1 (MQ=255)
                             gttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGATTGGTCGGCACTGGCGCTGTTCGAAgcgc                                                                                                                   <  1:99910/149‑1 (MQ=255)
                             gttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTg                                                                                                                             <  2:334917/139‑1 (MQ=255)
                             gttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTg                                                                                                                             >  1:334917/1‑139 (MQ=255)
                              ttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGcc                                                                                                                  >  5:357227/1‑149 (MQ=255)
                                   tgggTGCCAGTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGaa                                                                                                              >  8:176546/4‑148 (MQ=255)
                                         ccACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGa                                                                                                       <  7:87414/149‑1 (MQ=255)
                                              ggCAGCCACCTGGG‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTGGATGGGGGGATGGGGCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGGCCGAATACGTGACGCAg                                                                                                  >  1:269453/1‑149 (MQ=255)
                                               gCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAgg                                                                                                 >  8:226835/1‑149 (MQ=255)
                                                   ccAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCa                                                                                             >  4:340358/1‑149 (MQ=255)
                                                    cAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCAc                                                                                            <  3:340358/149‑1 (MQ=255)
                                                             t‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGt                                                                                     >  4:184379/1‑148 (MQ=255)
                                                             t‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGt                                                                                     >  2:182547/1‑148 (MQ=255)
                                                                        cTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTgccgcc                                                                          <  3:355268/149‑1 (MQ=255)
                                                                         tCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGccgccg                                                                         >  7:203847/1‑149 (MQ=255)
                                                                         tCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGccgccg                                                                         >  7:181834/1‑149 (MQ=255)
                                                                           ggCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgg                                                                        <  3:266092/148‑1 (MQ=255)
                                                                           ggCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgg                                                                        <  7:176546/148‑1 (MQ=255)
                                                                            gCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgggg                                                                      >  8:221425/1‑149 (MQ=255)
                                                                                 ccTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGTGCCAgg                                                                 <  3:184379/149‑1 (MQ=255)
                                                                                 ccTGGACTTTGCTCGATGGTGGGAGGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAgg                                                                 <  2:242178/149‑1 (MQ=255)
                                                                                 ccTGGACGTTGCTCGATGGTGGTATGGGACGAGAACTCAAGCGGAGCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGACGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAgg                                                                 <  4:33671/149‑1 (MQ=255)
                                                                                  cTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAggg                                                                >  5:10316/1‑149 (MQ=255)
                                                                                       cTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATc                                                           <  1:182547/149‑1 (MQ=255)
                                                                                         acgCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcac                                                         <  7:64999/147‑1 (MQ=255)
                                                                                           gCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcacc                                                        <  7:221425/148‑1 (MQ=255)
                                                                                            cTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcacca                                                       <  4:351022/148‑1 (MQ=255)
                                                                                            cTCGATGGTGGGATGGGACAAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATc                                                           >  5:257693/1‑144 (MQ=255)
                                                                                            cTCGATGGTGGGATGGGACAAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATc                                                           <  6:257693/144‑1 (MQ=255)
                                                                                             tCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcaccac                                                      >  4:339069/1‑148 (MQ=255)
                                                                                                    tgggatgggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGc                                               >  1:99608/1‑148 (MQ=255)
                                                                                                      ggagcggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCCGCCCGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACTTTACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGTTACTCACCACAAACAGCTAc                                            >  8:174552/6‑149 (MQ=255)
                                                                                                       gatgggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTa                                                              >  8:263911/1‑130 (MQ=255)
                                                                                                       gatgggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTa                                                              <  7:263911/130‑1 (MQ=255)
                                                                                                               gagaACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGgtggtg                                   <  6:357227/149‑1 (MQ=255)
                                                                                                                 gaACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGcc                                 >  8:10411/1‑149 (MQ=255)
                                                                                                                   aCTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGt                               <  7:10411/149‑1 (MQ=255)
                                                                                                                     tCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcaccac                                                      >  7:214459/1‑124 (MQ=255)
                                                                                                                     tCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcaccac                                                      <  8:214459/124‑1 (MQ=255)
                                                                                                                         gCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTcaca                         >  2:328963/1‑149 (MQ=255)
                                                                                                                            gATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTg                      <  7:226835/149‑1 (MQ=255)
                                                                                                                                ggagcACCCTTCCGGCATCCGGAATGGTGGGGACTGGCGCTGTTGGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACGATGATCTACGCGGTGGTGCCGTTTCACATTGGCGa                  <  7:174552/146‑1 (MQ=255)
                                                                                                                                   cgaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAgg                >  7:18688/3‑148 (MQ=255)
                                                                                                                                    gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                  >  1:158149/2‑145 (MQ=255)
                                                                                                                                    gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                  <  2:163779/144‑1 (MQ=255)
                                                                                                                                    gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                  <  2:158149/144‑1 (MQ=255)
                                                                                                                                    gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                  >  1:163779/2‑145 (MQ=255)
                                                                                                                                    cACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGaa              >  6:348644/1‑149 (MQ=255)
                                                                                                                                         ttCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGttt                             <  5:35729/129‑1 (MQ=255)
                                                                                                                                         ttCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGttt                             >  6:35729/1‑129 (MQ=255)
                                                                                                                                         ttCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGACCGCtt         >  6:215431/1‑149 (MQ=255)
                                                                                                                                         ttCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCtt         >  1:237791/1‑149 (MQ=255)
                                                                                                                                               cAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCTTTgcgcg   >  4:190438/1‑149 (MQ=255)
                                                                                                                                                aGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCTTTgcgcgc  >  6:218164/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCTTTGCGCGC  >  NC_002947/3161064‑3161354

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑