Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,022 A→G 5.9% intergenic (+209/+23) 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,918,0220AG 5.9% 220.1 / 3.7 68intergenic (+209/+23)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base A (38/26);  new base G (2/2);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

ACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAATCCAGTGCGTATCAATAGC  >  NC_000913/2917759‑2918296
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
aCTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                                                                                                                                                          >  2:224284/1‑275 (MQ=255)
                                                                ttgtcggtggaaaTCGGCTTTTTTTGACGCTGGTTTAAGATTGGTTACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGa                                                                                                                                                                                                          <  4:293274/263‑1 (MQ=255)
                                                                                                                          tgtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAGATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGTCGGTTTTGGTTTTCTGTCTTAACGCCCCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCCTTCCCATCCCATCGTCCCGCCCGCCCATCCCGTGTTTCTTATACAGTTTATg                                                                                                                                               >  2:93117/1‑275 (MQ=255)
                                                                                                                             ttATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCAGAACAGGCTATTTTGTGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttt                                                                                                                                                                                                                                                                           <  6:388427/149‑1 (MQ=255)
                                                                                                                               aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttg                                                                                                                                                                                                                                                                         >  7:194670/1‑149 (MQ=255)
                                                                                                                                     gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                                                                                                                                                   >  5:90477/6‑149 (MQ=255)
                                                                                                                                     gATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                                                                                                                                                   <  5:54402/149‑1 (MQ=255)
                                                                                                                                      ataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                                                  <  6:90477/145‑1 (MQ=255)
                                                                                                                                          ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTACTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt                                                                                                                                                                                                                                                              >  5:340947/1‑149 (MQ=255)
                                                                                                                                               tAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                                            >  2:236670/1‑241 (MQ=255)
                                                                                                                                                   gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATa                                                                                                                                                        <  4:236666/241‑1 (MQ=255)
                                                                                                                                                     cTTATCCGGCCTAGAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCAt                                                                                                                                                                                                                                                   <  8:301776/149‑1 (MQ=255)
                                                                                                                                                             gCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGACTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgc                                                                                                                                                                                                                                           <  8:380885/149‑1 (MQ=255)
                                                                                                                                                              ccTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCCCATGGGCTTTCATCACTTTATCCATATCGATTTCTACAcc                                                                                                           >  2:492612/1‑275 (MQ=255)
                                                                                                                                                                          aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTGCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGACCAGCCAGGAATCAGATAc                                                                                                                                                                                                       >  2:515898/1‑171 (MQ=255)
                                                                                                                                                                          aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAt                                                                                                                                                                                                              >  2:343932/1‑164 (MQ=255)
                                                                                                                                                                              aaTTCAATGAATANAAAANTGATGTAGACCTGATACATAGCACAACATGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCGCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGAACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCa                                                                                                                                                                                                   <  4:515893/171‑1 (MQ=255)
                                                                                                                                                                               aTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGaa                                                                                                                                                                                                                             >  6:199510/1‑144 (MQ=255)
                                                                                                                                                                               aTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGaa                                                                                                                                                                                                                             <  5:199510/144‑1 (MQ=255)
                                                                                                                                                                                     tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAAc                                                                                                                                                                                                                                                          <  8:22545/110‑1 (MQ=255)
                                                                                                                                                                                     tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAAc                                                                                                                                                                                                                                                          >  7:22545/1‑110 (MQ=255)
                                                                                                                                                                                     tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                                                                                                                                                  >  6:106466/1‑149 (MQ=255)
                                                                                                                                                                                            aaaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAg                                                                                                                                                                                                           >  7:429570/1‑149 (MQ=255)
                                                                                                                                                                                                   atgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCCTCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGGTCCATATCGATTTCTACACCCAGCCCCGTTTTTTCTGGCCACTTTACCAGCACCCCt                                                                      >  4:347279/1‑275 (MQ=255)
                                                                                                                                                                                                         aCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATTCTGCATTccc                                                                                                                                                                                              >  7:10590/1‑149 (MQ=255)
                                                                                                                                                                                                         aCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                                      >  2:120120/1‑173 (MQ=255)
                                                                                                                                                                                                            tGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTACAGCCAGGAATCAGATACTGCATTCCcatc                                                                                                                                                                                           >  7:23030/1‑149 (MQ=255)
                                                                                                                                                                                                             gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaa                                                                                                                                                                                                               >  1:71443/1‑128 (MQ=255)
                                                                                                                                                                                                             gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgt                                                                                                                                                                  <  4:120119/173‑1 (MQ=255)
                                                                                                                                                                                                              aTACATATCACAACAGACTATTTTGCGTTTAGCATCTGTCTCCAACCGGCTCCAGATAGAGCCGGTTTTGGTTTGCTGTCTTAACTCACCATGCACGGGCGCTTGTTATCGCAAGTCCAGCCAGCAATCCAATACTGCATTCCCATCCAATCGTCACTCCCGCCCCGCCCGGTTTTCTGATAGACCTCATGTGCTTCAATACCTTAATCCTTATCGTTTTCTACACCTACCCCGGGTTTTTACGGCACTGGTACCATCCCCCCCTTgccaccaa                                                            >  4:42479/1‑267 (MQ=255)
                                                                                                                                                                                                                aCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                                                                                                                                                    >  6:424938/1‑73 (MQ=255)
                                                                                                                                                                                                                aCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                                                                                                                                                    <  5:424938/73‑1 (MQ=255)
                                                                                                                                                                                                                 cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAg                                                                                                                                                                                                           <  3:71443/128‑1 (MQ=255)
                                                                                                                                                                                                                    aGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                                      >  2:137499/1‑162 (MQ=255)
                                                                                                                                                                                                                     gCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCGGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATTCCAAACGGtttt                                                    >  2:400630/1‑273 (MQ=255)
                                                                                                                                                                                                                     gCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTAGTTTTCTGGCTTAACGCACCATGCACGGGCGCTTGTTATCGAAGGTCCAGCCCGGAATCAGATACTGCATTCCCATGGCAACGTCACGACCGCCAAGCCCGTGTTTCTGCTACCGCTAATGGCCTTTCATTACTTTCTCCATTTCGATTTCCCCAACCAACCCCGGTTTTCGTGGCCACTGCAACAGCCCCCCTTTGACCTCAACCGGTTCt                                                    >  4:399589/1‑275 (MQ=255)
                                                                                                                                                                                                                        caacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgt                                                                                                                                                                  <  4:137498/162‑1 (MQ=255)
                                                                                                                                                                                                                        caacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             >  4:413247/1‑167 (MQ=255)
                                                                                                                                                                                                                         aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                                                                                                                                                  >  6:210602/1‑97 (MQ=255)
                                                                                                                                                                                                                         aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                                                                                                                                                  <  5:210602/97‑1 (MQ=255)
                                                                                                                                                                                                                         aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACATCTCATGGGCTTTCATCACTTCATCCATATCGATTGCTACACCCAGCCCCGGTTTTTCTGGCATCTGTCCCCGCCCCACTTTTATCCCAAACGGTTATTTgg                                                >  1:75256/1‑275 (MQ=255)
                                                                                                                                                                                                                          acaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcgc                                                                                                                                                                             <  8:194670/149‑1 (MQ=255)
                                                                                                                                                                                                                            aGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGAt                                                                                                                                                         <  2:413252/167‑1 (MQ=255)
                                                                                                                                                                                                                              gCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg                                                                                                                                                                                                                                                                        >  2:468818/1‑55 (MQ=255)
                                                                                                                                                                                                                                  ttttGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                                                                                                                                                    <  4:468813/55‑1 (MQ=255)
                                                                                                                                                                                                                                   tttGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGccc                                                                                                                                                                    >  7:255775/1‑149 (MQ=255)
                                                                                                                                                                                                                                      gCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgtg                                                                                                                                                                 <  8:429570/149‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCCTTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATa                                                                                                                                                        >  7:420823/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aTCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATa                                                                                                                                                        >  6:89257/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  aGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                                            >  2:361820/1‑78 (MQ=255)
                                                                                                                                                                                                                                                   gTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgtg                                                                                                                                                                 >  1:116324/1‑136 (MQ=255)
                                                                                                                                                                                                                                                      tcAAACCGGCTCCAGAGAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcc                                                                                                                                                                                                                        <  4:361815/78‑1 (MQ=255)
                                                                                                                                                                                                                                                      caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgca                                                                                                                                                                                        >  8:145998/3‑110 (MQ=39)
                                                                                                                                                                                                                                                      caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgca                                                                                                                                                                                        <  7:145998/108‑1 (MQ=39)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATTCTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCAt                                                                                                                                                <  8:10590/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             <  3:116324/136‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                                                                                                                                                          >  5:94855/1‑139 (MQ=255)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                                                                                                                                                          <  6:94855/139‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             >  8:313483/1‑135 (MQ=255)
                                                                                                                                                                                                                                                        aaaCCGGCTCTATCTGGAGCCGGTTATGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             <  7:313483/135‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCCCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTTCTGCCAAATCCa                 >  2:466826/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgtg                                                                                                                                                                 >  2:298031/1‑125 (MQ=255)
                                                                                                                                                                                                                                                              gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTTCCTTCCTGCCAAATCCAGTGcct           >  2:269083/1‑273 (MQ=255)
                                                                                                                                                                                                                                                               cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtc                                                                                                                                        <  5:106466/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                tCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                                       >  5:220426/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  cAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             <  4:298027/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACGGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACg                                                         <  2:157974/225‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtcatca                                                                                                                                    >  6:249322/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       agagTCGGTTTATGTTTTCTGTCTTAACTACCCATACACGGGCGGCTTGTCTCGAACGGCGAGCCAGGAACCAGAGTCTGAATTCCCATCGCATCGTCACGAGCGCCAAGCCCGTGTTTCTTATACAACTCATGGGCTTTCATCACTTGTTCCATATCGCTTGCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACGAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAAACCAGTGCGTATCAATAGc  <  4:400625/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
ACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAATCCAGTGCGTATCAATAGC  >  NC_000913/2917759‑2918296

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: