Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,078 G→A 6.5% intergenic (+87/+320) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,0780GA 6.5% 975.7 / 37.8 315intergenic (+87/+320)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (125/166);  new base A (4/16);  total (131/184)
Fisher's exact test for biased strand distribution p-value = 5.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

CCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTA  >  NC_000913/1551929‑1552227
                                                                                                                                                     |                                                                                                                                                       
ccGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATg                                                                                                                                                         >  1:1086484/1‑150 (MQ=255)
  gcgcTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGGGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcg                                                                                                                                                       >  2:4291258/1‑150 (MQ=255)
   cgcTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                                                                                      <  2:6740855/150‑1 (MQ=255)
   cgcTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                                                                                      >  1:2946575/1‑150 (MQ=255)
     cTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                                                                                    <  1:2350244/150‑1 (MQ=255)
         tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                <  2:6364779/150‑1 (MQ=255)
           tCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                                                                                              <  1:6589794/150‑1 (MQ=255)
            cGTGAAATCGAGAAGGCAGACGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                                             <  2:125777/150‑1 (MQ=255)
             gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGACGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                            <  1:5720480/150‑1 (MQ=255)
             gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                            >  2:442441/1‑150 (MQ=255)
             gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                            <  2:5041800/150‑1 (MQ=255)
              tGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                           <  1:5048307/150‑1 (MQ=255)
               gAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                          <  2:3371251/150‑1 (MQ=255)
                aaaTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAACGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                         <  1:6068827/150‑1 (MQ=255)
                    cGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                     <  2:6104909/150‑1 (MQ=255)
                      agaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                   <  1:2038800/150‑1 (MQ=255)
                      agaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                   <  1:1386660/150‑1 (MQ=255)
                       gtaGGCAGTAGCGAAAAATGAGAGTGAAGACGGGAGATAATCTCTTTTTAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGTCAACGTAGGCTTTGTTCATGCCGGACGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                  <  1:4400490/148‑1 (MQ=21)
                       gaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                  <  2:819215/150‑1 (MQ=255)
                       gaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                  <  2:668168/150‑1 (MQ=255)
                           gCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGa‑‑a                                                                                                                            >  1:660359/1‑150 (MQ=255)
                           gCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa‑‑a                                                                                                                            >  1:3857016/1‑150 (MQ=255)
                           gCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                    >  1:1208129/1‑144 (MQ=25)
                               aaGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAAGTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAAcc                                                                                                                        <  2:4903076/150‑1 (MQ=255)
                                  cGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCtt                                                                                                                     <  1:2272958/150‑1 (MQ=255)
                                   gAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTc                                                                                                                    <  2:245629/150‑1 (MQ=255)
                                   gAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTc                                                                                                                    >  1:2201276/1‑150 (MQ=255)
                                       aaaGAGAGTGAAGAAGGGCTATCATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAAt                                                                                                                <  1:3298131/150‑1 (MQ=255)
                                       aaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAAt                                                                                                                >  1:2020197/1‑150 (MQ=255)
                                        aaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATc                                                                                                               <  1:4858944/150‑1 (MQ=255)
                                         agagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATcc                                                                                                              >  1:5572962/1‑150 (MQ=255)
                                         agagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATcc                                                                                                              >  1:476094/1‑150 (MQ=255)
                                          gagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCa                                                                                                             <  1:805443/150‑1 (MQ=255)
                                          gagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCa                                                                                                             <  1:6035946/150‑1 (MQ=255)
                                             agTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATa                                                                                                          >  1:1071540/1‑150 (MQ=255)
                                                  agaagGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTg                                                                                                     <  2:6228649/150‑1 (MQ=255)
                                                          tataATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAAc                                                                                             <  1:6458319/150‑1 (MQ=255)
                                                           ataATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACg                                                                                            >  1:4034088/1‑150 (MQ=255)
                                                           ataATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCATCAAATCCAATAGATTGCAGTGAACg                                                                                            >  1:6149841/1‑150 (MQ=255)
                                                             aaTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACgtg                                                                                          >  1:5322854/1‑150 (MQ=255)
                                                              aTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                                                                                         <  2:4544021/150‑1 (MQ=255)
                                                              aTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                                                                                         >  1:4167103/1‑150 (MQ=255)
                                                               tctcTTTTGAACTTTACGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                                                                                        <  1:3404379/150‑1 (MQ=255)
                                                                ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                       <  1:1811213/150‑1 (MQ=255)
                                                                ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                       <  2:2134899/150‑1 (MQ=255)
                                                                ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                       >  1:3353021/1‑150 (MQ=255)
                                                                  ctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGc                                                                                     <  1:1755647/150‑1 (MQ=255)
                                                                   ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                    >  1:4566630/1‑150 (MQ=255)
                                                                   ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                    >  2:1624811/1‑150 (MQ=255)
                                                                   ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                    >  1:4577470/1‑150 (MQ=255)
                                                                    tttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                                                                                   <  1:3939552/150‑1 (MQ=255)
                                                                    tttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                                                                                   <  2:166304/150‑1 (MQ=255)
                                                                    tttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAAACCGATAGATTGCAGTGACCGTGTAGGCCt                                                                                   >  1:5452404/1‑150 (MQ=255)
                                                                     ttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                                                                                  >  2:3366421/1‑150 (MQ=255)
                                                                      tGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGa                                                                                 >  2:4512269/1‑150 (MQ=255)
                                                                      tGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGa                                                                                 <  2:3300710/150‑1 (MQ=255)
                                                                      tGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCAGGCATGA‑‑AAACCCTTCAAATCCAATATATTGCAGTGAACGTGTAGGCCTGa                                                                                 <  1:4456384/150‑1 (MQ=255)
                                                                       gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGTGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                                                                                <  2:3673889/150‑1 (MQ=255)
                                                                       gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                                                                                <  2:660359/150‑1 (MQ=255)
                                                                       gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                                                                                >  2:4180252/1‑150 (MQ=255)
                                                                           tttAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCGTCACAGCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                            >  1:6687644/1‑150 (MQ=255)
                                                                              aaGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                         <  1:2108407/150‑1 (MQ=38)
                                                                                gCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                       >  2:5556199/1‑150 (MQ=35)
                                                                                gCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                       >  2:228561/1‑150 (MQ=35)
                                                                                    aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGACCGTGTAGGCCTGATAAGCGTAGCGCAt                                                                   >  2:2078193/1‑150 (MQ=17)
                                                                                    aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt                                                                   >  2:6032272/1‑150 (MQ=32)
                                                                                     aaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATc                                                                  >  1:5930349/1‑150 (MQ=32)
                                                                                      aaaggCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGAGAACGTGTAGGCCTGATAAGCGTAGCGCATCa                                                                 <  1:2595563/147‑1 (MQ=17)
                                                                                          gcgcTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                             >  1:3667614/1‑150 (MQ=32)
                                                                                          gcgcTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCACGCCGCATCCGGCATGa                                                                                                                                                 >  1:5531690/1‑68 (MQ=11)
                                                                                            gcggTAAACGGCGCCTATGTCAGAGTGTAGACAACGTAGGCCTTGGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAAGAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt                                                           <  1:6257088/147‑1 (MQ=11)
                                                                                             cTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAAt                                                                                                                >  1:1522996/1‑96 (MQ=33)
                                                                                               gTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACgt                                                                                           >  2:3884087/1‑115 (MQ=32)
                                                                                               gTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                        <  2:5772906/150‑1 (MQ=32)
                                                                                                 aaaaGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                      <  2:6184990/150‑1 (MQ=32)
                                                                                                        gcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                <  1:3142141/133‑1 (MQ=31)
                                                                                                        gcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtc                                               <  2:1595819/150‑1 (MQ=31)
                                                                                                         ccATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtca                                              <  1:872598/150‑1 (MQ=11)
                                                                                                                    ttGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaaa                                   <  1:3170844/150‑1 (MQ=33)
                                                                                                                      gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                 >  1:4226744/1‑150 (MQ=34)
                                                                                                                       tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATgg                                <  1:3011033/150‑1 (MQ=34)
                                                                                                                        aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                               <  2:4034088/150‑1 (MQ=34)
                                                                                                                        aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                               <  1:516418/150‑1 (MQ=34)
                                                                                                                         gACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGACCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTGTGTCATCAGTTTCAAATGgcg                              >  2:5176600/1‑150 (MQ=21)
                                                                                                                         gACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                              <  2:1874961/150‑1 (MQ=34)
                                                                                                                          aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                             >  1:1659780/1‑150 (MQ=34)
                                                                                                                           cAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCt                            <  1:4180252/150‑1 (MQ=34)
                                                                                                                                    cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                   <  1:980357/150‑1 (MQ=255)
                                                                                                                                    cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                   <  2:5490078/150‑1 (MQ=255)
                                                                                                                                     tttGTTCATGCCGGATGCGGCGTGAACGCCATATCCGGCAAGA‑‑AAAACCTTCAAAAACAATAGATTGCAGAGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                  >  1:4467896/1‑150 (MQ=14)
                                                                                                                                        gTTCATGCCGGATGCGGCGTGAAAGCCTTATCCGGCATTA‑‑AAACCCATCAAATCCAATAGATTGCAGTGTACGTGTAGGCCTGATAAGCGTATCGCATCAGGCAAAGTTGCGTTTGTAATCAGTTTCAAATGGCGCTGTTAAAAGCggcc               >  1:2873305/1‑147 (MQ=11)
                                                                                                                                            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAGAAATGGCGCTGTAGAAGGGGGCAtttt           >  2:6442021/1‑150 (MQ=11)
                                                                                                                                            aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCGGTTTCAAATGGCGCTGTAAAAGGCGTCAtt             >  2:5457878/1‑150 (MQ=14)
                                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc          <  1:1748931/150‑1 (MQ=255)
                                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCATCAAATCCAAAAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATAAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc          >  1:1004049/1‑150 (MQ=255)
                                                                                                                                                 ggATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCCGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat      >  2:2884815/1‑150 (MQ=255)
                                                                                                                                                 ggATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat      >  2:1206405/1‑150 (MQ=255)
                                                                                                                                                  gATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATAtt     <  2:1181370/150‑1 (MQ=255)
                                                                                                                                                   aTGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg    <  2:3045730/150‑1 (MQ=255)
                                                                                                                                                     gcggcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa  <  1:2760370/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTA  >  NC_000913/1551929‑1552227

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