Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 819,246 Δ1 bp 6.3% coding (400/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,2460G. 6.3% 202.6 / 7.1 63coding (400/627 nt)PP_0705DNA‑3‑methyladenine glycosylase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base . (2/2);  total (31/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 = 2.89e-01

GCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCCGA  >  NC_002947/819101‑819393
                                                                                                                                                   |                                                                                                                                                    
gCTGGCGGTATCCCCGGAAACCCTGCATGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAAc                                                                                                                                                   >  1:11166/1‑149 (MQ=255)
 cTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAAc                                                                                                                                                   <  8:159384/148‑1 (MQ=255)
    gCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATCAACGCCt                                                                                                                                               <  4:202326/149‑1 (MQ=255)
    gCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATCAACGCCt                                                                                                                                               <  4:311638/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                                                                                                                                <  3:280130/149‑1 (MQ=255)
                   aCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgg                                                                                                                                <  4:177051/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                                                                                                                     >  3:56784/1‑64 (MQ=255)
                                                                                                                   gCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggt                                                                                                                     <  4:56784/64‑1 (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)
                                                                                                                                               gATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGgcgc     >  4:149237/1‑149 (MQ=255)
                                                                                                                                                 tCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGcc    >  7:192155/1‑148 (MQ=255)
                                                                                                                                                  cGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGAAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCCGa  >  8:181388/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCCGA  >  NC_002947/819101‑819393

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: