Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,602,002 T→A 6.2% intergenic (+15/‑152) PP_1404 → / → PP_5469 hypothetical protein/SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,602,0020TA 6.2% 199.7 / 2.6 64intergenic (+15/‑152)PP_1404/PP_5469hypothetical protein/SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base A (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

GGTGCGCCAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCC  >  NC_002947/1601858‑1602147
                                                                                                                                                 |                                                                                                                                                 
ggTGCGCCAGGGCCGTACCTTGCAGGGGTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCt                                                                                                                                               <  4:173458/149‑1 (MQ=255)
ggTGCGCCAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCt                                                                                                                                               <  6:157921/149‑1 (MQ=255)
ggTGCGCCAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTACCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCt                                                                                                                                               <  7:65091/149‑1 (MQ=255)
ggTGCGCCAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTACCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCt                                                                                                                                               <  7:60664/149‑1 (MQ=255)
   gcgcCAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTg                                                                                                                                              >  1:229729/1‑148 (MQ=255)
     gcCAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCAAGGCTGGacg                                                                                                                                          <  1:204497/149‑3 (MQ=255)
      ccAGGGCCGTACCTTACAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCAAGGCTGGacgc                                                                                                                                         >  5:132992/1‑146 (MQ=255)
       cAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTc                                                                                                                                        <  4:242795/149‑1 (MQ=255)
       cAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTc                                                                                                                                        <  5:16567/149‑1 (MQ=255)
          ggCCGTACCTTTCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCAAGGCTGGacgctcct                                                                                                                                     <  6:132992/149‑8 (MQ=255)
            ccGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCtgtg                                                                                                                                   >  4:15727/1‑149 (MQ=255)
                  cTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGcaca                                                                                                                             <  1:14452/149‑1 (MQ=255)
                       aGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCtt                                                                                                                        >  5:123885/1‑149 (MQ=255)
                            tttGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCtt                                                                                                                   >  4:58095/1‑149 (MQ=255)
                                tCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGACTTGCCTTGGa                                                                                                                <  6:285869/148‑1 (MQ=255)
                                  gAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGaaa                                                                                                              >  4:261975/1‑148 (MQ=255)
                                    gTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTGGCGTTGGAAACAg                                                                                                           >  7:123083/1‑149 (MQ=255)
                                       cATCGAAGCGTGTGGCGAGCTGCTAAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGc                                                                                                         >  7:203022/1‑148 (MQ=255)
                                        aTCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTg                                                                                                                  >  7:101097/1‑138 (MQ=255)
                                        aTCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTg                                                                                                                  <  8:101097/138‑1 (MQ=255)
                                        aTCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGa                                                                                                       >  6:194257/1‑149 (MQ=255)
                                        aTCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGa                                                                                                       >  8:45455/1‑149 (MQ=255)
                                          cGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCt                                                                                                                                               <  7:8460/107‑1 (MQ=255)
                                          cGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCt                                                                                                                                               >  8:8460/1‑107 (MQ=255)
                                            aaGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAgg                                                                                                   <  2:229729/149‑1 (MQ=255)
                                            aaGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAgg                                                                                                   <  3:58095/149‑1 (MQ=255)
                                                        gAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCAc                                                                                        <  8:203022/149‑1 (MQ=255)
                                                            tgctTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGc                                                                                    <  3:15727/149‑1 (MQ=255)
                                                             gctTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCt                                                                                   >  5:181550/1‑149 (MQ=255)
                                                              ctTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGc                                                                                    <  2:86599/147‑1 (MQ=255)
                                                                  gCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGGGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGc                                                                               <  8:155568/148‑1 (MQ=255)
                                                                   ccAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTcc                                                                                                                                       <  1:166119/91‑1 (MQ=255)
                                                                   ccAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTcc                                                                                                                                       >  2:166119/1‑91 (MQ=255)
                                                                   ccAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCg                                                                             <  1:209364/149‑1 (MQ=255)
                                                                   ccAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCg                                                                             <  3:261975/149‑1 (MQ=255)
                                                                           tggcgccACCCCATGCACGCAAGGAGGTGCGCAGCCGCTTGGTGATTCTCGACTGACGGAGGGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGTCAGGGTGTCTGGCACCGGCTTTGCCGGTGTTcgc                                                                     <  8:123083/146‑1 (MQ=255)
                                                                           tGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTcgc                                                                     >  3:9647/1‑149 (MQ=255)
                                                                           tGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTcgc                                                                     >  6:332385/1‑149 (MQ=255)
                                                                             ccgccAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCgg                                                                   <  4:9647/149‑1 (MQ=255)
                                                                             ccgccAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCgg                                                                   <  7:45455/149‑1 (MQ=255)
                                                                                    cccATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGcc                                                            >  3:257037/1‑149 (MQ=255)
                                                                                        tGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCt                                                        >  7:163164/1‑149 (MQ=255)
                                                                                           acacAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCaa                                                               <  3:248571/136‑1 (MQ=255)
                                                                                           acacAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCaa                                                               >  4:248571/4‑139 (MQ=255)
                                                                                            cGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCca                                                    >  8:56971/1‑149 (MQ=255)
                                                                                             gCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCca                                                    >  3:126211/1‑148 (MQ=255)
                                                                                              cAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCcaca                                                  <  5:332385/149‑1 (MQ=255)
                                                                                                 cGAGTTGCCCAACCGCTTGGTGATTCACGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAggg                                               >  3:233668/1‑149 (MQ=255)
                                                                                                         ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTg                                                                                <  4:254496/108‑1 (MQ=255)
                                                                                                         ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTg                                                                                >  3:254496/1‑108 (MQ=255)
                                                                                                          cAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                      <  4:126211/149‑1 (MQ=255)
                                                                                                             ccGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCAc                                   >  4:262341/1‑149 (MQ=255)
                                                                                                                       aTTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATa                          <  8:163164/148‑1 (MQ=255)
                                                                                                                        ttCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAgg                        <  7:56971/149‑1 (MQ=255)
                                                                                                                        ttCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAgg                        <  4:233668/149‑1 (MQ=255)
                                                                                                                           tcGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                      >  6:1500/1‑132 (MQ=255)
                                                                                                                           tcGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                      <  5:1500/132‑1 (MQ=255)
                                                                                                                             gACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGtta                   >  6:228496/1‑148 (MQ=255)
                                                                                                                             gACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTc                   >  2:256108/1‑149 (MQ=255)
                                                                                                                                          tCCAGACTTGCTCGATTTCCTCTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGAAAGCCCCATCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGc      >  2:113709/1‑149 (MQ=38)
                                                                                                                                           ccAGCCATGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGcc     >  2:279524/1‑149 (MQ=255)
                                                                                                                                              gCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGccccc  >  3:131801/1‑149 (MQ=255)
                                                                                                                                                cTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                      >  6:195264/1‑111 (MQ=255)
                                                                                                                                                cTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                      <  5:195264/111‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTGCGCCAGGGCCGTACCTTGCAGGGCTTTGTCGAGTGCATCGAAGCGTGTGGCGAGCTGCTTAGCCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCC  >  NC_002947/1601858‑1602147

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: