Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,250,528 T→A 6.9% I64F (ATC→TTC)  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,5280TA 6.9% 215.6 / 4.9 72I64F (ATC→TTC) nthDNA glycosylase/apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (3/2);  total (37/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 = 8.13e-01

CCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGGGGTTGTC  >  NC_002947/1250381‑1250672
                                                                                                                                                     |                                                                                                                                                
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                                                                                                       >  7:309364/1‑149 (MQ=255)
                                          gcTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGcc                                                                                                       >  1:372263/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                                                                                                      <  3:139126/128‑1 (MQ=255)
                                                                ttCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCt                                                                                                      >  4:139126/1‑128 (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                                         <  6:290865/148‑1 (MQ=255)
                                                                                                          tCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgg                                         <  2:45544/148‑1 (MQ=255)
                                                                                                           cTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGCTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgggg                                       >  6:250484/1‑149 (MQ=255)
                                                                                                            ttGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGt                                      >  1:33856/1‑149 (MQ=255)
                                                                                                             tGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgggg                                       <  5:280168/147‑1 (MQ=255)
                                                                                                             tGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgggg                                       >  6:280168/1‑147 (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)
                                                                                                                                                 tAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGGGGTTGTc  >  4:88515/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGGGGTTGTC  >  NC_002947/1250381‑1250672

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: