Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,096,069 T→A 7.7% L134* (TTG→TAG) ‡ PP_5607 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,096,0690TA 7.7% 181.0 / 3.5 62L134* (TTG→TAG) ‡PP_5607membrane protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base A (3/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGGCCTG  >  NC_002947/4095928‑4096214
                                                                                                                                             |                                                                                                                                                  
cGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAgc                                                                                                                                             <  5:425931/149‑1 (MQ=255)
  ggCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAgcgc                                                                                                                                           <  1:258416/149‑1 (MQ=255)
   gCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcgag                                                                                                                                          >  3:7979/1‑147 (MQ=255)
   gCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcga                                                                                                                                           <  4:7979/148‑2 (MQ=255)
    ccTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGccc                                                                                                                                         >  3:390453/1‑149 (MQ=255)
        gtgtTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATg                                                                                                                                     >  7:396485/1‑149 (MQ=255)
         tgtTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATgg                                                                                                                                    <  1:265267/149‑1 (MQ=255)
             gctgcATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCg                                                                                                                                 <  4:389103/148‑1 (MQ=255)
                 catggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTgg                                                                                                                              <  5:188475/147‑1 (MQ=255)
                 catggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCt                                                                                                                            <  1:394376/149‑1 (MQ=255)
                   tggcatggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAg                                                                                                                          >  5:363527/1‑149 (MQ=255)
                        tggGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAAc                                                                                                                     <  4:323063/149‑1 (MQ=255)
                                  gggCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGcc                                                                                                            <  8:97371/148‑1 (MQ=255)
                                      aCCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTc                                                                                                       >  6:437798/1‑149 (MQ=255)
                                        cGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccc                                                                                                     >  6:164589/1‑149 (MQ=255)
                                        cGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccc                                                                                                     >  4:114670/1‑149 (MQ=255)
                                        cGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccc                                                                                                     >  6:121616/1‑149 (MQ=255)
                                          gCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTccccc                                                                                                   >  2:167494/1‑149 (MQ=255)
                                                 ccccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcga                                                                                                                                           >  6:381189/1‑101 (MQ=255)
                                                 ccccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGGCCAAGGTCAgcga                                                                                                                                           <  5:381189/102‑2 (MQ=255)
                                                   ccTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGgcg                                                                                          >  3:58524/1‑149 (MQ=255)
                                                        gACCAGGGGTGTGCACAGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTACTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGtt                                                                                     <  1:298918/149‑1 (MQ=255)
                                                           cAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTc                                                                                  <  5:437798/149‑1 (MQ=255)
                                                               ggTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCTTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTg                                                                             >  1:292208/1‑149 (MQ=255)
                                                                gtgtGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTgg                                                                            <  6:378978/149‑1 (MQ=255)
                                                                gtgtGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTgg                                                                            <  7:360570/149‑1 (MQ=255)
                                                                     cACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGc                                                                       <  8:41812/149‑1 (MQ=255)
                                                                      aCCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCg                                                                      >  8:210966/1‑149 (MQ=255)
                                                                       ccGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGc                                                                       <  6:193798/147‑1 (MQ=255)
                                                                       ccGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGc                                                                       >  5:193798/1‑147 (MQ=255)
                                                                          gTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCgtggt                                                                  >  1:121646/1‑149 (MQ=255)
                                                                             cAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGtt                                                               <  6:363527/149‑1 (MQ=255)
                                                                                ccATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCgtgg                                                                   >  6:181349/1‑142 (MQ=255)
                                                                                ccATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCgtgg                                                                   <  5:181349/142‑1 (MQ=255)
                                                                                 cATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGc                                                            >  1:57023/1‑148 (MQ=255)
                                                                                       ccTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGatcatc                                                     <  7:210966/149‑1 (MQ=255)
                                                                                       ccTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGatc                                                        <  5:121616/146‑1 (MQ=255)
                                                                                        cTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGAtcatca                                                    >  5:275333/1‑149 (MQ=255)
                                                                                                 aCTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGt                                           <  3:114670/149‑1 (MQ=255)
                                                                                                  cTGCATAACCTGTTTGGTCTGCGGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGtt                                          >  7:271667/1‑149 (MQ=255)
                                                                                                     cATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGAt                                       <  2:426907/149‑1 (MQ=255)
                                                                                                      aTAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATg                                      <  2:292208/149‑1 (MQ=255)
                                                                                                         aCCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGc                                   <  2:57023/149‑1 (MQ=255)
                                                                                                            tGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGa                                                          >  8:164548/1‑123 (MQ=255)
                                                                                                            tGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGa                                                          <  7:164548/123‑1 (MQ=255)
                                                                                                                 ggTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCACTTAGGTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcgg                            >  1:141234/1‑149 (MQ=255)
                                                                                                                 ggTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcgg                            <  1:156210/148‑1 (MQ=255)
                                                                                                                    ctgctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggcg                          >  5:182078/1‑147 (MQ=255)
                                                                                                                     tgctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTGTGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAgcggcg                          >  4:72788/1‑146 (MQ=255)
                                                                                                                       ctgATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCaga                      <  5:95363/148‑1 (MQ=255)
                                                                                                                          aTGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGt                    <  4:58524/147‑1 (MQ=255)
                                                                                                                          aTGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTg                   <  6:182078/148‑1 (MQ=255)
                                                                                                                            gTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGc                 >  2:14172/1‑148 (MQ=255)
                                                                                                                             ttCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGc                 >  2:170081/1‑147 (MQ=255)
                                                                                                                             ttCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCa               <  4:390453/149‑1 (MQ=255)
                                                                                                                              atagggcGCTGGCCAAGGTCAGCGACCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGAAGAGTGGCCAg              >  6:367050/7‑149 (MQ=255)
                                                                                                                               ccTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAg              <  6:275333/148‑1 (MQ=255)
                                                                                                                               ccTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAg              <  5:164589/148‑1 (MQ=255)
                                                                                                                                     gCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTc        <  2:121646/148‑1 (MQ=255)
                                                                                                                                      cTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCgg      <  1:167494/149‑1 (MQ=255)
                                                                                                                                           cTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATTGGCCAAGCGGCGCAGAGTGGCCTGCTGTTCGGCCTg  >  3:272994/1‑148 (MQ=255)
                                                                                                                                           cTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGGCCTg  >  7:342428/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CGGGCCTGGTGTTGCTGCATGGCATGGGCAGTTTGGGCACCGGCTGGTTCCCCTGCGACCAGGGGTGTGCACCGGTACAGCCATCTACCTCGCAGCAACTGCATAACCTGTTTGGTCTGCTGATGTTCCTCTCGCTGACCTTGGCCAGCGCCCTATGGGCGTGGCTAGGCAACCGCATTGCCGGCTCCCCCGCCCTGGCGCTGTTTTCCCTGGCTTGCGTGGTGTTGGCGATCATCACCGTCGCGTTGATGGGCCAAGCGGCGCAGAGTGGCCAGCTGTTCGGCCTG  >  NC_002947/4095928‑4096214

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: