Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 819,241:1 +C 6.7% coding (395/627 nt) PP_0705 → DNA‑3‑methyladenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947819,2411.C 6.7% 234.2 / 9.9 60coding (395/627 nt)PP_0705DNA‑3‑methyladenine glycosylase
Reads supporting (aligned to +/- strand):  ref base . (26/30);  new base C (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

GCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGA  >  NC_002947/819101‑819382
                                                                                                                                              |                                                                                                                                              
gCTGGCGGTATCCCCGGAAACCCTGCATGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAAc                                                                                                                                        >  1:11166/1‑149 (MQ=255)
 cTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAAc                                                                                                                                        <  8:159384/148‑1 (MQ=255)
    gCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATCAACGCCt                                                                                                                                    <  4:311638/149‑1 (MQ=255)
    gCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATCAACGCCt                                                                                                                                    <  4:202326/149‑1 (MQ=255)
      ggTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCtgg                                                                                                                                  >  8:246150/1‑149 (MQ=255)
          tCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGgc                                                                                                                              <  5:105774/149‑1 (MQ=255)
          tCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGCCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCtggtgg                                                                                                                               <  1:72480/148‑1 (MQ=255)
             ccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCt                                                                                                                           >  3:41407/1‑149 (MQ=255)
             ccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCt                                                                                                                           >  8:24046/1‑149 (MQ=255)
             ccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCt                                                                                                                           <  2:351143/149‑1 (MQ=255)
                  aaCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTg                                                                                                                      <  4:41407/149‑1 (MQ=255)
                   aCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgg                                                                                                                     <  4:177051/149‑1 (MQ=255)
                   aCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgg                                                                                                                     <  3:280130/149‑1 (MQ=255)
                     ccTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgggg                                                                                                                   >  4:193824/1‑149 (MQ=255)
                      cTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGGCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGt                                                                                                                  <  3:59378/149‑1 (MQ=255)
                                       ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                  <  8:127012/148‑1 (MQ=255)
                                           gCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAt                                                                                             >  7:69716/1‑149 (MQ=255)
                                             gAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgc                                                                                           >  2:268756/1‑149 (MQ=255)
                                                cAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc                                                                                        <  3:264947/149‑1 (MQ=255)
                                                  aGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATgctgct                                                                                       >  6:358223/1‑148 (MQ=255)
                                                   gCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGa                                                                                     >  7:129997/1‑149 (MQ=255)
                                                       gTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATc                                                                                   <  6:50270/147‑1 (MQ=255)
                                                       gTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGCACCGACGCGTGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc                                                                                        <  5:288053/142‑1 (MQ=255)
                                                       gTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGCACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc                                                                                        >  6:288053/1‑142 (MQ=255)
                                                         gACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCtata                                                                               <  7:96564/149‑1 (MQ=255)
                                                          aCGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCtata                                                                               <  2:11166/148‑1 (MQ=255)
                                                          aCGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCtata                                                                               <  1:90369/148‑1 (MQ=255)
                                                          aCGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCtata                                                                               <  8:129997/148‑1 (MQ=255)
                                                           cGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGc                                                                             <  3:193824/149‑1 (MQ=255)
                                                           cGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGc                                                                             >  2:364403/1‑149 (MQ=255)
                                                              tCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGa                                                                                     <  3:56874/138‑1 (MQ=255)
                                                              tCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGa                                                                                     >  4:56874/1‑138 (MQ=255)
                                                                cAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTgg                                                                         <  7:24046/148‑1 (MQ=255)
                                                                       ttGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCt                                                                <  8:69716/149‑1 (MQ=255)
                                                                       ttGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCt                                                                >  1:308928/1‑149 (MQ=255)
                                                                           tCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACa                                                            >  2:376687/1‑149 (MQ=255)
                                                                                      tGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTg                                                 >  4:340721/1‑149 (MQ=255)
                                                                                       ggAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTg                                                 >  1:283397/1‑149 (MQ=255)
                                                                                             ccTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGAtt                                          >  6:105097/1‑149 (MQ=255)
                                                                                              cTGCTACAGACGCGGGAGCAGGCGTTAACCATGCCTGACGAGGCGTTGATAGAACCCATGGTGGCGCTGAGTGGGGTATGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg                                                                        >  3:213479/1‑119 (MQ=255)
                                                                                               tGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGAtttt                                        >  2:274620/1‑149 (MQ=255)
                                                                                                 gTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTg                                       <  5:358223/148‑1 (MQ=255)
                                                                                                      gACGCGGGAGCAGGCGATAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTTGGGTAGGGCGGTGGACGGAGGAGATGCTGCAGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTgc                                 >  2:336118/1‑149 (MQ=255)
                                                                                                        cgcgGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTgcgc                               <  2:308928/149‑1 (MQ=255)
                                                                                                          cgGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGa                             <  7:246150/149‑1 (MQ=255)
                                                                                                             gAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGCGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAgg                           <  1:268756/149‑1 (MQ=255)
                                                                                                                   gCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggt                                                                                                          <  4:56784/64‑1 (MQ=255)
                                                                                                                   gCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggt                                                                                                          >  3:56784/1‑64 (MQ=255)
                                                                                                                      ttAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATcggcg                 >  4:263340/1‑149 (MQ=255)
                                                                                                                          ccATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCt              >  3:350442/1‑148 (MQ=255)
                                                                                                                            aTGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCaa           <  5:190625/149‑1 (MQ=255)
                                                                                                                               gCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAgg         <  5:105097/148‑1 (MQ=255)
                                                                                                                                   accaGGCGTTCGATCAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg    >  1:325889/4‑149 (MQ=255)
                                                                                                                                   aCGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg    <  3:340721/149‑1 (MQ=255)
                                                                                                                                   aCGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg    <  1:364403/149‑1 (MQ=255)
                                                                                                                                    ccaGGCGTTCGATCAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg    >  1:260471/3‑148 (MQ=255)
                                                                                                                                     gAGGCGTTGATCGACCGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGga  <  2:43446/149‑1 (MQ=255)
                                                                                                                                     gAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg    >  2:211440/1‑147 (MQ=255)
                                                                                                                                      aGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg    >  7:247246/1‑146 (MQ=255)
                                                                                                                                      aGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg    >  3:355159/1‑146 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGA  >  NC_002947/819101‑819382

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: