Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,694 G→C 5.7% intergenic (‑308/+19) 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,6940GC 5.7% 235.8 / 4.4 71intergenic (‑308/+19)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (40/27);  new base C (2/2);  total (42/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 = 7.53e-01

CATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAG  >  NC_000913/1764424‑1764964
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                               
cATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGaca                                                                                                                                                                                                                                                                              >  2:508085/1‑274 (MQ=255)
                                                       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)
                                                                                                                  aCGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTACGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCg                                                                                                                                                                                                                                                                                 >  4:522123/1‑157 (MQ=255)
                                                                                                                    nttAAGAATTAGGTATCTTATTGTTAGATAATGGTTATTATCTAAGGAGCATTANGGCTNTGATGACACAGAAAAAGTTGTATNGAAAATGACTATTTAAGAGATAGGTACAAAAAGGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGGGCGACACAGGCTTAAGGGTTTCTACCCCCTCCGGCGCTTATCTCCGGCACTCTCa                                                                                                                                                                                                                            <  4:508080/205‑1 (MQ=255)
                                                                                                                      tAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGacac                                                                                                                                                                                                                                                                             <  2:522128/157‑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)
                                                                                                                                                      ttattaTCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAg                                                                                                                                                                                                                                                                           <  6:200943/127‑1 (MQ=255)
                                                                                                                                                      ttattaTCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAg                                                                                                                                                                                                                                                                           >  5:200943/1‑127 (MQ=255)
                                                                                                                                                       tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                                                                                                                                                         >  7:65289/1‑128 (MQ=255)
                                                                                                                                                       tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                                                                                                                                                         <  8:65289/128‑1 (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)
                                                                                                                                                                                    gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGaca                                                                                                                                                                                                                                                                              <  4:258278/94‑1 (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                                                                                                                                                                                          <  8:105913/133‑1 (MQ=255)
                                                                                                                                                                                                                                tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                                                                                                          >  7:105913/1‑133 (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)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                               
CATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAG  >  NC_000913/1764424‑1764964

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