Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,173,658 G→A 5.8% G193S (GGC→AGC)  birA → bifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,173,6580GA 5.8% 221.5 / 2.1 69G193S (GGC→AGC) birAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base A (2/2);  total (37/32)
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.43e-01

GCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGC  >  NC_000913/4173514‑4173805
                                                                                                                                                |                                                                                                                                                   
gcggcgATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCaa                                                                                                                                                 <  7:248588/149‑1 (MQ=255)
      aTTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTgg                                                                                                                                           >  4:38276/1‑149 (MQ=255)
          gTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGAt                                                                                                                                       >  1:274269/1‑149 (MQ=255)
            ttAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATgc                                                                                                                                     >  7:202505/1‑149 (MQ=255)
                 tCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATgcggcg                                                                                                                                 >  7:368043/1‑148 (MQ=255)
                   tGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCa                                                                                                                               >  8:496/1‑148 (MQ=255)
                        aTCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAg                                                                                                                          <  3:38276/148‑1 (MQ=255)
                           ggTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCa                                                                                                                       <  7:496/148‑1 (MQ=255)
                           ggTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCAt                                                                                                                      >  7:306983/1‑149 (MQ=255)
                                  tGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTAGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGcc                                                                                                               >  2:37846/1‑149 (MQ=255)
                                    aTGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCgg                                                                                                             <  8:36058/149‑1 (MQ=255)
                                    aTGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCgg                                                                                                             <  8:202505/149‑1 (MQ=255)
                                       gCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGa                                                                                                           <  7:295826/148‑1 (MQ=255)
                                       gCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGAt                                                                                                          >  8:308080/1‑149 (MQ=255)
                                         ggAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTAGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATc                                                                                                         >  8:4959/1‑148 (MQ=255)
                                             gTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATGCTGGGGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACAt                                                                                                    <  1:273768/149‑1 (MQ=255)
                                                ttACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATgg                                                                                                  >  8:189423/1‑148 (MQ=255)
                                                ttACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATgg                                                                                                  >  5:149460/1‑148 (MQ=255)
                                                   cGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCaa                                                                                               <  2:43669/148‑1 (MQ=255)
                                                   cGCAAGCTGGGTGCAGATAAAGTTAGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAAt                                                                                              >  2:113127/1‑149 (MQ=255)
                                                    gCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATg                                                                                             >  7:129971/1‑149 (MQ=255)
                                                      aaGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTCGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATgcg                                                                                           >  2:12056/1‑149 (MQ=255)
                                                                cAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGaa                                                                                 <  7:61957/149‑1 (MQ=255)
                                                                  gATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAaga                                                                               <  8:129971/149‑1 (MQ=255)
                                                                    tAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAagaga                                                                             <  8:306983/149‑1 (MQ=255)
                                                                     aaaGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAagagag                                                                            <  7:189423/149‑1 (MQ=255)
                                                                        gTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAgtgt                                                                         >  1:11632/1‑149 (MQ=255)
                                                                        gTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCAGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAgtgt                                                                         >  7:294401/1‑149 (MQ=255)
                                                                          tCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCg                                                                       >  3:160154/1‑149 (MQ=255)
                                                                              gtTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCAGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCgtgt                                                                                     <  7:51725/131‑1 (MQ=255)
                                                                              gtTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCAGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCgtgt                                                                                     >  8:51725/1‑131 (MQ=255)
                                                                              gtTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCGGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGGTAATAGTCATTGGGGGCGGGATCAACATGGCTATGCGCCGTGTTGATGGGAGTGTGGGTaa                                                                   >  8:234432/1‑149 (MQ=255)
                                                                               ttAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAAt                                                                  >  3:80021/1‑149 (MQ=255)
                                                                                  aaTGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAg                                                               >  1:211307/1‑149 (MQ=255)
                                                                                  aaTGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAg                                                               <  4:160154/149‑1 (MQ=255)
                                                                                  aaTGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAg                                                               <  8:368043/149‑1 (MQ=255)
                                                                                   aTGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTAGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAgg                                                              <  1:37846/149‑1 (MQ=255)
                                                                                    tGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAggg                                                             <  1:266738/149‑1 (MQ=255)
                                                                                           aTGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGAt                                                       <  7:308080/148‑1 (MQ=255)
                                                                                                ctctATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCt                                                 >  4:213730/1‑149 (MQ=255)
                                                                                                ctctATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCt                                                 >  7:271894/1‑149 (MQ=255)
                                                                                                     tCTGCTTGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCATAACTGGCCTTGCGTCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGc                                               <  4:272812/146‑1 (MQ=255)
                                                                                                     tCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAg                                             >  8:124598/1‑148 (MQ=255)
                                                                                                            gATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCa                                                                >  2:193565/1‑122 (MQ=255)
                                                                                                            gATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCa                                                                <  1:193565/122‑1 (MQ=255)
                                                                                                            gATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCggg                                     <  1:113127/149‑1 (MQ=255)
                                                                                                              tCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGa                                   <  7:124598/149‑1 (MQ=255)
                                                                                                                  aaGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCaa                               <  7:109423/149‑1 (MQ=255)
                                                                                                                   aGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCa                                                                >  7:143445/1‑115 (MQ=255)
                                                                                                                   aGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCa                                                                <  8:143445/115‑1 (MQ=255)
                                                                                                                    gCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAtc                             >  1:135295/1‑149 (MQ=255)
                                                                                                                      tGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTTGATCACGCTGCAGGAAGCGGGGATCAAtctc                           >  3:274904/1‑149 (MQ=255)
                                                                                                                        gcaggcaTTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGa                         <  7:234432/149‑1 (MQ=255)
                                                                                                                        gcaggcaTTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGAGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGa                         >  8:140901/1‑149 (MQ=255)
                                                                                                                          aggcaTTCTGGTGGAGCTGACTAGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCATGGGGGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATc                       <  7:4959/149‑1 (MQ=255)
                                                                                                                           ggcaTTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCg                      >  5:308881/1‑149 (MQ=255)
                                                                                                                               ttCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTaa                   >  4:302942/1‑148 (MQ=255)
                                                                                                                               ttCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAAt                  <  3:302942/149‑1 (MQ=255)
                                                                                                                                 cTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCgg                                                                                                             <  5:88878/56‑1 (MQ=255)
                                                                                                                                 cTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCgg                                                                                                             >  6:88878/1‑56 (MQ=255)
                                                                                                                                  tggtggAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAg                                                               >  8:18112/1‑101 (MQ=255)
                                                                                                                                  tggtggAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAg                                                               <  7:18112/101‑1 (MQ=255)
                                                                                                                                    gtggAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGtt             >  4:8662/1‑149 (MQ=255)
                                                                                                                                    gtggAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGtt             >  6:233714/1‑149 (MQ=255)
                                                                                                                                             aCTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCa     <  8:271894/148‑1 (MQ=255)
                                                                                                                                             aCTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCAt    >  3:323340/1‑149 (MQ=255)
                                                                                                                                             aCTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGGGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGcc      <  8:294401/147‑1 (MQ=255)
                                                                                                                                              cTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATg   >  2:192438/1‑149 (MQ=255)
                                                                                                                                               tGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGc  <  2:135295/149‑1 (MQ=255)
                                                                                                                                                ggCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGc  >  2:335035/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGC  >  NC_000913/4173514‑4173805

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: