Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,250,523 A→G 7.0% Y65Y (TAT→TAC) ‡ nth ← DNA glycosylase/apyrimidinic (AP) lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,250,5230AG 7.0% 200.6 / 4.7 71Y65Y (TAT→TAC) ‡nthDNA glycosylase/apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base G (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TAGCGCCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGGGGTTGT  >  NC_002947/1250376‑1250671
                                                                                                                                                     |                                                                                                                                                    
tAGCGCCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCAt                                                                                                                                                     <  8:34022/149‑1 (MQ=255)
   cgcCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCa                                                                                                                                                      >  3:51892/1‑145 (MQ=255)
     ccTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATg                                                                                                                                                >  6:305677/1‑149 (MQ=255)
          aGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCGTCTATGACATTCTTGGCCTTGGTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCACCGCCCAGGGCATAGATGGCCTc                                                                                                                                           >  1:240531/1‑149 (MQ=255)
            cgcTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCgg                                                                                                                                         <  5:305677/149‑1 (MQ=255)
               ttCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCggggg                                                                                                                                      >  1:261767/1‑149 (MQ=255)
               ttCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCggggg                                                                                                                                      >  5:119668/1‑149 (MQ=255)
                 cgcgGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTTGCTGAGGCCATCTATGCCCTGGGCGTTGAAGGCCTcagcgaa                                                                                                                                     <  1:73295/148‑7 (MQ=37)
                 cgcgGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGgt                                                                                                                                     <  1:178142/148‑1 (MQ=255)
                   cgGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGtt                                                                                                                                  <  3:157345/149‑1 (MQ=255)
                    gggTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTg                                                                                                                                 <  5:49094/149‑1 (MQ=255)
                            ggTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACggg                                                                                                                         <  1:21094/149‑1 (MQ=255)
                                 cTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGaa                                                                                                                      >  3:42968/1‑147 (MQ=255)
                                        tCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCggg                                                                                                              <  8:180135/148‑1 (MQ=255)
                                         cGTGGCGCTCGAGCACCAGCCGTCACGCTTCTAGGACAGTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCa                                                                                                            <  8:318557/149‑1 (MQ=255)
                                          gTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAg                                                                                                           >  1:2662/1‑149 (MQ=255)
                                              cgcTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCGTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGc                                                                                                       >  3:347000/1‑149 (MQ=255)
                                               gcTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACa                                                                                                                    <  6:311951/135‑1 (MQ=255)
                                               gcTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACa                                                                                                                    >  5:311951/1‑135 (MQ=255)
                                               gcTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGcc                                                                                                      >  1:372263/1‑149 (MQ=255)
                                               gcTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGcc                                                                                                      >  7:309364/1‑149 (MQ=255)
                                                   gATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCTACGGGGAACAGGCGGGCAGTGGCCtgtt                                                                                                  >  7:322762/1‑146 (MQ=255)
                                                   gATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCtt                                                                                                    <  4:51892/147‑1 (MQ=255)
                                                    aTCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCttgtt                                                                                                 <  4:329318/149‑1 (MQ=255)
                                                       agcagcCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGa                                                                                               <  7:250719/148‑1 (MQ=255)
                                                           gcCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAACGCCAGATTGTCTTGATGTATTCGCTGAGGCGTTCAACGCCCAGGTCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCg                                                                                          <  2:240531/149‑1 (MQ=255)
                                                            ccGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGa                                                                                         <  2:2662/149‑1 (MQ=255)
                                                                   gCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGt                                                                                  >  5:295979/1‑149 (MQ=255)
                                                                     ttCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCt                                                                                                     >  4:139126/1‑128 (MQ=255)
                                                                     ttCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCt                                                                                                     <  3:139126/128‑1 (MQ=255)
                                                                       cTATGACATTCTTGGCCTTGCTGTTGTAAAGCGCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGa                                                                               >  1:45544/1‑149 (MQ=255)
                                                                         aTGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTg                                                                            >  5:338384/1‑149 (MQ=255)
                                                                          tGACATTCTTGGCCTTGCTGTTTTAAAGACCTATTGTCTTGATGTATTCTCTGCGGCCATCAACGCCCAGGGCATAGATGGCATCTGGGGTGTTGGCGCCGGGGAACAGTAGGGCAGTGGCATTCTTGACGAGGACGTCGGTGGATAgc                                                                           >  4:294579/1‑149 (MQ=255)
                                                                          tGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGAtt                                                                             >  4:330090/1‑147 (MQ=255)
                                                                          tGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGAtt                                                                             <  4:42968/147‑1 (MQ=255)
                                                                            aCATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTgcgc                                                                         <  3:330090/149‑1 (MQ=255)
                                                                            aCATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTgcg                                                                          <  5:45835/148‑1 (MQ=255)
                                                                               ttCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGa                                                                      <  8:229039/149‑1 (MQ=255)
                                                                                    ggCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACa                                                                    <  3:22896/146‑1 (MQ=255)
                                                                                     gCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAgtg                                                                 >  3:342980/1‑148 (MQ=255)
                                                                                          gCTGTTGTAAAGCCGGATTGTCTTGATGTATTCGCTGAGGGCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGTGGGGGAACAGGCTGGCAGTGGCCTTGTTGACGCCGAAGTCGGTGGATTGCGCGGACAGTGTCACg                                                            <  5:192375/148‑1 (MQ=255)
                                                                                           cTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGc                                                          <  3:87002/149‑1 (MQ=255)
                                                                                            tgttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCGACCGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGc                                                          >  1:72451/1‑149 (MQ=255)
                                                                                            tgttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCAAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGc                                                          >  5:268523/1‑148 (MQ=255)
                                                                                              ttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGa                                                        <  2:372263/148‑1 (MQ=255)
                                                                                              ttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGa                                                        <  4:342980/148‑1 (MQ=255)
                                                                                                 tAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACcag                                                    <  2:261767/149‑1 (MQ=255)
                                                                                                 tAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACcag                                                    <  6:119668/149‑1 (MQ=255)
                                                                                                          gATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCa                                             >  5:290865/1‑148 (MQ=255)
                                                                                                          gATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCa                                             >  8:17459/1‑148 (MQ=255)
                                                                                                           aTTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCaaa                                           >  3:68391/1‑149 (MQ=255)
                                                                                                               tCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGCTGGCCTCGTGGGTGTTGGCGACGGGGAACAGCCGGGCAGTGGCCTAGTTGACGCCGACATCGTTGGATTGCCCGGACAGTGTCACGGCGACCACCAGTTCAAACggg                                       >  4:222127/1‑149 (MQ=255)
                                                                                                               tCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgg                                        <  2:45544/148‑1 (MQ=255)
                                                                                                               tCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgg                                        <  6:290865/148‑1 (MQ=255)
                                                                                                                cTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGCTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgggg                                      >  6:250484/1‑149 (MQ=255)
                                                                                                                 ttGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGt                                     >  1:33856/1‑149 (MQ=255)
                                                                                                                  tGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgggg                                      >  6:280168/1‑147 (MQ=255)
                                                                                                                  tGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgggg                                      <  5:280168/147‑1 (MQ=255)
                                                                                                                    aTGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCgt                                  >  7:85068/1‑149 (MQ=255)
                                                                                                                     tGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCgtg                                 <  6:340644/149‑1 (MQ=255)
                                                                                                                           ccccgAGGCCATCTATGCCCTGGGCGTTGAAGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGc                            >  4:189283/5‑148 (MQ=38)
                                                                                                                            gCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCa                          <  6:295979/149‑1 (MQ=255)
                                                                                                                            cccgAGGCCATCTATGCCCTGGGCGTTGAAGGCCGCGGGGGTGTCGGCCACGGGGAACAGGAGTGCAGTGGCCTTGTTGCCGCCGACGTCGGGGGCTTGCGCGGCCAGTGTCACGGCGACCAGCAGTTCAACCGGGGTCGGGTAGGc                            >  3:54474/4‑147 (MQ=255)
                                                                                                                             ccgAGGCCATCTATGCCCTGGGCGTTGAAGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAg                         >  3:373056/3‑149 (MQ=38)
                                                                                                                                aGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGt                        <  2:200031/147‑1 (MQ=255)
                                                                                                                                     ttCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCt                 >  5:235254/1‑149 (MQ=255)
                                                                                                                                     ttCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCt                 >  5:327798/1‑149 (MQ=255)
                                                                                                                                         aaCGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCttt               <  6:268523/148‑1 (MQ=255)
                                                                                                                                          aCGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTgg             >  5:92398/1‑149 (MQ=255)
                                                                                                                                               ctgggCGTTGAAGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCgg        <  3:189283/147‑1 (MQ=255)
                                                                                                                                                aGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCggg       <  6:338384/149‑1 (MQ=255)
                                                                                                                                                  ggCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGGGGt     >  1:164824/1‑149 (MQ=255)
                                                                                                                                                     aTAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGGGGTTGt  >  2:346131/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGCGCCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGGGGTTGT  >  NC_002947/1250376‑1250671

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: