Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 354,883 G→T 8.2% Q97K (CAG→AAG)  cbcW ← choline/betaine/carnitine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947354,8830GT 8.2% 188.5 / 6.9 61Q97K (CAG→AAG) cbcWcholine/betaine/carnitine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base T (3/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GGTAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGA  >  NC_002947/354737‑355029
                                                                                                                                                  |                                                                                                                                                   
ggtAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTgcg                                                                                                                                                   >  5:269366/1‑149 (MQ=255)
     ggATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCCTGGCGCAGGTCACCTTCGCCacg                                                                                                                                              <  8:204486/149‑3 (MQ=255)
        tCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGtc                                                                                                                                           >  2:121773/1‑149 (MQ=255)
          aGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGtctc                                                                                                                                         <  6:269366/149‑1 (MQ=255)
          aGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGtctc                                                                                                                                         <  8:16281/149‑1 (MQ=255)
          aGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGtct                                                                                                                                          <  6:293853/148‑1 (MQ=255)
          aGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTCACCTGCGCCAGGGtctc                                                                                                                                         >  8:79170/1‑149 (MQ=255)
             tAGTCGAAGGGGGGCACCGTCTGCATCAGGTCGATCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGACCAGTGGCACGCCGATCAACACACATACCACCGTGGCGAAGGTGACATGCGCCAGGGTCTCCa                                                                                                                                       <  8:175481/148‑1 (MQ=255)
               gACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑Tg                                                                                                                                    <  5:57570/148‑1 (MQ=255)
                  gAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGACCAGTGGCACGCCGATCACCTCACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGttt                                                                                                                                >  5:309343/1‑149 (MQ=255)
                       tgtgCACCGTGTGCATCAGGTAGAGCACCGGGCGCATGGCGGTATAGAAGATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATAGCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGc                                                                                                                           <  6:288407/146‑1 (MQ=255)
                         ggCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCa                                                                                                                         <  6:183110/149‑1 (MQ=255)
                          gCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAg                                                                                                                        >  3:173022/1‑149 (MQ=255)
                           cACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAg                                                                                                                        >  3:129858/1‑148 (MQ=255)
                           cACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAg                                                                                                                        >  3:135072/1‑148 (MQ=255)
                                        aGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAg                                                                                                          >  4:93224/1‑149 (MQ=255)
                                        aGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAg                                                                                                          <  3:93224/149‑1 (MQ=255)
                                         ggTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAgg                                                                                                         <  6:139024/149‑1 (MQ=255)
                                         ggTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAgg                                                                                                         <  1:271689/149‑1 (MQ=255)
                                             gAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                              <  6:270253/109‑1 (MQ=255)
                                             gAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                              >  5:270253/1‑109 (MQ=255)
                                                cACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATcagca                                                                                                  <  1:121773/149‑1 (MQ=255)
                                                      gcgcATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAAc                                                                                             >  1:358662/1‑148 (MQ=255)
                                                          aTGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAg                                                                                        >  2:230826/1‑149 (MQ=255)
                                                                     gAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATg                                                                             >  1:168859/1‑149 (MQ=255)
                                                                     gAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATg                                                                             >  1:174329/1‑149 (MQ=255)
                                                                     gAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATg                                                                             >  3:149705/1‑149 (MQ=255)
                                                                       aCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACGTCACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                           <  6:309343/149‑1 (MQ=255)
                                                                          tCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGt                                                                         >  5:152205/1‑148 (MQ=255)
                                                                           cGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATcag                                                                                                    >  3:256243/1‑120 (MQ=255)
                                                                           cGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATcag                                                                                                    <  4:256243/120‑1 (MQ=255)
                                                                            ggCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAg                                                                      <  4:173022/149‑1 (MQ=255)
                                                                                   gcgcGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCAc                                                               >  5:99464/1‑149 (MQ=255)
                                                                                    cgcgGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACt                                                              <  2:174329/149‑1 (MQ=255)
                                                                                    cgcgGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACt                                                              <  2:168859/149‑1 (MQ=255)
                                                                                     gcggcgAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCgg                                                                          <  7:287464/136‑1 (MQ=255)
                                                                                     gcggcgAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCgg                                                                          >  8:287464/1‑136 (MQ=255)
                                                                                          gAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCt                                                        <  6:88611/149‑1 (MQ=255)
                                                                                            ggATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCAc                                                               <  7:199980/140‑1 (MQ=255)
                                                                                            ggATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCAc                                                               >  8:199980/1‑140 (MQ=255)
                                                                                             gATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGa                                                     <  2:358662/149‑1 (MQ=255)
                                                                                                    aGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATg                                              >  4:328176/1‑149 (MQ=255)
                                                                                                          aCGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAg                                        >  3:336492/1‑149 (MQ=255)
                                                                                                               gATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCAGTGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGggcg                                    >  1:89269/1‑149 (MQ=255)
                                                                                                                aTCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGggcggc                                  <  7:79170/149‑1 (MQ=255)
                                                                                                                 tCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGGTGAAATGCGCCAGGgcggcg                                 <  4:135072/149‑1 (MQ=255)
                                                                                                                 tCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGGTGAAATGCGCCAGGgcggcg                                 <  4:129858/149‑1 (MQ=255)
                                                                                                                   accacACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGaa                               <  1:230826/149‑1 (MQ=255)
                                                                                                                        acaTACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAggg                          <  3:114351/149‑1 (MQ=255)
                                                                                                                         caTACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGc                         <  6:99464/149‑1 (MQ=255)
                                                                                                                                 cctGGCGCAGGTCACCTTCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAgaa                 >  3:333267/3‑149 (MQ=255)
                                                                                                                                 cctGGCGCAGGTCACCTTCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAgaa                 >  3:306009/3‑149 (MQ=255)
                                                                                                                                  ctGGCGCAGGTCACCTTCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAgaac                >  2:179449/2‑149 (MQ=255)
                                                                                                                                   tGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAAt                              >  7:242933/1‑134 (MQ=255)
                                                                                                                                   tGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAAt                              <  8:242933/134‑1 (MQ=255)
                                                                                                                                       gAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGaacgaacg           >  2:337363/1‑149 (MQ=255)
                                                                                                                                         aGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACggg         <  4:149705/149‑1 (MQ=255)
                                                                                                                                             gACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTg      >  2:2330/1‑148 (MQ=255)
                                                                                                                                              aCCTTCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGCAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTga     <  4:306009/148‑1 (MQ=255)
                                                                                                                                                cTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTgaag   >  6:256626/1‑148 (MQ=255)
                                                                                                                                                cTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTgaag   >  5:150214/1‑148 (MQ=255)
                                                                                                                                                cTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACACAGGTTGAGGATCAGCAGGAACCACAGGGCGAAAAATGCGGTCAGTGCCCACATGCGCTGGATGAAATGCGACAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTcaag   >  5:209920/1‑148 (MQ=255)
                                                                                                                                                  gcgcCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCAACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTgaaga  >  4:335821/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGTAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCA‑TGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGA  >  NC_002947/354737‑355029

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