Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,764,710 T→C intergenic (‑324/+224) sufA ← / ← ydiH Fe‑S cluster assembly protein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7100TC90.8% 196.8 / ‑3.1 65intergenic (‑324/+224)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base C (37/22);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 7.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAA  >  NC_000913/1764567‑1764855
                                                                                                                                               |                                                                                                                                                  
aTAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTAcc                                                                                                                                            >  1:701697/1‑151 (MQ=255)
 tAATGGTTATTATCTAAGGAGCATTAAGGCTGCGATGAAAAAGAAAAAGTTGTATTGAAAATAACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTAccc                                                                                                                                           >  1:58205/1‑151 (MQ=255)
  aaTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTAcccc                                                                                                                                          >  1:383000/1‑151 (MQ=255)
    tGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCAt                                                                                                                                        >  2:251840/1‑151 (MQ=255)
        tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCgg                                                                                                                                    >  1:297144/1‑151 (MQ=255)
        tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCg                                                                                                                                     >  2:953571/1‑150 (MQ=255)
         attatCTAAGGAGCATTAAGTCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGAGCAGCGTTCAGAAATAAGAAAACTCTTAAGTCTTTGCGACACAGGCTTAAGGGTTTCTACCCCAt                                                                                                                                        >  2:893964/1‑146 (MQ=255)
         attatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGgc                                                                                                                                   >  2:769786/1‑151 (MQ=255)
                       aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGc                                                                                                                      >  2:337098/1‑150 (MQ=255)
                        ttAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                                                                                    >  1:66409/1‑151 (MQ=255)
                         tAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCAct                                                                                                                   >  1:453172/1‑151 (MQ=255)
                                gtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGt                                                                                                             <  1:829256/150‑1 (MQ=255)
                                gtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                                                                                            >  1:457436/1‑151 (MQ=255)
                                 tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTgg                                                                                                           >  2:1026224/1‑151 (MQ=255)
                                     gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCtt                                                                                                        >  1:167073/1‑150 (MQ=255)
                                      aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTa                                                                                                       <  1:701362/150‑1 (MQ=255)
                                      aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAg                                                                                                      <  2:345266/151‑1 (MQ=255)
                                       aaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                                                                                                     <  1:999240/151‑1 (MQ=255)
                                        aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGct                                                                                                    >  2:596118/1‑151 (MQ=255)
                                          agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctct                                                                                                  <  2:66536/151‑1 (MQ=255)
                                          agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctct                                                                                                  >  2:247333/1‑151 (MQ=255)
                                           gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCtt                                                                                                 >  1:668291/1‑151 (MQ=255)
                                                    gTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCg                                                                                        <  1:797292/151‑1 (MQ=255)
                                                      aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCCTCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGAGGCGAt                                                                                      >  2:719632/1‑151 (MQ=255)
                                                        tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATa                                                                                     <  1:251865/150‑1 (MQ=255)
                                                         gAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATa                                                                                     >  2:434157/1‑149 (MQ=255)
                                                           aaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                                                                                 <  2:456599/151‑1 (MQ=255)
                                                                aCTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc                                                                            <  2:541226/151‑1 (MQ=255)
                                                                   aTTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                                          >  2:50960/1‑150 (MQ=255)
                                                                   aTTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCa                                                                         <  2:807714/151‑1 (MQ=255)
                                                                    tttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt                                                                        >  2:49991/1‑151 (MQ=255)
                                                                    tttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt                                                                        >  1:493169/1‑151 (MQ=255)
                                                                     ttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt                                                                        <  1:730052/150‑1 (MQ=255)
                                                                       aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                      <  1:998780/150‑1 (MQ=255)
                                                                         gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                                   >  1:921831/1‑151 (MQ=255)
                                                                         gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                                   <  2:1066804/151‑1 (MQ=255)
                                                                             tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                                               >  2:284067/1‑151 (MQ=255)
                                                                                gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccg                                                            <  1:759782/151‑1 (MQ=255)
                                                                                gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccg                                                            <  1:405689/151‑1 (MQ=255)
                                                                                    aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGt                                                        >  2:880737/1‑151 (MQ=255)
                                                                                    aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGt                                                        >  2:565745/1‑151 (MQ=255)
                                                                                     aaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtt                                                       >  1:293791/1‑151 (MQ=255)
                                                                                                  cAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTg                                          >  1:846599/1‑151 (MQ=255)
                                                                                                   aaaaaaaaaaaaCCCCTTAAGTCTGTGCGACACAGTTTAAAGGGCTTCTCCCCCATCCGGCTTTTATCTCGGTCTCTCTCAGGGTCTTATCTTTTGAAGGGGAGATAAAAATAATCTCAAAAATCAAACCCGCCGTTTTAACACaa                                              <  2:365119/144‑1 (MQ=255)
                                                                                                    gAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCg                                        >  1:383941/1‑151 (MQ=255)
                                                                                                    gAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCAATCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCg                                        >  2:788505/1‑151 (MQ=255)
                                                                                                     aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCg                                        >  2:730447/1‑150 (MQ=255)
                                                                                                       aTAAAAAACCCCTTAAGTTGGTGCGACACAGGCTAAGGGTCTTCTCCCCCACCGGGCTCTTATCCCCGGCACTCTCGGTGGTTTAGCTTTTGAAGGGGCAATAAGAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                     <  2:959120/151‑1 (MQ=255)
                                                                                                         aaGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                     >  2:276305/1‑149 (MQ=255)
                                                                                                            aaaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtat                                >  2:707862/1‑151 (MQ=255)
                                                                                                             aaccccTTAAGTCTGTGCGACACAGGTTTAAGGGCTTCTCCCCCATCCGGCGCTTATCTCCGGCCCTCTCAGTGGTTTAGCTCTTGAAGGGGCGAAAAAAATAATCTAATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGtatt                               <  1:831860/148‑1 (MQ=255)
                                                                                                                  cTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                          >  1:819608/1‑151 (MQ=255)
                                                                                                                   ttagttcggtgCGAACCAGGCTTAAGGCTTTCTCCCCCACCCGGCGCTTATCTCCGGCCCTCTCAGTGTTTTAGCTTTTGAAGGGGCAAAAAAAATAATCTCAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                          <  2:862097/142‑1 (MQ=255)
                                                                                                                    tAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                        >  2:905852/1‑151 (MQ=255)
                                                                                                                     aattCTGTGCGACGCAGGCTAAGGGTCTTCTCCCCCATCCGGCGCTTACCTCCGCCACTCTCAGTGTCTTAGCTCTTGAGGGGGCGATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt                         <  2:484672/146‑1 (MQ=255)
                                                                                                                       gtcggggCGACACGGGCTTAAGGGCTTCTCCCCCATCGGGCTTTTATCCCCGGCTCTCTCAGGGCTTTAGCTTTTGAGGGGGCAAAAAAAATAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGtatta                              <  1:829678/136‑1 (MQ=255)
                                                                                                                        tCTGTGCGACACAGGCTTAAGGGTTTCTCCCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTTTTGAAGGGGCAATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                      <  2:325300/149‑1 (MQ=255)
                                                                                                                        tCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAAc                     >  1:967415/1‑150 (MQ=255)
                                                                                                                               gACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc               >  2:1050719/1‑149 (MQ=255)
                                                                                                                               gACACAGGCTAAAGGGCTTCTACCCCACCCGGCGCTTACCCCCGGCACTCTCGGTGGCTTAGCTCTTGAGGGGGCGATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAata                 <  1:127561/147‑1 (MQ=255)
                                                                                                                                  actggCTTAAGGATTTCCCCCCCAGCCGGCGTTTACCGCGGTGCCCCGCAGTGGCTTATCTTTTGAAGGGGTTATAAAAATAATCTAAAAATCCAATCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc           <  2:1053774/147‑1 (MQ=255)
                                                                                                                                  acaGGCTTAAGGTCTTCCCCCCCACCCGGCGTTTATCCCCGCCACTCTCAGTGTTTTAGCTCTAAAAGGGGCAATAAGAATACTCTAATAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc           <  1:565135/150‑1 (MQ=255)
                                                                                                                                  acaGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCg          >  2:787198/1‑151 (MQ=255)
                                                                                                                                    aGGCTAAAGGGCTCCGCCCCCACCCGCCGTTTATCCCCGACCCCCTCAGTGGCTTAGTTTTTAAGGGGGCAATAAGAATAATCCCAAAAATCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAgg            <  1:1036170/147‑1 (MQ=255)
                                                                                                                                          aaGGGCTTCTCCCCCTCCCGGCGCTTATCTCCGGCACTCACAGTGGTTTAGCTTTTGAAGGGGCGATAAAAAAAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATaa  <  2:637856/151‑1 (MQ=255)
                                                                                                                                            ggctttCTACCCCACCGGGCGCTTATCCCCGCCCCTCACATGGGCTTAGCTTTTGAAGGGGCAAAAAAAATATTCACAAAAAGATAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGAt    <  2:405307/144‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAA  >  NC_000913/1764567‑1764855

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: