Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,019 G→T 5.9% intergenic (+206/+26) 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,0190GT 5.9% 237.9 / 4.1 68intergenic (+206/+26)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base G (38/26);  new base T (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 = 8.01e-01

ACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAATCCAGTGCGT  >  NC_000913/2917759‑2918287
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                            
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                                                                                                                                                                                                                                                >  7:22545/1‑110 (MQ=255)
                                                                                                                                                                                     tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAAc                                                                                                                                                                                                                                                <  8:22545/110‑1 (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                                                                                                                                                                                                                         <  5:210602/97‑1 (MQ=255)
                                                                                                                                                                                                                         aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                                                                                                                                         >  6:210602/1‑97 (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                                                                                                                                                                               <  7:145998/108‑1 (MQ=39)
                                                                                                                                                                                                                                                      caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgca                                                                                                                                                                               >  8:145998/3‑110 (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)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                            
ACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAATCCAGTGCGT  >  NC_000913/2917759‑2918287

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