Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 370,752 T→G 7.2% L78R (CTG→CGG) ‡ dgcA → dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947370,7520TG 7.2% 216.5 / 8.0 70L78R (CTG→CGG) ‡dgcAdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base G (3/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGC  >  NC_002947/370604‑370892
                                                                                                                                                    |                                                                                                                                            
gTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCt                                                                                                                                              <  1:81413/149‑1 (MQ=255)
   tACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTc                                                                                                                                           >  3:170353/1‑149 (MQ=255)
     cGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTcga                                                                                                                                         <  7:358472/149‑1 (MQ=255)
        cACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTcgaccg                                                                                                                                      >  3:373998/1‑149 (MQ=255)
        cACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTcgacc                                                                                                                                       >  4:170068/1‑148 (MQ=255)
          cTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTcgacc                                                                                                                                       <  6:182880/146‑1 (MQ=255)
          cTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTcgacc                                                                                                                                       >  5:182880/1‑146 (MQ=255)
                     aTGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATccc                                                                                                                         <  4:85963/149‑1 (MQ=255)
                          gaccgaccGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGc                                                                                                                       >  1:16764/1‑146 (MQ=255)
                          gaccgaccGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGc                                                                                                                       <  2:16764/146‑1 (MQ=255)
                          gaccgaccGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACt                                                                                                                    <  5:47046/149‑1 (MQ=255)
                            ccgaccgCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCACAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTc                                                                                                                  >  5:155610/1‑149 (MQ=255)
                             cgaccgCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTcc                                                                                                                 <  4:373998/149‑1 (MQ=255)
                                 cgCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCGACAGGTtgaccgatgc                                                                                                                                  >  3:134182/1‑125 (MQ=255)
                                 cgCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCGACAGGTtgaccgatgc                                                                                                                                  <  4:134182/128‑4 (MQ=255)
                                  gCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAAc                                                                                                            >  4:320761/1‑149 (MQ=255)
                                   cTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAAcc                                                                                                           >  6:238984/1‑149 (MQ=255)
                                   cTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAAcc                                                                                                           >  5:382082/1‑149 (MQ=255)
                                      tgtgAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTTGGCTTCGTTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCtt                                                                                                         <  1:234096/148‑1 (MQ=255)
                                       gtgAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCTGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTg                                                                                                        <  4:319397/148‑1 (MQ=255)
                                         gAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGcc                                                                                                      >  1:270352/1‑148 (MQ=255)
                                                tacGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTcgaccgac                                                                                                                                    <  4:44310/111‑1 (MQ=255)
                                                tacGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTcgaccgac                                                                                                                                    <  6:406713/111‑1 (MQ=255)
                                                     agagAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGcac                                                                                          >  2:179438/1‑148 (MQ=255)
                                                     agagAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGcac                                                                                          >  7:176976/1‑148 (MQ=255)
                                                      gagaAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACaa                                                                                        <  1:179438/149‑1 (MQ=255)
                                                      gagaAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACaa                                                                                        >  7:336742/1‑149 (MQ=255)
                                                          aGGCCAAGGGCGGCATCGGGCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCAt                                                                                    <  5:238984/149‑1 (MQ=255)
                                                          aGGCCAAGCGCGGCATCGTGCTGGCCTTCTGTGGCGGCTCGTCGGTGGTGTCCATCGACATCCCGGAGGAATGGTGGTCTTCGGTCTTCCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCAt                                                                                    <  7:256073/149‑1 (MQ=255)
                                                             ccAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATccc                                                                                                                         >  5:416466/1‑109 (MQ=255)
                                                             ccAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATccc                                                                                                                         <  6:416466/109‑1 (MQ=255)
                                                              cAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCAt                                                                                    <  7:314611/145‑1 (MQ=255)
                                                              cAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCAt                                                                                    >  8:314611/1‑145 (MQ=255)
                                                              cAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGgcg                                                                                >  2:367810/1‑149 (MQ=255)
                                                                    cggcATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGa                                                                          <  6:382082/149‑1 (MQ=255)
                                                                      gcATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTTGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATc                                                                        <  1:134259/149‑1 (MQ=255)
                                                                      gcATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGgc                                                                                 <  7:207532/140‑1 (MQ=255)
                                                                      gcATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGgc                                                                                 >  8:207532/1‑140 (MQ=255)
                                                                       cATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCa                                                                       >  6:226659/1‑149 (MQ=255)
                                                                           ggCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGa                                                                    <  2:270352/148‑1 (MQ=255)
                                                                                     tCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTAcc                                                         >  6:196561/1‑149 (MQ=255)
                                                                                         tgGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTAccc                                                        <  3:170068/146‑1 (MQ=255)
                                                                                           gcggcTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACAt                                                    <  4:170353/148‑1 (MQ=255)
                                                                                                       tggtggCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGtcgctcg                                       >  1:88431/1‑149 (MQ=255)
                                                                                                          tggCCATCGACAGCCCACAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCgcgc                                    <  6:155610/149‑1 (MQ=255)
                                                                                                            gCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGCCCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTg                                  >  7:282905/1‑149 (MQ=255)
                                                                                                              cATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGcac                                                                                          <  7:306633/91‑1 (MQ=255)
                                                                                                              cATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGcac                                                                                          >  8:306633/1‑91 (MQ=255)
                                                                                                                   aCAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGtcgctcgc                                      >  5:368640/1‑138 (MQ=255)
                                                                                                                   aCAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGtcgctcgc                                      <  6:368640/138‑1 (MQ=255)
                                                                                                                     aGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCtt                         <  5:226659/149‑1 (MQ=255)
                                                                                                                      gCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTc                        <  1:367810/149‑1 (MQ=255)
                                                                                                                       cccGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCCCATGGGCCGTCGCTCGCGCTGGGACGGCTTc                        >  7:49562/1‑148 (MQ=255)
                                                                                                                        ccGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATggg                                                 <  2:46807/122‑1 (MQ=255)
                                                                                                                        ccGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATggg                                                 >  1:46807/1‑122 (MQ=255)
                                                                                                                         cGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAAc                     <  8:336742/149‑1 (MQ=255)
                                                                                                                            aGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCtt                         >  5:95680/1‑142 (MQ=255)
                                                                                                                            aGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCtt                         <  6:95680/142‑1 (MQ=255)
                                                                                                                                aTGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGa              <  2:88431/149‑1 (MQ=255)
                                                                                                                                        ggtCGGTCGACAGGTTGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTg                   >  4:61530/3‑136 (MQ=255)
                                                                                                                                        ggtCGGTCGACAGGTTGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTg                   <  3:61530/134‑1 (MQ=255)
                                                                                                                                        cATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCtt                         >  4:5108/1‑130 (MQ=255)
                                                                                                                                        cATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCtt                         <  3:5108/130‑1 (MQ=255)
                                                                                                                                        cATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACt                    >  5:87560/1‑135 (MQ=255)
                                                                                                                                        cATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACt                    <  6:87560/135‑1 (MQ=255)
                                                                                                                                        cATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATg      >  3:361071/1‑149 (MQ=255)
                                                                                                                                        cATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGAATGGCGCCAAGAACATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATg      >  6:318764/1‑149 (MQ=255)
                                                                                                                                          tCGGTCGACAGGTTGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTc    >  5:345356/1‑149 (MQ=255)
                                                                                                                                             gTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGAt       >  6:118039/1‑143 (MQ=255)
                                                                                                                                             gTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGAt       <  5:118039/143‑1 (MQ=255)
                                                                                                                                             gTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGc  >  4:85696/1‑148 (MQ=255)
                                                                                                                                                    tGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCaa                                                                            >  4:120986/1‑67 (MQ=255)
                                                                                                                                                    tGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCaa                                                                            <  3:120986/67‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGCTCGTCGGTGGTGGCCATCGACAGCCCGCAGGAATGGTGGGCATCGGTCAACCTGTCGACCGACCGCATCATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGC  >  NC_002947/370604‑370892

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: