Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,979,525 T→C 6.2% intergenic (‑310/‑228) insA ← / → uspC IS1 repressor TnpA/universal stress protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,979,5250TC 6.2% 212.1 / 3.5 64intergenic (‑310/‑228)insA/uspCIS1 repressor TnpA/universal stress protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

ATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATA  >  NC_000913/1979276‑1979660
                                                                                                                                                                                                                                                         |                                                                                                                                       
atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAATAATGAAACAGGAATACTCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTA‑‑TTTAAACAGACCCTTACCGCACAACAACCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACTTATATTAACCa                                                                                                              >  2:6304‑M2/4‑275 (MQ=255)
                                                                   agcaaAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTACAGAAATGTTTTACATATATTAAACAATTGTTTAAGGCGGAGAATAAGTACAATATTTTATNATAATNNTAATTTAAATtaaaataatttt                                             >  3:687/1‑265 (MQ=255)
                                                                       aaaaTGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTg                                                             <  1:687/255‑1 (MQ=255)
                                                                                             aaTGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGttt                                                                                                                             >  2:544839/1‑169 (MQ=255)
                                                                                                 ggATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGtt                                                                                                                              >  2:274780/1‑164 (MQ=255)
                                                                                                 ggATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACa                                                                                                                         <  4:544834/169‑1 (MQ=255)
                                                                                                     gTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTAc                                                                                                                          <  4:274776/164‑1 (MQ=255)
                                                                                                       tGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACttaatgttaa                                                                                                                                       <  8:24404/149‑1 (MQ=255)
                                                                                                            acaGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGaaa                                                                                                                                  >  8:465557/1‑149 (MQ=255)
                                                                                                             caGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAAt                                                                                                                                 <  5:52034/149‑1 (MQ=255)
                                                                                                                 cAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGttt                                                                                                                             <  6:26841/149‑1 (MQ=255)
                                                                                                                  aaaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGtttt                                                                                                                            >  7:6291/1‑149 (MQ=255)
                                                                                                                  aaaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGtttt                                                                                                                            >  5:75577/1‑149 (MQ=255)
                                                                                                                   aaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTa                                                                                                                           <  5:220609/149‑1 (MQ=255)
                                                                                                                   aaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTa                                                                                                                           <  7:353578/149‑1 (MQ=255)
                                                                                                                    aaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTAc                                                                                                                          >  6:117686/1‑149 (MQ=255)
                                                                                                                      tAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACat                                                                                                                        <  8:372158/149‑1 (MQ=255)
                                                                                                                          tCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGtt                                                                                                                              >  2:541759/1‑139 (MQ=255)
                                                                                                                          ccACCGTAAATAATCATCTGCTATAATTAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATGAAATACAAAGAAAGGCAACAATCAACTTAATGCTNATGAAATGTATTACATATATTAACCACTTGTTTAATGCGGATAATCAGt                                                                                       >  3:45987/2‑178 (MQ=255)
                                                                                                                              cGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTAc                                                                                                                          <  4:541754/139‑1 (MQ=255)
                                                                                                                              cGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAAc                                                                                                                <  7:465557/149‑1 (MQ=255)
                                                                                                                              cGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACACACACTTACCGCACATCAAACTAAGTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACGTATATTAACCACTTGTTTAATGCGGAGAATCAGTacac                                                                                   <  1:45987/178‑1 (MQ=255)
                                                                                                                                   ataatCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACtt                                                                                                           >  6:371811/1‑149 (MQ=255)
                                                                                                                                       tcatcTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAGTACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGttt                                                                                                       >  5:138244/1‑149 (MQ=255)
                                                                                                                                          tcTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTa                                                                                                                           >  7:300429/1‑126 (MQ=255)
                                                                                                                                          tcTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACCAATTAAATACAAAGAAAGGCAACAATCAACTTCATGTTAATGAAATGTTTTa                                                                                                                           <  8:300429/126‑1 (MQ=255)
                                                                                                                                           cTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATg                                                                                                   <  7:54349/149‑1 (MQ=255)
                                                                                                                                                      aatCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAg                                                                                        <  7:260490/149‑1 (MQ=255)
                                                                                                                                                        tCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTa                                                                                      <  6:75577/149‑1 (MQ=255)
                                                                                                                                                            tttCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACt                                                                                  >  5:85240/1‑149 (MQ=255)
                                                                                                                                                             ttCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTGATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGt                                       >  2:319270/1‑191 (MQ=255)
                                                                                                                                                                 tGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTGATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGtgatt                                   <  4:319266/191‑1 (MQ=255)
                                                                                                                                                                  gCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            <  6:151478/149‑1 (MQ=255)
                                                                                                                                                                   cAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAATAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCt                                                                           <  6:85240/149‑1 (MQ=255)
                                                                                                                                                                        ccAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATaa                                                                      <  5:309919/149‑1 (MQ=255)
                                                                                                                                                                                          aaCAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAAt                                                         >  8:245085/1‑144 (MQ=255)
                                                                                                                                                                                          aaCAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAAt                                                         <  7:245085/144‑1 (MQ=255)
                                                                                                                                                                                          aaCAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTatat                                                    <  5:117686/149‑1 (MQ=255)
                                                                                                                                                                                           aCAGACACTTAGCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAAAGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGCATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATg                                                   >  5:388060/1‑149 (MQ=255)
                                                                                                                                                                                             aGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATtat                                                      >  7:84255/1‑144 (MQ=255)
                                                                                                                                                                                             aGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATtat                                                      <  8:84255/144‑1 (MQ=255)
                                                                                                                                                                                                 acTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCAcc                                             >  8:283451/1‑149 (MQ=255)
                                                                                                                                                                                                 acTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCAcc                                             >  6:105299/1‑149 (MQ=255)
                                                                                                                                                                                                               aaaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            >  7:120624/1‑104 (MQ=255)
                                                                                                                                                                                                               aaaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            <  8:120624/104‑1 (MQ=255)
                                                                                                                                                                                                               aaaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTgattgatt                               >  5:321211/1‑149 (MQ=255)
                                                                                                                                                                                                                aaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGatta                              >  8:419958/1‑149 (MQ=255)
                                                                                                                                                                                                                 aCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGattaa                             <  7:283451/149‑1 (MQ=255)
                                                                                                                                                                                                                 aCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGattaa                             <  5:105299/149‑1 (MQ=255)
                                                                                                                                                                                                                     aTTAAATACAAAGAAAGGCAACAATCAACTTAATGTTGATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTAt                         >  5:286374/1‑149 (MQ=255)
                                                                                                                                                                                                                      ttAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            >  5:188954/1‑97 (MQ=255)
                                                                                                                                                                                                                      ttAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            <  6:188954/97‑1 (MQ=255)
                                                                                                                                                                                                                           tACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGga                                                                                               >  1:101771/1‑73 (MQ=255)
                                                                                                                                                                                                                           tACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACt                                                                                  <  6:458359/86‑1 (MQ=255)
                                                                                                                                                                                                                           tACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACt                                                                                  >  5:458359/1‑86 (MQ=255)
                                                                                                                                                                                                                           tACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTAc                   >  8:8130/1‑149 (MQ=255)
                                                                                                                                                                                                                            aCAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACt                  <  7:419958/149‑1 (MQ=255)
                                                                                                                                                                                                                               aagaaagGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGTAt                                                                                           <  3:101771/73‑1 (MQ=255)
                                                                                                                                                                                                                                 gaaagGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACtaataa             >  8:424825/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aagGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATc           <  6:354087/149‑1 (MQ=255)
                                                                                                                                                                                                                                    agGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCa          <  6:388060/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGa        <  8:6291/149‑1 (MQ=255)
                                                                                                                                                                                                                                       ctacaaTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGAc       <  6:138244/147‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAat   <  6:286374/149‑1 (MQ=255)
                                                                                                                                                                                                                                            aTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAata  >  8:240168/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                       
ATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATA  >  NC_000913/1979276‑1979660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑