Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,917,910 T→C 6.0% intergenic (+97/+135) barA → / ← gudD hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,917,9100TC 6.0% 191.1 / 3.0 67intergenic (+97/+135)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base T (23/40);  new base C (2/2);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATG  >  NC_000913/2917766‑2918057
                                                                                                                                                |                                                                                                                                                   
gACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                                 >  1:974014/1‑149 (MQ=255)
gACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                                 <  2:537874/149‑1 (MQ=255)
gACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                                 >  1:650462/1‑149 (MQ=255)
 aCGAGATGGATAATGGCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                                                                <  2:887104/149‑1 (MQ=255)
  cGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                                                               <  1:365612/149‑1 (MQ=255)
   gagaTGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                                                              <  1:537145/149‑1 (MQ=255)
   gagaTGGATAATGTCGCGCGCGAAGCCAGCAAAATACTCGCGTAATGGATGTCTAGTCGAGGAGCAATCGGCTTTTTATGACGCTGGGTTCAGATTGGTGACAAAGTGCCTGGGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                                                              <  1:38571/149‑1 (MQ=38)
     gaTGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                                                            >  2:25829/1‑149 (MQ=255)
     gaTGGATAATGTCGAGCGCGAAGCCATCAAAATTCTCGGGGAATGGATGTCTTGTCGAGTAGCAATCGGCGTTTTATGCCGCTGGATTCAGAAGGGCGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                                                            <  1:768322/149‑1 (MQ=37)
       tGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                                                          <  1:720187/149‑1 (MQ=255)
        ggATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                                                                         >  1:404188/1‑149 (MQ=255)
        ggATAAGGTCGCGCGCGAAGCCAGCATACTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGGGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                                                                         <  2:950097/149‑1 (MQ=255)
                              aGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAAt                                                                                                                   <  2:226205/149‑1 (MQ=255)
                               gCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAAtt                                                                                                                  >  1:766217/1‑149 (MQ=255)
                                cAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTa                                                                                                                 >  1:644337/1‑149 (MQ=255)
                                   aaTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaa                                                                                                              >  2:50287/1‑149 (MQ=255)
                                    aTTCTCGTGTAATGGATGTCTTGTCGAGTAGCAATCGGATTTTTATGACGCTGGTTTCAGATTGGTGACAAAGAGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaaa                                                                                                             <  1:296019/149‑1 (MQ=255)
                                    aTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaaa                                                                                                             <  2:316162/149‑1 (MQ=255)
                                       ctcGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatg                                                                                                          <  2:93303/149‑1 (MQ=255)
                                        tcGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatga                                                                                                         <  1:546184/149‑1 (MQ=255)
                                         cGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatgat                                                                                                        <  2:622031/149‑1 (MQ=255)
                                         cGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatgat                                                                                                        <  2:209953/149‑1 (MQ=255)
                                          gggTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatgatg                                                                                                       <  1:422882/149‑1 (MQ=255)
                                            gTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGtt                                                                                                     >  1:813465/1‑148 (MQ=255)
                                              aaTGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGa                                                                                                   <  2:810293/149‑1 (MQ=255)
                                                 ggATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCt                                                                                                <  2:306786/149‑1 (MQ=255)
                                                 ggATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCt                                                                                                <  1:202242/149‑1 (MQ=255)
                                                   aTGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGa                                                                                              <  1:800161/149‑1 (MQ=255)
                                                        ttGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACAt                                                                                         <  2:650462/149‑1 (MQ=255)
                                                         tGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                                        <  2:626051/149‑1 (MQ=255)
                                                          ggctaGTAGCAAGCGGCTTGTTATGCCGCTGGTTTCAGCTTGGTGACAACGTGCCTGTGGTTAGGCCGGATGCGGCGTTAACGCCGTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAg                                                                                       <  1:870776/145‑1 (MQ=38)
                                                          gTCGAGTAGCAATAGTCTTTTTATGAAGCTTGTTTCAGATTGGTGACAAAGTGCCTGTTTTTATTCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAg                                                                                       <  1:109005/149‑1 (MQ=38)
                                                             gAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGcac                                                                                    >  2:646964/1‑149 (MQ=255)
                                                                tAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATCTAGACCTGATACATAGCacaac                                                                                 <  2:979141/149‑1 (MQ=255)
                                                                  gCAATCTGCTTTTTATGACGCAGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAg                                                                               <  2:659263/149‑1 (MQ=255)
                                                                   cAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAgg                                                                              <  2:541759/149‑1 (MQ=255)
                                                                          cTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCATACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTAtttt                                                                       <  2:63393/149‑1 (MQ=255)
                                                                             tttATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCg                                                                    <  2:404188/149‑1 (MQ=255)
                                                                                  gACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACATGCTATTTTGCGTTTAg                                                               >  2:81068/1‑149 (MQ=255)
                                                                                   aCGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGc                                                              <  1:853278/149‑1 (MQ=255)
                                                                                              aGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCaa                                                   >  2:335975/1‑149 (MQ=255)
                                                                                               gATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCaaa                                                  <  1:217606/149‑1 (MQ=255)
                                                                                                aTTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAAc                                                 <  1:713206/149‑1 (MQ=255)
                                                                                                   ggTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCgg                                              <  1:224386/149‑1 (MQ=255)
                                                                                                     tGACAAAGTGCCTGTGTTTATGCAGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCt                                            <  2:470456/149‑1 (MQ=255)
                                                                                                       aCAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCtac                                          >  2:124029/1‑147 (MQ=255)
                                                                                                         aaaGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATGCGGCCTACAAAACCAATTAAATTCAATGAATTTAAAAATGATGTAGCCCAGATACATAGCTCAACAGGCTAATTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAg                                        >  1:918756/1‑149 (MQ=38)
                                                                                                           aGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACACCAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAt                                      >  1:670770/1‑149 (MQ=255)
                                                                                                                cTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGc                                 <  2:309200/149‑1 (MQ=255)
                                                                                                                    gtTTAAGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCATTGAATTAAAAAATGATGGAGACCTGATACATATCACAACAGGCGATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggt                             <  1:685694/149‑1 (MQ=38)
                                                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCGggtt                            >  1:65311/1‑149 (MQ=255)
                                                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtt                            <  2:974014/149‑1 (MQ=255)
                                                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggttt                           <  1:482097/149‑1 (MQ=255)
                                                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg                        >  2:995474/1‑149 (MQ=255)
                                                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACAAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAAAAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggt                       <  2:61442/149‑1 (MQ=255)
                                                                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtt                      >  2:725655/1‑149 (MQ=255)
                                                                                                                             tgATGCGGCGTTAACGCCTTATCCCGCCTTCAAAACCAATTAAATTCAATGAATTAAATAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTTCTTTAACCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                    <  2:171392/148‑1 (MQ=255)
                                                                                                                             ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAAATAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATATTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                    >  2:218841/1‑149 (MQ=255)
                                                                                                                              gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                   >  1:420856/6‑149 (MQ=255)
                                                                                                                              gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                   >  2:359547/6‑149 (MQ=255)
                                                                                                                                  cggcgTTACCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAAGGATGTAGCCCTGATCCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTc               <  1:707710/149‑1 (MQ=255)
                                                                                                                                   ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt              <  2:352260/149‑1 (MQ=255)
                                                                                                                                   ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAAGGATGTAGACCTGCTACATAGCACACCAGGCTATTTTGCGTTTAGCATCAGTCTCACACCGGGTCCAGATAGCGCCGGTTTTGGTTgtctgtct              >  2:771630/1‑149 (MQ=255)
                                                                                                                                        tAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGCATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTAGTTTGCGTTTAGCATCAGTCTCAAACCTGCTCCAGATAGAGCAGCTTTTGGTTGTCTGTATTAACg         >  2:676186/1‑149 (MQ=255)
                                                                                                                                         aaCGCCTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGcc       >  2:878075/1‑148 (MQ=255)
                                                                                                                                           cGCCTTATCCGGCCTACAAAACCAATTAAATTAAATGAATTAAAAAATGATGTAGAACTGATACATAGCACAACAGGCTATTATGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTTGTTTTCTGTATTAACGCAc      >  2:929679/1‑149 (MQ=255)
                                                                                                                                            gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGCCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAcc     >  2:782442/1‑149 (MQ=255)
                                                                                                                                               ttATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCAAAACAGGCTATTTTGCTTTTAGAATCAGTCTCAAAACGGCTCAAGATAGAGCAGGTTTTGGTTTTCTTTATTAACTCACCAta  >  2:651403/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATG  >  NC_000913/2917766‑2918057

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: