Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 781,067 T→A 5.6% Q260L (CAG→CTG)  glyA‑II ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947781,0670TA 5.6% 215.1 / 2.5 72Q260L (CAG→CTG) glyA‑IIserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (2/2);  total (35/37)
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.73e-01

CCACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTCGGGTACACGC  >  NC_002947/780921‑781214
                                                                                                                                                  |                                                                                                                                                     
ccACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGTCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTgg                                                                                                                                                     >  2:289456/1‑149 (MQ=255)
ccACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTgg                                                                                                                                                     <  4:76380/149‑1 (MQ=255)
    gTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGgcgc                                                                                                                                                 >  6:303061/1‑149 (MQ=255)
     tCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGcc                                                                                                                                                >  3:62096/1‑149 (MQ=255)
       gAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTTGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGccc                                                                                                                                               <  3:70448/148‑1 (MQ=255)
       gAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATAGCCTTGGCCGCGATGACGTGCTCCAGCGGGCGCCCAGGGCGgccc                                                                                                                                              <  8:162492/148‑2 (MQ=255)
       gAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATAGCCTTGGCCGCGATGACGTGCTCCAGCGGGCGCCCAGGGCGgccc                                                                                                                                              >  7:162492/1‑147 (MQ=37)
         aaCCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTTAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTTGGCGCCCggg                                                                                                                                           <  2:173813/149‑1 (MQ=255)
         aaCCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCggg                                                                                                                                           <  6:308568/149‑1 (MQ=255)
         aaCCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCggg                                                                                                                                           <  3:208137/149‑1 (MQ=255)
               gcTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGa                                                                                                                                       <  8:102548/147‑1 (MQ=255)
               gcTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGa                                                                                                                                       >  7:102548/1‑147 (MQ=255)
                  cGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGAcggcg                                                                                                                                  <  3:17551/149‑1 (MQ=255)
                    aTGAACACGCTGGCCATGGCCTGGGCTTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGa                                                                                                                                <  8:260360/149‑1 (MQ=255)
                                  cATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTtcttct                                                                                                                  >  7:66356/1‑149 (MQ=255)
                                  cATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTtcttct                                                                                                                  >  3:362781/1‑149 (MQ=255)
                                      gCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGAt                                                                                                              <  8:157556/149‑1 (MQ=255)
                                        cTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGt                                                                                                            <  3:36152/149‑1 (MQ=255)
                                            gCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGc                                                                                                        <  6:347328/149‑1 (MQ=255)
                                                 cttcACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCt                                                                                                                                                       <  1:123697/98‑1 (MQ=255)
                                                 cttcACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCt                                                                                                                                                       >  2:123697/1‑98 (MQ=255)
                                                       gACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTg                                                                                                    <  6:251953/142‑1 (MQ=255)
                                                       gACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTg                                                                                                    >  5:251953/1‑142 (MQ=255)
                                                         ccTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGc                                                                                            >  8:219177/1‑148 (MQ=255)
                                                             ctgctgGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGa                                                                                       <  4:163931/149‑1 (MQ=255)
                                                                ctgGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGAt                                                                                      >  2:144881/1‑147 (MQ=255)
                                                                     aaGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCGCCCAggg                                                                                                                                                    >  3:92642/1‑80 (MQ=255)
                                                                     aaGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCGCCCAggg                                                                                                                                                    <  4:92642/80‑1 (MQ=255)
                                                                      aGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGgccg                                                                               <  4:283180/148‑1 (MQ=255)
                                                                        ccTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCAAGCGGGCCGCCCGGGGCGCCCGGGAAGCCGGCGGAGGTCAGCTTCTTCTCGATGTCGGAGTTGGCACGGGCCAGGATCATGccgccg                                                                            >  6:360286/1‑149 (MQ=255)
                                                                          ttGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCa                                                                                    <  1:183351/139‑1 (MQ=255)
                                                                          ttGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCa                                                                                    >  2:183351/1‑139 (MQ=255)
                                                                           tGAACTCTGGCTGCAGGGCTTCCTTGAAGAAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCACCTTCTTCTCGATGTCGGCGTTGGCACGGGCCTGGATCAGGCCTCcgcgc                                                                         >  7:122114/1‑149 (MQ=255)
                                                                            gAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCcgcgc                                                                         <  1:144881/148‑1 (MQ=255)
                                                                            gAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGA‑‑‑CGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGc                                                                     >  8:254231/1‑149 (MQ=255)
                                                                                tctGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCCCGGGCACGTATCAGGCCCCCCCCCGGGcg                                                                    >  7:236562/1‑148 (MQ=255)
                                                                                    gCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGc                                                                                                        <  8:109496/109‑1 (MQ=255)
                                                                                    gCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGc                                                                                                        >  7:109496/1‑109 (MQ=255)
                                                                                    gCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAg                                                                                   >  5:7999/1‑130 (MQ=255)
                                                                                    gCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAg                                                                                   <  6:7999/130‑1 (MQ=255)
                                                                                        cAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAgg                                                            <  4:62096/148‑1 (MQ=255)
                                                                                        cAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGA‑‑‑CGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGt                                                          <  7:254231/147‑1 (MQ=255)
                                                                                         aGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGt                                                          >  6:370470/1‑149 (MQ=255)
                                                                                             cTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTtg                                                      >  8:290171/1‑149 (MQ=255)
                                                                                              ttCCTTGAAGCAGATGGCCTTGGCCGGGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTtgt                                                     <  4:362781/149‑1 (MQ=255)
                                                                                              ttCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTtgt                                                     >  3:191238/1‑149 (MQ=255)
                                                                                                 cTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCACGCGGGCCGCGCAGGGTCTTgtggg                                                  >  4:334858/1‑149 (MQ=255)
                                                                                                  ttGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAgg                                                                                  >  5:42914/1‑117 (MQ=255)
                                                                                                  ttGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAgg                                                                                  <  6:42914/117‑1 (MQ=255)
                                                                                                   tGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTgtgggtg                                                >  1:180571/1‑149 (MQ=255)
                                                                                                    gAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTgtgggtgg                                               >  3:376133/1‑149 (MQ=255)
                                                                                                    gAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTgtgggtg                                                <  5:303061/148‑1 (MQ=255)
                                                                                                     aaGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGggtggt                                              >  3:325545/1‑149 (MQ=255)
                                                                                                       gCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGggtggtgg                                            <  8:66356/149‑1 (MQ=255)
                                                                                                        cAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGggtggtggt                                           <  7:219177/149‑1 (MQ=255)
                                                                                                         aGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGCCGGCGGAGGGCAGCTTCTTCTCTATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGggttgtggtg                                          <  5:360286/149‑1 (MQ=25)
                                                                                                         aGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCACGCGGGCCGCGCAGGGTCTTGTGggtggtggtg                                          <  3:334858/149‑1 (MQ=255)
                                                                                                               ccTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGAc                                     >  5:339388/1‑148 (MQ=255)
                                                                                                                   ggCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACCGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGt                                >  6:83766/1‑149 (MQ=255)
                                                                                                                    gCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGGGGTGGTGACCACGTc                               <  8:181317/149‑1 (MQ=255)
                                                                                                                                 gCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGaa                    >  1:146855/1‑147 (MQ=255)
                                                                                                                                   tCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTtgt                                                     >  2:1744/1‑112 (MQ=255)
                                                                                                                                   tCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTtgt                                                     <  1:1744/112‑1 (MQ=255)
                                                                                                                                    ccAGCGGGCCGCCCTGGGCGCCCGGGGAAGACGGCGGAGTTCAGGTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCCGGCCGCCGCGCCGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAAc                        >  1:196626/1‑141 (MQ=21)
                                                                                                                                    ccAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAAc                        <  2:196626/140‑1 (MQ=255)
                                                                                                                                     cAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAAc                        <  4:48882/139‑1 (MQ=255)
                                                                                                                                     cAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAAc                        >  3:48882/1‑139 (MQ=255)
                                                                                                                                       gCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTc            <  2:146855/149‑1 (MQ=255)
                                                                                                                                       gCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGGTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTc            <  4:297172/149‑1 (MQ=255)
                                                                                                                                         gggCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTCgg          >  7:131420/1‑149 (MQ=255)
                                                                                                                                           gccgccCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGCTGTCTTGTGTGTGGTGGTTGTGACCACGTCGGCGAACGGAACAGGGTTCggg         <  1:289456/149‑1 (MQ=25)
                                                                                                                                             cgccCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTCGGGTa       >  3:120315/1‑148 (MQ=255)
                                                                                                                                                 cTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGTTGGTGACAAAGTCGGCGAAAGGCACAGGGTTAGGGTACACGc  >  6:268424/1‑149 (MQ=25)
                                                                                                                                                 cTGGGCGACAGGGAAGACTGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTCGGGTACACGc  <  8:122114/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTCGGGTACACGC  >  NC_002947/780921‑781214

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: