Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 781,061 Δ1 bp 5.5% coding (785/1254 nt) glyA‑II ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947781,0610C. 5.5% 238.5 / 6.2 73coding (785/1254 nt)glyA‑IIserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base . (2/2);  total (36/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.86e-01

GAAACCACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTCGGG  >  NC_002947/780917‑781207
                                                                                                                                                |                                                                                                                                                    
gAAACCACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGgccgcc                                                                                                                                                  >  4:70448/1‑149 (MQ=255)
 aaaCCACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGccc                                                                                                                                                 >  3:163931/1‑149 (MQ=255)
    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)
           gAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGccc                                                                                                                                                 >  5:91406/1‑139 (MQ=255)
           gAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGccc                                                                                                                                                 <  6:91406/139‑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                                                                                                                                    <  3:208137/149‑1 (MQ=255)
             aaCCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCggg                                                                                                                                    <  6:308568/149‑1 (MQ=255)
                   gcTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGa                                                                                                                                >  7:102548/1‑147 (MQ=255)
                   gcTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGa                                                                                                                                <  8:102548/147‑1 (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                 >  3:48882/1‑139 (MQ=255)
                                                                                                                                         cAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAAc                 <  4:48882/139‑1 (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)
                                                                                                                                                |                                                                                                                                                    
GAAACCACGTCGAAACCGCGCTCGATGAACACGCTGGCCATGGCCTGGGCGTTCTTCACGACCTGCTGCTGGTAAGCCTTGAACTCTGGCTGCAGGGCTTCCTTGAAGCAGATGGCCTTGGCCGCGATGACGTGCTCCAGCGGGCCGCCCTGGGCGCCCGGGAAGACGGCGGAGTTCAGCTTCTTCTCGATGTCGGCGTTGGCACGGGCCAGGATCAGGCCGCCGCGCGGGCCGCGCAGGGTCTTGTGGGTGGTGGTGGTGACCACGTCGGCGAACGGAACAGGGTTCGGG  >  NC_002947/780917‑781207

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: