Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,887 G→T 6.4% A239D (GCC→GAC)  mutY ← A/G‑specific adenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,8870GT 6.4% 196.3 / 2.4 62A239D (GCC→GAC) mutYA/G‑specific adenine glycosylase
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base T (2/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTTGCTGCGCGTACACAGC  >  NC_002947/344742‑345033
                                                                                                                                                 |                                                                                                                                                  
ggcgtggcgGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGcc                                                                                                                                                  >  8:202395/1‑148 (MQ=255)
 gcgtggcgGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCaa                                                                                                                                                >  5:15118/1‑149 (MQ=255)
        ggcCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGGCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGc                                                                                                                                         <  5:188631/149‑1 (MQ=255)
          cccgccagGCGCAGGCCGTGCTGGGAGGCCAGGGCGGCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCCCAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCAt                                                                                                                                       <  6:203570/146‑1 (MQ=255)
             gccagGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCaa                                                                                                                                    <  1:252032/149‑1 (MQ=255)
              ccagGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAAt                                                                                                                                   <  6:15118/149‑1 (MQ=255)
                       ggCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGt                                                                                                                          >  1:182378/1‑149 (MQ=255)
                       ggCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGt                                                                                                                          >  1:40397/1‑149 (MQ=255)
                                 gTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAgcg                                                                                                                <  2:182378/149‑1 (MQ=255)
                                          gtcgtcGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAgcg                                                                                                                >  8:234936/1‑140 (MQ=255)
                                          gtcgtcGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAgcg                                                                                                                <  7:234936/140‑1 (MQ=255)
                                           tcccggAGCTGGACGACATCGCCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCgg                                                                                                      >  2:216220/6‑149 (MQ=255)
                                               cgAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTg                                                                                                  >  7:39701/1‑149 (MQ=255)
                                                gAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTg                                                                                                  <  3:195064/148‑1 (MQ=255)
                                                  gCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGc                                                                                               <  3:137022/149‑1 (MQ=255)
                                                      ggCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCgg                                                                                           <  5:139400/149‑1 (MQ=255)
                                                      ggCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCg                                                                                            <  2:40397/148‑1 (MQ=255)
                                                         gacatcgccCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCTATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGt                                                                                        <  1:216220/145‑1 (MQ=255)
                                                            gtcgtcCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTg                                                                                                  <  3:45101/136‑1 (MQ=255)
                                                            gtcgtcCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTg                                                                                                  >  4:45101/1‑136 (MQ=255)
                                                               gtcCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTa                                                                                       <  5:1129/144‑1 (MQ=255)
                                                               gtcCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTa                                                                                       >  6:1129/1‑144 (MQ=255)
                                                               gtcCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCggg                                                                                  >  8:234738/1‑149 (MQ=255)
                                                                tcCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGt                                                                                 >  7:232848/1‑149 (MQ=255)
                                                                   aGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTc                                                                              >  3:28293/1‑149 (MQ=255)
                                                                          ggAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCg                                                                       >  3:80458/1‑149 (MQ=255)
                                                                             ggCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTg                                                                     >  3:192617/1‑148 (MQ=255)
                                                                             ggCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGc                                                                    <  8:39701/149‑1 (MQ=255)
                                                                              gCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGTGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTg                                                                     >  2:183796/1‑147 (MQ=255)
                                                                              gCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATCGCGCCTTCGTGCTTGGCCAACAGTGGCATCAATGTGCGCCGATGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCa                                                                   >  8:89689/1‑149 (MQ=255)
                                                                                tCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTgc                                                                                                          >  3:249941/1‑108 (MQ=255)
                                                                                tCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTgc                                                                                                          <  4:249941/108‑1 (MQ=255)
                                                                                   acaAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTgc                                                              >  4:260847/1‑149 (MQ=255)
                                                                                   acaAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTgc                                                              >  1:264792/1‑149 (MQ=255)
                                                                                       accaccCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTATCCTCGCCGTGCAGGTGCGCtt                                                          >  1:120666/1‑149 (MQ=255)
                                                                                         caccCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCg                                                        <  7:202395/149‑1 (MQ=255)
                                                                                              cacaGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGTCCAACAGTGGCATCAATGTGCGCCGTTGCGTCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGcc                                                   >  2:275301/1‑149 (MQ=255)
                                                                                              cacaGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGTCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGcc                                                   >  2:271094/1‑149 (MQ=255)
                                                                                              cacaGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCg                                                        <  3:76920/144‑1 (MQ=255)
                                                                                              cacaGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCg                                                        >  4:76920/1‑144 (MQ=255)
                                                                                                caGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCcgc                                                 <  2:264792/149‑1 (MQ=255)
                                                                                                 aGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCcgcg                                                >  3:225138/1‑149 (MQ=255)
                                                                                                   gCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGTGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCt                                              <  1:183796/149‑1 (MQ=255)
                                                                                                        tGGAAGGGCGGCGGTAAAGCAGGATGGCGCGTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAAc                                         >  1:25505/1‑149 (MQ=255)
                                                                                                          gAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACgg                                       >  7:287021/1‑149 (MQ=255)
                                                                                                               gcggcgGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACggg                                      >  6:227943/1‑145 (MQ=255)
                                                                                                                 ggcggTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCcgcg                                                <  1:148078/133‑1 (MQ=255)
                                                                                                                 ggcggTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCcgcg                                                >  2:148078/1‑133 (MQ=255)
                                                                                                                    ggTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAgg                             >  8:87408/1‑149 (MQ=255)
                                                                                                                    ggTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAgg                             <  5:291488/149‑1 (MQ=255)
                                                                                                                          gCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTg                       >  6:255425/1‑149 (MQ=255)
                                                                                                                           cAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTgg                      <  7:87408/149‑1 (MQ=255)
                                                                                                                           cAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTgg                      <  8:287021/149‑1 (MQ=255)
                                                                                                                              gATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCt                   <  7:234738/149‑1 (MQ=255)
                                                                                                                                  gcgcCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTtgct               >  3:244540/1‑149 (MQ=255)
                                                                                                                                   cgcCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACggg                                      >  2:240516/1‑125 (MQ=255)
                                                                                                                                   cgcCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACggg                                      <  1:240516/125‑1 (MQ=255)
                                                                                                                                    gcagtCGTGGGTTGGCAACATTGGCATCTATGTGCGCCGTTGCGTGAGCGGCTTGCGCGGTTTGGACTCGGGGTAGCGGGTTTCCTCGCCGTGCTGGTTCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTtgctgc             <  2:132060/145‑1 (MQ=255)
                                                                                                                                    gcCTTCGTGGTTGTCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTtgctgc             <  1:271094/149‑1 (MQ=255)
                                                                                                                                    gcCTTCGTGGTTGTCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTtgctgc             <  1:275301/149‑1 (MQ=255)
                                                                                                                                           tGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTTGCTGCGCGTaca      >  8:13273/1‑149 (MQ=255)
                                                                                                                                            ggTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTTGCTGCGCGTaca      >  7:198468/1‑148 (MQ=255)
                                                                                                                                                ggCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTTGCTGCGCGTACACAGc  <  3:260847/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCTCCACAAACCACCCCACAGGCCGCTGGAAGGGCGGCGGTAAAGCAGGATGGCGCCTTCGTGGTTGGCCAACAGTGGCATCAATGTGCGCCGTTGCGGCAGCGCCTTGCGCGGTTTGGGCTCGGGGTAGCGGGTTTCCTCGCCGTGCAGGTGCGCTTCGCAGCCGCGCTGCAACGGGCAGATCAGGCAACTGGGCTTGCTGCGCGTACACAGC  >  NC_002947/344742‑345033

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: ‑