Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 221,954 G→A 6.2% L242L (CTG→CTA)  PP_0169 → TauD/TfdA family dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947221,9540GA 6.2% 195.8 / 2.0 64L242L (CTG→CTA) PP_0169TauD/TfdA family dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/2);  total (33/31)
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.03e-01

CCGAGCAACTGGCGCAGGTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGGCGCAAGCTGTACCGC  >  NC_002947/221812‑222098
                                                                                                                                              |                                                                                                                                                 
ccGAGCAACTGGCGCAGGTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTgcagc                                                                                                                                               <  3:225050/147‑1 (MQ=255)
  gAGCAACTGGCGCAGGTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCt                                                                                                                                           >  5:169584/1‑149 (MQ=255)
   aGCAACTGGCGCAGGTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTg                                                                                                                                          >  6:141239/1‑149 (MQ=255)
        cTGGCGCAGGTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑c                                                                                                                                    >  4:188542/1‑149 (MQ=255)
           gcgcAGGTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑ccc                                                                                                                                  <  2:265846/148‑1 (MQ=255)
             gcAGGTAAAGGAGGTCATTCACCAGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCt                                                                                                                                 >  3:265221/1‑147 (MQ=255)
                 gTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGt                                                                                                                           >  2:261779/1‑149 (MQ=255)
                          gTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAg                                                                                                                    >  2:102510/1‑147 (MQ=255)
                          gTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCa                                                                                                                  <  5:334670/149‑1 (MQ=255)
                          gTCATTCACCCGGTAGTCCGCACCAACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCa                                                                                                                  <  6:322882/149‑1 (MQ=255)
                            cATTCACCCGGTTGTCCGCACCGACGCGGAGAACGGCCGGAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCGTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGc                                                                                                                         <  6:328506/140‑1 (MQ=255)
                            cATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGc                                                                                                                         >  5:328506/1‑140 (MQ=255)
                             aTTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGaa                                                                                                               >  5:200676/1‑149 (MQ=255)
                              ttCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTACAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGaa                                                                                                               >  7:302502/1‑148 (MQ=255)
                              ttCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTACAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGaa                                                                                                               >  1:37426/1‑148 (MQ=255)
                                cACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACa                                                                                                             <  3:204881/148‑1 (MQ=255)
                                     ggTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATc                                                                                                       <  1:334273/149‑1 (MQ=255)
                                           ccGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCAGCATTGTCTGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCATCAGCTGTATG‑CCCTTAGCGTGCTGGAGCATAACATTTATCGccac                                                                                                  <  7:29989/148‑1 (MQ=255)
                                             gCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAg                                                                                               >  2:181319/1‑149 (MQ=255)
                                               acccaaccGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTg                                                                                             <  4:160616/149‑1 (MQ=255)
                                                     ccGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTg                                                                                                                       >  2:209510/1‑117 (MQ=255)
                                                     ccGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTg                                                                                                                       <  1:209510/117‑1 (MQ=255)
                                                     ccGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGcc                                                                                       <  1:261779/149‑1 (MQ=255)
                                                      cGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTAt                                                                                                        >  4:75045/1‑131 (MQ=255)
                                                      cGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTAt                                                                                                        <  3:75045/131‑1 (MQ=255)
                                                      cGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGccc                                                                                      >  5:262576/1‑149 (MQ=255)
                                                       ggAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGccc                                                                                      >  3:36849/1‑148 (MQ=255)
                                                        gagaACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTa                                                                                    <  7:209300/149‑1 (MQ=255)
                                                            aCGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGAc                                                                                <  6:200676/149‑1 (MQ=255)
                                                               gCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTg                                                                             <  3:15872/149‑1 (MQ=255)
                                                                      ggCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCt                                                                      <  4:115151/149‑1 (MQ=255)
                                                                      ggCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCt                                                                      <  6:169584/149‑1 (MQ=255)
                                                                         ccTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTACAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTggg                                                                   <  2:37426/149‑1 (MQ=255)
                                                                         ccTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTACAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTggg                                                                   <  8:302502/149‑1 (MQ=255)
                                                                             ttCGTCAGTGAACGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATGCCCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGa                                                                  <  1:102510/147‑1 (MQ=255)
                                                                                 tCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGt                                                           <  4:36849/149‑1 (MQ=255)
                                                                                   aGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTc                                                         >  6:310038/1‑149 (MQ=255)
                                                                                          gCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATc                                                  >  2:271906/1‑149 (MQ=255)
                                                                                          gCTTCACAACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATc                                                  >  2:78990/1‑149 (MQ=255)
                                                                                            ttCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGCCAACCGTTCGCTGATCCa                                                >  5:109916/1‑149 (MQ=255)
                                                                                              caccacCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACt                                              >  6:84953/1‑149 (MQ=255)
                                                                                               accaccCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACtt                                             >  5:65986/1‑149 (MQ=255)
                                                                                                ccaccCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTg                                            >  3:103395/1‑149 (MQ=255)
                                                                                                  accCGCATTGTCGGCCTGCCTGATGAGGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGc                                          >  1:9046/1‑149 (MQ=255)
                                                                                                     cGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGATAACCGTTCGCTGATCCACTTGGCCAc                                       <  1:181319/149‑1 (MQ=255)
                                                                                                      gCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCAcc                                      >  8:185031/1‑149 (MQ=255)
                                                                                                          tGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCt                                  >  6:48550/1‑149 (MQ=255)
                                                                                                           gTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCt                                  <  6:65986/148‑1 (MQ=255)
                                                                                                               gcctgcctGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTgg                                                                                            <  1:200101/86‑1 (MQ=255)
                                                                                                               gcctgcctGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTgg                                                                                            >  2:200101/1‑86 (MQ=255)
                                                                                                                 ctgcctgATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCcgc                           <  3:188542/149‑1 (MQ=255)
                                                                                                                    cctgATGATGAAAGCCGGGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCa                        <  4:265221/149‑1 (MQ=255)
                                                                                                                       gatgatgaAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCt                     >  6:36698/1‑149 (MQ=255)
                                                                                                                        atgatgaAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTg                    >  5:251356/1‑149 (MQ=255)
                                                                                                                         tgatgaAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTgc                   >  3:102723/1‑149 (MQ=255)
                                                                                                                            tgaAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACtt                                             >  2:285208/1‑120 (MQ=255)
                                                                                                                            tgaAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACtt                                             <  1:285208/120‑1 (MQ=255)
                                                                                                                            tgaAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTgcggc                <  5:36698/149‑1 (MQ=255)
                                                                                                                              aaaGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTgcggcg               >  4:15945/1‑148 (MQ=255)
                                                                                                                              aaaGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTgcggcg               >  6:187268/1‑148 (MQ=255)
                                                                                                                               aaGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGgcgc              >  6:292419/1‑148 (MQ=255)
                                                                                                                                 gCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGGCGCAAg           <  5:310038/149‑1 (MQ=255)
                                                                                                                                 gCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGGCGCAAg           <  6:251356/149‑1 (MQ=255)
                                                                                                                                 gCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGGCGCAAg           <  5:141239/149‑1 (MQ=255)
                                                                                                                                       aCGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGGCGCAAGCTGTAc     <  5:84953/149‑1 (MQ=255)
                                                                                                                                          tgctgcAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGGCGCAAGCTGTACCGc  <  4:102723/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CCGAGCAACTGGCGCAGGTAAAGGAGGTCATTCACCCGGTAGTCCGCACCCACCCGGAGAACGGCCGCAAGGCCCTGTTCGTCAGTGAAGGCTTCACCACCCGCATTGTCGGCCTGCCTGATGATGAAAGCCGCGACGTGCTGCAGCAGCTGTATG‑CCCTTAGCGTGCTGGAGCAGAACATTTATCGCCACCAGTGGCAGCCCTATGACCTGGTGTTCTGGGACAACCGTTCGCTGATCCACTTGGCCACCGGCTGCCCCGCGCACCTGCGGCGCAAGCTGTACCGC  >  NC_002947/221812‑222098

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