Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,547,239 T→A 6.3% intergenic (‑70/‑162) yddG ← / → fdnG aromatic amino acid exporter/formate dehydrogenase N subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,547,2390TA 6.3% 216.6 / 4.4 63intergenic (‑70/‑162)yddG/fdnGaromatic amino acid exporter/formate dehydrogenase N subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GGCGAGTAAATAGCCTTTCGGGATTTGCCGAATACGCGGAAATCCAACCGTGAAGATTAACAGCAGCCCGCTTAATGAATAGATAGCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAG  >  NC_000913/1546969‑1547461
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                               
cggaggtaaaagcccTTTCGGGATTTACCGAAAACGCGGAAATCCAACCGGGGAGATTAACAGCAGCCCGCTTAATGAAACGATGGCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTc                                                                                                                                                                                                                                <  2:451891/259‑1 (MQ=255)
                     gATTTGCCGAATACGCGGAAATCCAACCGTGAAGATTAACAGCAGCCCGCTTAATGAATAGATAGCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTc                                                                                                                                                                                                             <  1:298323/270‑1 (MQ=255)
                        ttGCCGAATACGCGGAAATCCAACCGTGAAGATTAACAGCAGCCCGCTTAATGAATAGATAGCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCAGCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCg                                                                                                                                                                                                        <  1:468985/272‑1 (MQ=255)
                                                               cagcCCGCTTAATGAATAGATAGCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGGTCGTTTTTCTACCGCTATTGAGGTAGGTTCAATTTCGAAgggcggtgattttt                                                                                                                                                                >  2:287279/1‑260 (MQ=255)
                                                                          aTGAATAGATAGCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCCCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCCGTGGTGCGTTTTTCTACCGCCATTGTGGTAGGTGTATTTGATGAGGCGGCTCATTTTTTggtcaataaa                                                                                                                                                     >  2:531156/1‑264 (MQ=255)
                                                                                  atagCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGc                                                                                                                                                                                                                            >  1:279483/1‑194 (MQ=255)
                                                                                       agcTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTTCGTTTTTTTACCGCTATTTAGGTAGGTCAATTTGc                                                                                                                                                                                 >  2:785784/1‑231 (MQ=255)
                                                                                        gcTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCg                                                                                                                                                                                <  1:785785/231‑1 (MQ=255)
                                                                                          tGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGttttt                                                                                                                                                                                                              >  1:740040/1‑200 (MQ=255)
                                                                                           gcgcCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTc                                                                                                                                                                                                             <  2:740039/200‑1 (MQ=255)
                                                                                                                             cacaaatcactcCTCCCCTCTTGTGCCACAGGACGTTCGTGTTCAGCCCTTTGAGCTTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTcc                                                                                                  <  2:253413/268‑1 (MQ=255)
                                                                                                                              gcgAATCAATCCTACCATCGTGCTCCACAGGACGTTCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTcc                                                                                                  <  1:287280/271‑1 (MQ=255)
                                                                                                                               cgAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTTAGGCAGTGGTGCGTTTTTCTACCGCTATTTAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTAATGTTAGTTTCCt                                                                                                 >  1:743864/1‑271 (MQ=255)
                                                                                                                                  atcaatcCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGttttt                                                                                                                                                                                                              >  2:58138/1‑160 (MQ=255)
                                                                                                                                  atcaatcCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGGTGTTCTTTTTGATTTCACGCAAAAAAATTTAAAATTTTAATTTTAGTTTTCtt                                                                                                >  2:461391/1‑268 (MQ=255)
                                                                                                                                  atcaatcCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTc                                                                                                   >  1:253413/1‑266 (MQ=255)
                                                                                                                                   tcaatcCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTc                                                                                                                                                                                                             <  1:58138/160‑1 (MQ=255)
                                                                                                                                          tACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAgtgg                                                                                                                                                                                                                       >  2:724658/1‑143 (MQ=255)
                                                                                                                                          tACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGGGCAAAAAATTTCAGCAATTTAATGTTAGTTTCCTCCCTCCGTaaaa                                                                                     >  2:1070381/1‑270 (MQ=255)
                                                                                                                                           aCCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAgtggt                                                                                                                                                                                                                      <  1:724659/143‑1 (MQ=255)
                                                                                                                                                  tGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTGTTTCTACCGCTATTGAGGCAGGTCAATTTGCGATGGCGGATTATTTTTTGGCAAACAGATGTTCTTTTTGATTTCGCGCCAAAAGATTACGAATTTTATTTTTAATTTTCTCTGGCAGTAAAAcaccgtaa                                                                             >  1:1037231/1‑268 (MQ=255)
                                                                                                                                                       cacaGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAATCGGTGTTCTTTTTGATTTCGCGCAAAAAATTTCAGAATTTTAATTTTAGTTTTCTCCGCCAGTACAACCCCTTaaaaaaa                                                                        >  2:494228/1‑271 (MQ=255)
                                                                                                                                                         caGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAAtt                                                                                                                                                                                    >  1:418693/1‑162 (MQ=255)
                                                                                                                                                          aGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAAttt                                                                                                                                                                                   <  2:418692/162‑1 (MQ=255)
                                                                                                                                                            gACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCg                                                                                                                                                                                                        >  2:975326/1‑140 (MQ=255)
                                                                                                                                                            gACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAgg                                                                                                                                                                            >  2:427979/1‑167 (MQ=255)
                                                                                                                                                             aCGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGc                                                                                                                                                                                                       <  1:975327/140‑1 (MQ=255)
                                                                                                                                                             aCGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGc                                                                                                                                                                           <  1:427980/167‑1 (MQ=255)
                                                                                                                                                               gATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCt                                                                                                 >  1:123579/1‑239 (MQ=255)
                                                                                                                                                                    cTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTtgt                                                                                                                                                               >  2:820015/1‑172 (MQ=255)
                                                                                                                                                                    cTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTc                                                                                                   >  2:119864/1‑232 (MQ=255)
                                                                                                                                                                    cTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTc                                                                                                   >  2:1016943/1‑232 (MQ=255)
                                                                                                                                                                    cTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGgaaagaaaag                                                            >  1:1055859/1‑269 (MQ=255)
                                                                                                                                                                     tATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTtgtg                                                                                                                                                              <  1:820016/172‑1 (MQ=255)
                                                                                                                                                                     tATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTcc                                                                                                  <  1:1016944/232‑1 (MQ=255)
                                                                                                                                                                     tATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTcc                                                                                                  <  1:119864/232‑1 (MQ=255)
                                                                                                                                                                         aGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGtt                                                                                                                                                                                                                 >  2:214102/1‑118 (MQ=255)
                                                                                                                                                                          gCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGttt                                                                                                                                                                                                                <  1:214102/118‑1 (MQ=255)
                                                                                                                                                                              tATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTTGAAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTTAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATAACCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCg                                                 >  1:141190/1‑272 (MQ=255)
                                                                                                                                                                              tATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTAAATTTGCGAAGGCGGATTAttt                                                                                                                                                                  >  1:740000/1‑159 (MQ=255)
                                                                                                                                                                               aTGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTAtttt                                                                                                                                                                 <  2:739999/159‑1 (MQ=255)
                                                                                                                                                                                  aGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGt                                                                                                                                                                                                                                 >  2:233177/1‑93 (MQ=255)
                                                                                                                                                                                   gctttGTTTTTTTTCGTGTCTTGTTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTGGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGAc                                            <  1:590374/269‑1 (MQ=255)
                                                                                                                                                                                   gCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTc                                                                                                                                                                                                                                <  1:233177/93‑1 (MQ=255)
                                                                                                                                                                                        gCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACa                                                                                                                                                       >  2:1033954/1‑160 (MQ=255)
                                                                                                                                                                                         cTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAg                                                                                                                                                      <  1:1033955/160‑1 (MQ=255)
                                                                                                                                                                                                   tgtCATGCTCGTTGTTTTGTCTCTCTTTCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCTTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCt                             <  1:531157/271‑1 (MQ=255)
                                                                                                                                                                                                    tttaTGCTCGCTGTTTTGTCTTTCTTTCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTTCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTgc                           <  1:1070382/269‑1 (MQ=255)
                                                                                                                                                                                                                  tttGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTc                                                                                                                                                                                                             >  1:824878/1‑81 (MQ=255)
                                                                                                                                                                                                                   ttGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCt                                                                                                                                                                                                            <  2:824877/81‑1 (MQ=255)
                                                                                                                                                                                                                      tctctcTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTAcc                                                                                                                                                                                                         >  2:922571/1‑81 (MQ=255)
                                                                                                                                                                                                                       ctctctTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCg                                                                                                                                                                                                        <  1:922572/81‑1 (MQ=255)
                                                                                                                                                                                                                        tttctTGCTGTTAAAAATTATGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTTGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAATAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaac       <  1:489226/270‑1 (MQ=255)
                                                                                                                                                                                                                            ttGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGAAAATTTTTTAAAAACAGCACGGGAGGTATGGCTgggaaaaaa   >  2:489225/1‑265 (MQ=255)
                                                                                                                                                                                                                             ttccGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCGTTGTGATTTCGCGCAAAAAGATGCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAg  <  1:494229/270‑1 (MQ=255)
                                                                                                                                                                                                                              gCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGACAGGCAGTGGTTCGTTTTTCTACCGCTATTGAGGTAGGTCAAtt                                                                                                                                                                                    >  2:636779/1‑93 (MQ=255)
                                                                                                                                                                                                                               ccGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGACAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAAttt                                                                                                                                                                                   <  1:636780/93‑1 (MQ=255)
                                                                                                                                                                                                                                   tAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACa                                                                                                                                                       >  2:544056/1‑117 (MQ=255)
                                                                                                                                                                                                                                    aaaaaTTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAg                                                                                                                                                      <  1:544057/117‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTAAGCTGAATTTTATAGCATTTTTTTAACGGGCCTGTCAGGCAGTTGTGCGTTTTGCTACCGCTATTGAGGTAGGTCTATTGGCGAATGCGGATTATTTTTGGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTc                                                                                                <  1:422070/167‑1 (MQ=255)
                                                                                                                                                                                                                                               tGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTAttt                                                                                                                                                                  >  2:495677/1‑94 (MQ=255)
                                                                                                                                                                                                                                                gAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTAtttt                                                                                                                                                                 <  1:495678/94‑1 (MQ=255)
                                                                                                                                                                                                                                                           cATTTTTTTAACTGGCCTGACAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGaaaaa                                                                          >  1:584491/1‑170 (MQ=255)
                                                                                                                                                                                                                                                            aTTTTTTTAACTGGCCTGACAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGaaaaaa                                                                         <  2:584490/170‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                               
GGCGAGTAAATAGCCTTTCGGGATTTGCCGAATACGCGGAAATCCAACCGTGAAGATTAACAGCAGCCCGCTTAATGAATAGATAGCAGCTGCGCCGCCGACCGGGCCGAGCCCCTCACTGACACCGCGAATCAATCCTACCATCGTGCTCCACAGGACGATCGCTATCAGCCCTATGAGCGTTGCTTTTTGTCGTGTCATGCTCGCTGTTTTGTCTCTCTTGCCGTTAAAAATTAAGCTGAATTTTATAGCATTTTTTTAACTGGCCTGTCAGGCAGTGGTGCGTTTTTCTACCGCTATTGAGGTAGGTCAATTTGCGAAGGCGGATTATTTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAG  >  NC_000913/1546969‑1547461

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