Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,953,408 T→C 7.3% T45A (ACC→GCC)  mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,953,4080TC 7.3% 187.9 / 4.4 68T45A (ACC→GCC) mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base C (3/2);  total (34/34)
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

TGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGGTTG  >  NC_002947/1953261‑1953554
                                                                                                                                                    |                                                                                                                                                    
tGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTc                                                                                                                                                    >  2:61772/1‑149 (MQ=255)
tGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTc                                                                                                                                                    >  5:130031/1‑149 (MQ=255)
  aCGAAGTGCAGGCGCGAGTCGAGGTGTGAGTGCTGGTCGAACAGCTGCCCGGCGCAGAAAAAATGCATGCCGGTGCCAAAACCGGTTGCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTg                                                                                                                                                  <  6:309992/149‑1 (MQ=255)
    gAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGc                                                                                                                                                >  3:233569/1‑149 (MQ=255)
       gtacAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGa                                                                                                                                             <  7:21925/146‑1 (MQ=255)
       gTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGa                                                                                                                                             <  5:255686/149‑1 (MQ=255)
          cAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAACACGCTCGGCCAGACCCGCCTGGCCGAg                                                                                                                                            <  1:226045/148‑1 (MQ=255)
          cAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGAg                                                                                                                                            >  2:226045/1‑147 (MQ=255)
           aGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAgg                                                                                                                                           >  6:122923/1‑147 (MQ=255)
           aGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAgg                                                                                                                                           <  5:122923/147‑1 (MQ=255)
               gcgcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACg                                                                                                                                     <  4:233569/149‑1 (MQ=255)
                           gtgcgtgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTttct                                                                                                                         >  5:375263/1‑149 (MQ=255)
                              cgtgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGGTCGGCCAGGGGGGTCCTGGCCGAGGAACACGTGCTTGGTttctt                                                                                                                        <  1:239052/148‑1 (MQ=255)
                                     gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGa                                                                                                                                             <  1:173053/119‑1 (MQ=255)
                                     gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGa                                                                                                                                             >  2:173053/1‑119 (MQ=255)
                                     gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCtt                                                                                                              >  7:196132/1‑149 (MQ=255)
                                       cGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAg                                                                                                                                            <  6:333163/118‑1 (MQ=255)
                                       cGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAg                                                                                                                                            >  5:333163/1‑118 (MQ=255)
                                       cGAACAGCTGCCAGGCGCAGAAAAAATTCACGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCg                                                                                                            <  4:299559/149‑1 (MQ=255)
                                        gAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGt                                                                                                           >  6:152288/1‑149 (MQ=255)
                                         aaCAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                                          <  6:375263/149‑1 (MQ=255)
                                         aaCAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                                          >  3:195762/1‑149 (MQ=255)
                                            aGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGTCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGAc                                                                                                       >  8:77098/1‑149 (MQ=255)
                                            aGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGAc                                                                                                       >  8:267345/1‑149 (MQ=255)
                                              cTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                                          >  3:399358/1‑144 (MQ=255)
                                              cTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                                          <  4:399358/144‑1 (MQ=255)
                                                  cAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGTCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGaa                                                                                                 <  7:77098/149‑1 (MQ=255)
                                                   aGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGaaa                                                                                                >  1:374258/1‑149 (MQ=255)
                                                         aGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGAt                                                                                                                   >  6:380744/1‑124 (MQ=255)
                                                         aGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGAt                                                                                                                   <  5:380744/124‑1 (MQ=255)
                                                               aaTTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGt                                                                                                                             >  2:311619/1‑109 (MQ=255)
                                                               aaTTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGt                                                                                                                             <  1:311619/109‑1 (MQ=255)
                                                                   cATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACAAGTGCTTGGTTTCTTCGCATGCCTTCGTTGACTGCGAAATAGACGTCGTCATa                                                                                  <  1:61772/148‑1 (MQ=255)
                                                                       ccGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCcg                                                                            <  6:130031/149‑1 (MQ=255)
                                                                             ccAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATg                                                                      >  1:231243/1‑149 (MQ=255)
                                                                               aaaaCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggg                                                                    >  2:164795/1‑149 (MQ=255)
                                                                               aaaaCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggg                                                                    >  7:248429/1‑149 (MQ=255)
                                                                                  aCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACATCGTCATAGTGCCGCGAATgg                                                                     <  2:156706/145‑1 (MQ=255)
                                                                                  aCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACATCGTCATAGTGCCGCGAATgg                                                                     >  1:156706/1‑145 (MQ=255)
                                                                                    cGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggcgg                                                               <  7:139319/149‑1 (MQ=255)
                                                                                       tttCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTttct                                                                                                                         <  6:293986/89‑1 (MQ=255)
                                                                                       tttCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTttct                                                                                                                         >  5:293986/1‑89 (MQ=255)
                                                                                       tttCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGccc                                                            <  8:204871/149‑1 (MQ=255)
                                                                                           gCCGATCACCGTGCAGGCGTGCGGTGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGt                                                                                                                                                      >  2:146713/1‑57 (MQ=255)
                                                                                           gCCGATCACCGTGCAGGCGTGCGGTGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGt                                                                                                                                                      <  1:146713/57‑1 (MQ=255)
                                                                                           gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgt                                                        >  2:260012/1‑149 (MQ=255)
                                                                                                 cACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTcc                                                   >  1:239130/1‑148 (MQ=255)
                                                                                                  aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCa                                                 <  1:164795/149‑1 (MQ=255)
                                                                                                  aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCAGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCa                                                 <  8:196132/149‑1 (MQ=255)
                                                                                                    cGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggc                                                                 >  5:309932/1‑131 (MQ=255)
                                                                                                    cGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggc                                                                 <  6:309932/131‑1 (MQ=255)
                                                                                                     gTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTttct                                                                                                                         <  4:57645/75‑1 (MQ=255)
                                                                                                     gTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTttct                                                                                                                         >  3:57645/1‑75 (MQ=255)
                                                                                                        cAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGATTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGa                                            <  7:267345/148‑1 (MQ=255)
                                                                                                         aGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATc                                          >  6:365302/1‑149 (MQ=255)
                                                                                                            cGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTgg                                       <  4:250398/149‑1 (MQ=255)
                                                                                                                 ggcgCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATgggg                                                                   >  5:118713/1‑116 (MQ=255)
                                                                                                                 ggcgCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATgggg                                                                   <  6:118713/116‑1 (MQ=255)
                                                                                                                       aGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTggag                            >  7:338328/1‑149 (MQ=255)
                                                                                                                             gCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAggg                       <  1:260012/148‑1 (MQ=255)
                                                                                                                             gCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGt                      <  2:374258/149‑1 (MQ=255)
                                                                                                                               aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTgg                    <  8:248429/149‑1 (MQ=255)
                                                                                                                                     gCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTcc                                                   >  5:298404/1‑113 (MQ=255)
                                                                                                                                     gCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTcc                                                   <  6:298404/113‑1 (MQ=255)
                                                                                                                                     ccTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCAtg               >  1:12250/2‑149 (MQ=255)
                                                                                                                                       tCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTgg                    >  6:277869/1‑142 (MQ=255)
                                                                                                                                       tCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTgg                    <  5:277869/142‑1 (MQ=255)
                                                                                                                                          gCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATCGTggg          >  8:303861/1‑149 (MQ=255)
                                                                                                                                                  ggTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTgg                                       >  7:379538/1‑112 (MQ=255)
                                                                                                                                                  ggTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTgg                                       <  8:379538/112‑1 (MQ=255)
                                                                                                                                                    tCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGGTTg  >  2:203721/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGGTTG  >  NC_002947/1953261‑1953554

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: