Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,277,325 T→G 8.3% Q15H (CAA→CAC PP_1116 ← resolvase family site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,277,3250TG 8.3% 220.4 / 8.9 72Q15H (CAA→CACPP_1116resolvase family site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (2/4);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCAAA  >  NC_002947/1277180‑1277473
                                                                                                                                                 |                                                                                                                                                    
gagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGa                                                                                                                                                   <  8:39391/149‑1 (MQ=255)
gagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGa                                                                                                                                                   >  8:123886/1‑149 (MQ=255)
      ttGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTg                                                                                                                                             >  1:181278/1‑149 (MQ=255)
      ttGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCAAGGTGActctg                                                                                                                                              <  1:81148/148‑2 (MQ=255)
         cACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCAAGGTGActctggcc                                                                                                                                           <  5:90431/148‑5 (MQ=255)
          aCCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTgg                                                                                                                                            >  3:203090/1‑146 (MQ=255)
           ccTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGa                                                                                                                                        <  4:203090/149‑1 (MQ=255)
           ccTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGa                                                                                                                                        <  7:123886/149‑1 (MQ=255)
             ttCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGAt                                                                                                                                       >  3:157317/1‑148 (MQ=255)
             ttCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCAAGGTGActctg                                                                                                                                              >  8:127850/1‑140 (MQ=255)
             ttCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCAAGGTGActctg                                                                                                                                              <  7:127850/141‑2 (MQ=255)
               cTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATAc                                                                                                                                     <  4:157317/148‑1 (MQ=255)
                  gACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCg                                                                                                                                  >  1:88266/1‑148 (MQ=255)
                    cGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAAt                                                                                                                               <  6:237848/149‑1 (MQ=255)
                     gATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATg                                                                                                                              >  2:265131/1‑149 (MQ=255)
                       tCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGt                                                                                                                             >  4:205111/1‑148 (MQ=255)
                         ttCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAg                                                                                                                          >  1:326145/1‑149 (MQ=255)
                          tCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACAACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAg                                                                                                                          <  4:275839/148‑1 (MQ=255)
                           cAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGc                                                                                                                         >  4:175609/1‑148 (MQ=255)
                              gagcgagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGcaac                                                                                                                      >  1:232763/1‑148 (MQ=255)
                                  gagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTgg                                                                                                                                            <  4:73048/122‑1 (MQ=255)
                                  gagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTgg                                                                                                                                            >  3:73048/1‑122 (MQ=255)
                                   agcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGCGTCCAGCCCTAGGCCAGAGTCCCCTTGACGCTTGGTAGATACTCGAATGTAAGCAACAACGgtt                                                                                                               >  4:101464/1‑149 (MQ=255)
                                    gcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGttt                                                                                                               <  2:88266/149‑1 (MQ=255)
                                    gcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGttt                                                                                                               >  4:33685/1‑149 (MQ=255)
                                                  acaTCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGc                                                                                                 >  4:186312/1‑149 (MQ=255)
                                                   caTCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGAc                                                                                                                                                  >  3:198307/1‑99 (MQ=255)
                                                   caTCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGAc                                                                                                                                                  <  4:198307/99‑1 (MQ=255)
                                                        acaAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCt                                                                                           >  4:204360/1‑149 (MQ=255)
                                                           aaTTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGTATGTAATCACGCTGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTgtgt                                                                                       <  3:101464/149‑1 (MQ=255)
                                                               cACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTAt                                                                                     <  2:181278/148‑1 (MQ=255)
                                                                aCCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGTCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTaa                                                                                                                           <  5:245783/109‑1 (MQ=255)
                                                                aCCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTaa                                                                                                                           >  6:245783/1‑109 (MQ=255)
                                                                  cAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGActct                                                                                                                                               <  4:6776/87‑1 (MQ=255)
                                                                  cAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGActct                                                                                                                                               >  3:6776/1‑87 (MQ=255)
                                                                  cAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAgg                                                                                 >  1:101041/1‑149 (MQ=255)
                                                                  cAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAgg                                                                                 <  1:265131/149‑1 (MQ=255)
                                                                  cAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAgg                                                                                 <  3:186312/149‑1 (MQ=255)
                                                                    aCGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGt                                                                                                                                           >  7:2159/1‑89 (MQ=255)
                                                                    aCGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGt                                                                                                                                           <  8:2159/89‑1 (MQ=255)
                                                                            ccAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTTGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTAt                                                                        <  2:156769/148‑1 (MQ=255)
                                                                            ccAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTAt                                                                        <  3:175609/148‑1 (MQ=255)
                                                                                   ttCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGcc                                                                <  2:232763/149‑1 (MQ=255)
                                                                                      gCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCAt                                                              >  7:241576/1‑148 (MQ=255)
                                                                                        ttGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGACGCCAATTTATCAACTGCCATGAg                                                           >  5:157423/1‑149 (MQ=255)
                                                                                          gggCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTg                                                                                                              <  7:260830/96‑1 (MQ=255)
                                                                                          gggCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTg                                                                                                              >  8:260830/1‑96 (MQ=255)
                                                                                           ggCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTa                                                         <  3:205111/148‑1 (MQ=255)
                                                                                               gCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGt                                                    <  8:241576/149‑1 (MQ=255)
                                                                                                    ttGGATGTAATCACGCTGGGCTTCCATACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAAtt                                               <  3:33685/149‑1 (MQ=255)
                                                                                                    ttGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAAtt                                               >  6:165559/1‑149 (MQ=255)
                                                                                                    ttGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAAtt                                               <  2:101041/149‑1 (MQ=255)
                                                                                                        aTGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTg                                              >  1:12246/1‑146 (MQ=255)
                                                                                                        aTGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTg                                              <  2:12246/146‑1 (MQ=255)
                                                                                                            aaTCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGAc                                       >  7:218254/1‑149 (MQ=255)
                                                                                                            aaTCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGAc                                       >  7:154942/1‑149 (MQ=255)
                                                                                                            aaTCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGAc                                       >  4:337237/1‑149 (MQ=255)
                                                                                                             aTCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGa                                                                                   <  4:246980/104‑1 (MQ=255)
                                                                                                             aTCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGa                                                                                   >  3:246980/1‑104 (MQ=255)
                                                                                                                      ggCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTa                               <  5:165559/147‑1 (MQ=255)
                                                                                                                       gCTTCCAGACCTAGGCCAGATTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGAGCAGACGTTCGTTAGt                             <  3:204360/148‑1 (MQ=255)
                                                                                                                             aGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATg                      >  5:149262/1‑149 (MQ=255)
                                                                                                                                  tAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAAt                 >  5:149900/1‑149 (MQ=255)
                                                                                                                                   aGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGAcgttcgt                                 <  6:193206/132‑1 (MQ=255)
                                                                                                                                   aGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGAcgttcgt                                 >  5:193206/1‑132 (MQ=255)
                                                                                                                                      ccAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTa             <  6:149262/149‑1 (MQ=255)
                                                                                                                                       aagagTCAAGGTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTc               >  6:88345/2‑146 (MQ=255)
                                                                                                                                        agagTCAAGGTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGa           <  5:88345/149‑1 (MQ=255)
                                                                                                                                         gagTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGAt          >  3:269232/1‑149 (MQ=255)
                                                                                                                                              aCCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAGATGGAAATTCTAGATTTcc      <  8:218254/148‑1 (MQ=255)
                                                                                                                                              aCCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTcc      <  3:337237/148‑1 (MQ=255)
                                                                                                                                              aCCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTcc      >  8:182550/1‑148 (MQ=255)
                                                                                                                                                 ttGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCaaa  >  3:206924/1‑149 (MQ=255)
                                                                                                                                                 ttGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCaaa  >  2:336004/1‑149 (MQ=255)
                                                                                                                                                 ttGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATAAACTGCCAtac                                                            >  2:47911/1‑89 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCAAA  >  NC_002947/1277180‑1277473

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: