Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,701 G→A 5.7% intergenic (‑315/+12) sufA ← / ← rydB Fe‑S cluster assembly protein/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7010GA 5.7% 219.2 / 3.3 70intergenic (‑315/+12)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (39/27);  new base A (2/2);  total (41/29)
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.70e-01

AAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAG  >  NC_000913/1764479‑1764964
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                        
aaaaCATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCg                                                                                                                                                                                                                                     >  4:268524/1‑259 (MQ=255)
    cATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                                                                                                                                                                                                 <  2:268528/259‑1 (MQ=255)
                          taaGCATAGACCCTGATGGGTGGGTTAACAGGCACGCTAAGAATTAGGTAACTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa                                                                                                                                                                                           <  2:169067/275‑1 (MQ=255)
                                        gATGTGTGGGTTAACAGGGACGCTAAGAATTAGGTATCTCATTGGTAGATGATGGGTATTATCTAAGAAGCATTCAGGCTGTGTTGAAAAAGACAAAGTTTTATTGAAAATGACTATTTAAGACATAGGTAAAAAAGTCCAGCGTTCAGAAATACGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAAGGTTTCTACTCAATCCGGCTCTTATCCACGGCACCCTCAGTGGCTTGCATCTTGAAGGGGCGGAACTAATAATCCCAGAAAGCaaa                                                                                                                                                                             >  4:190518/1‑275 (MQ=255)
                                                             nttAAGAATTAGGTATCTTATTGTTAGATAATGGTTATTATCTAAGGAGCATTANGGCTNTGATGACACAGAAAAAGTTGTATNGAAAATGACTATTTAAGAGATAGGTACAAAAAGGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGGGCGACACAGGCTTAAGGGTTTCTACCCCCTCCGGCGCTTATCTCCGGCACTCTCa                                                                                                                                                                                                                            <  4:508080/205‑1 (MQ=255)
                                                                             ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                                                                                                        >  2:202752/1‑275 (MQ=255)
                                                                               cATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTaa                                                                                                                                                                                                                                                                     <  8:330053/149‑1 (MQ=255)
                                                                                  tGTTAGATAATGGTAATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTAATGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTCCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCAACCGGAGCGTATCTCCGGCACTCTCAGTGCCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAGTAGTATTATTTTTGaa                                                                                                                                   <  2:190520/275‑1 (MQ=255)
                                                                                    ttAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                                                                                                                                                                   <  7:357297/145‑1 (MQ=255)
                                                                                    ttAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                                                                                                                                                                   >  8:357297/1‑145 (MQ=255)
                                                                                          aaTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAcc                                                                                                                                                                                                                                                         >  5:277828/1‑149 (MQ=255)
                                                                                           aTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAccc                                                                                                                                                                                                                                                        >  8:259989/1‑149 (MQ=255)
                                                                                                tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                                                                                                                                                         <  8:65289/128‑1 (MQ=255)
                                                                                                tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                                                                                                                                                         >  7:65289/1‑128 (MQ=255)
                                                                                                  ttatCTAAGGAGCATTAAGGGTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGAGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtttt                                                                                                                                                                  <  1:85789/228‑1 (MQ=255)
                                                                                                         aGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGc                                                                                                             >  2:421084/1‑275 (MQ=255)
                                                                                                         aGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAt                                                                                                                                                                                                                                          >  8:237691/1‑149 (MQ=255)
                                                                                                               aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                                                                                                                                                                                                    >  2:533131/1‑149 (MQ=255)
                                                                                                               aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                                                                                                                                              >  4:263573/1‑219 (MQ=255)
                                                                                                                   aGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAct                                                                                                                                                                                                                                <  4:533126/149‑1 (MQ=255)
                                                                                                                   aGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaaa                                                                                                                                                          <  2:263577/219‑1 (MQ=255)
                                                                                                                        gtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCa                                                                                             >  2:320823/1‑275 (MQ=255)
                                                                                                                         tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                                                       >  2:150496/1‑232 (MQ=255)
                                                                                                                            tgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGtt                                                                                                                                                                                                                                                               <  5:218507/109‑1 (MQ=255)
                                                                                                                            tgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGtt                                                                                                                                                                                                                                                               >  6:218507/1‑109 (MQ=255)
                                                                                                                             gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                                                   <  4:150495/232‑1 (MQ=255)
                                                                                                                              aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTaa                                                                                                                                                                                                                                                                     >  6:119842/1‑102 (MQ=38)
                                                                                                                                 aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                                                                                                                                                                                                                  <  6:277828/149‑1 (MQ=255)
                                                                                                                                   gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                                                                                                                                                                                                                >  6:156941/1‑149 (MQ=255)
                                                                                                                                          ttGTATTGAAAATGACTATTTCAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                                                                                                                                                                                                 <  6:350638/125‑1 (MQ=255)
                                                                                                                                          ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                                                                                                                                                                                                 >  5:350638/1‑125 (MQ=255)
                                                                                                                                                   aaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAg                                                                                                                                                                                                >  8:241919/1‑149 (MQ=255)
                                                                                                                                                         cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                                                                                                          >  7:227148/1‑149 (MQ=255)
                                                                                                                                                               aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                                                    >  7:34521/1‑149 (MQ=255)
                                                                                                                                                                  agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt                                                                                                       >  2:327424/1‑223 (MQ=255)
                                                                                                                                                                    aTAGGTAAAAAAGTGCAGCGTTCAGACATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAg                                                                                                                                                                                                                           <  8:342650/105‑1 (MQ=255)
                                                                                                                                                                    aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAg                                                                                                                                                                                                                           >  7:342650/1‑105 (MQ=255)
                                                                                                                                                                    aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                                               >  8:359505/1‑149 (MQ=255)
                                                                                                                                                                      aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctct                                                                                                                                                                                                               >  2:364629/1‑115 (MQ=255)
                                                                                                                                                                       ggTAAAAAAGTGCAGCGTTCAGCCATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGTGTTTCTACCCCATCCGTCGCTTACCTCCGGCACTCTCACTGGCTTAGCTCTTGAAGTGGCGATAAGAATAATCTCATAATGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGTCGTTAGATCAGCAGCATTAAGCAGGTAATCGAACCTGGACGGAGTTGAATACAGGTCGATTTGGCGAACTTGCg                                               >  4:90911/1‑274 (MQ=255)
                                                                                                                                                                         tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                                                                                                          >  7:105913/1‑133 (MQ=255)
                                                                                                                                                                         tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                                                                                                          <  8:105913/133‑1 (MQ=255)
                                                                                                                                                                         tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc                                                                                                                                                                          >  5:449585/1‑149 (MQ=255)
                                                                                                                                                                          aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGaa                                                                                                                                                                                                           <  4:364624/115‑1 (MQ=255)
                                                                                                                                                                                     gCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                                                                                                                                              <  5:156941/149‑1 (MQ=255)
                                                                                                                                                                                     gCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                                                                                                                                              <  7:241919/149‑1 (MQ=255)
                                                                                                                                                                                      cGTTCATAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCTGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                                                                                                                                                      <  2:382731/92‑1 (MQ=255)
                                                                                                                                                                                        ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaa                                                                                                                                                           >  5:140827/1‑149 (MQ=255)
                                                                                                                                                                                        ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt                                                                                                       >  4:76067/1‑201 (MQ=255)
                                                                                                                                                                                            gAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga                                                                                                   <  2:76068/201‑1 (MQ=255)
                                                                                                                                                                                            acaaTAAGATAACCCTGAAGTCTGTGCGACACAGGTTTAAGTGGTTCTCCCCCCTCCGTCGTTTTGCTCCGGCACTCTCAGTTGTTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCa                         <  4:202749/273‑1 (MQ=255)
                                                                                                                                                                                             aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                                                      >  6:277682/1‑149 (MQ=255)
                                                                                                                                                                                             aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                                                      <  7:237691/149‑1 (MQ=255)
                                                                                                                                                                                                   gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                                                                >  5:381079/1‑149 (MQ=255)
                                                                                                                                                                                                           ttAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                                                                       <  7:192434/118‑1 (MQ=38)
                                                                                                                                                                                                           ttAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                                                                       >  8:192434/1‑118 (MQ=39)
                                                                                                                                                                                                            taggTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                                                       <  8:34521/146‑1 (MQ=255)
                                                                                                                                                                                                             aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTGGCAGGGGCGATAAGATTAATCTCATAAAGCTAACACGCCGTTTTAACACAAACGGCGATTAGTATTATTTTTGAACAATAACAGGCGGTAGATAAGCAGTATTAAGAAGGTAATTGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCCAACTTGCCGCATTAATTCAGGAGCAATGCACGCCCAAACGACaaggc        >  4:428628/1‑271 (MQ=255)
                                                                                                                                                                                                                tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                                                   <  6:381079/149‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                                                   >  5:207384/1‑149 (MQ=255)
                                                                                                                                                                                                                aggcacagACCCAGGTGTAAGGGTTCNNNCCCCGTCCGGGCTTTTTCTCCGGCACTCTCTTTGTCTTGACTCTTGCAGGGGCTATCAGAATAATCTCATAAATCTCACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAAAATCCGGCGGTAGATAAGCAGTATTAAGAAAGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCa        <  3:154004/265‑1 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                                                                                                                                                            >  1:33930/1‑105 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                                                                                                                                                >  2:308621/1‑117 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACACTATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAATTCAGGATCAATGCCCGCCAAGACTCCAGCCAAAGTAg  >  2:4249/1‑275 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                                                     >  1:38582/1‑192 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCgg                                              >  4:329982/1‑231 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTANTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTANATTTGGCGAATTTGCGGCATTAAGTTAGGANNAATGCACGCACAGCCCCCAGGCAAAATAg  >  3:57771/1‑275 (MQ=255)
                                                                                                                                                                                                                       gACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgc                                                                                                                                                                        <  3:33930/105‑1 (MQ=255)
                                                                                                                                                                                                                       gACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCAtt                                          <  2:329986/231‑1 (MQ=255)
                                                                                                                                                                                                                           caGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                                                        >  8:101641/1‑149 (MQ=255)
                                                                                                                                                                                                                              gCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                                                                     <  5:277682/149‑1 (MQ=255)
                                                                                                                                                                                                                              gCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCCCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGCACACNATCAGGCGGTAGa                                                                                                                  >  3:67127/1‑152 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                        
AAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAG  >  NC_000913/1764479‑1764964

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: