Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,256,157 C→T Q165Q (CAG→CAA cyuA ← putative L‑cysteine desulfidase CyuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,1570CT98.6% 227.1 / ‑4.8 70coding (495/1311 nt)cyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/31);  minor base . (0/1);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCA  >  NC_000913/3255912‑3256388
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
gATGAAGAGAGATCTTTCGCCAGCAACCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCCCTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGt                                                                                                                                                                                                                                    <  2:15527/251‑1 (MQ=255)
  tGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACGATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTa                                                                                                                                                                                                                                   >  1:329504/1‑250 (MQ=255)
                    cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCg                                                                                                                                                                                                                <  1:192716/251‑1 (MQ=255)
                      gCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGCTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTCACCACCACACCATCGTGCGtc                                                                                                                                                                                                              >  1:204658/1‑251 (MQ=255)
                         aaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCg                                                                                                                                                                                                           <  1:344762/251‑1 (MQ=255)
                              cgcTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGCTAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATAGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTACACACCACACCATCGTGCGTCTCGATATg                                                                                                                                                                                                      >  1:373470/1‑251 (MQ=255)
                                  cGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGcac                                                                                                                                                                                                   <  1:258242/250‑1 (MQ=255)
                                       tGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                                                                                                                                                                                             >  1:132558/1‑251 (MQ=255)
                                           tttCCAGCGGCGCGCCAATATGCAGCCCCCCCTTACCGCTCAAACCTTCCTGCGATAACGCCCAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGGCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGt                                                                                                                                                                                          <  1:172287/250‑1 (MQ=255)
                                            ttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                                                                                                                                                                                         <  1:21308/250‑1 (MQ=255)
                                            ttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTGCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                                                                                                                                                                                         <  1:75888/250‑1 (MQ=255)
                                                 gCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccg                                                                                                                                                                                   <  1:365633/251‑1 (MQ=255)
                                                     cgcgCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgccg                                                                                                                                                                                >  1:149922/1‑250 (MQ=255)
                                                      gcgcCAATATGCAGCCCCCACTTACCGCTCAAACCCTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgac                                                                                                                                                                              >  1:106750/1‑251 (MQ=255)
                                                       cgcCAATCTGCCGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGCATAAAGCGGATCGCCGAAAACGGGACTTCATTGTCGACCTTCAGCATCTCAGCCAGCGTCGCTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGACCTCTGCCACACACCCCTGTTGGGTAAACACCACACCATCGTGAGTCTCGATATGCACAATGTTGGTCTGCACGCcgacg                                                                                                                                                                             >  2:165983/1‑251 (MQ=255)
                                                       cgcCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacg                                                                                                                                                                             <  1:80353/251‑1 (MQ=255)
                                                             tatGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTg                                                                                                                                                                       >  1:165172/1‑251 (MQ=255)
                                                             tatGCACCCCCCACTTACCGCTCAACCCTTCATGCGCTAACGCACAATTTACCTTCGCGGAATCGAGAATAAAGCGCATCGCCGCAAACGGCACTTCATTGACGCACTTCAGGATCTGAGCCAGCGTCGTTCTGGATAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCGCAAACGCCTGTTGCGCAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACTATGGTg                                                                                                                                                                       >  2:107119/1‑251 (MQ=255)
                                                                        cACTTACCGCTCAAACCTTCCTGAGATAACGCACAACTTAGCTTCGCGGAATCGAGAATAAATCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACAATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGGCACACGccc                                                                                                                                                             >  2:238947/1‑250 (MQ=255)
                                                                         aCTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCACAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTCTGCCCGCCGCCGATGGCGCCCCACGCCCAc                                                                                                                                                           >  2:43602/1‑251 (MQ=255)
                                                                            tacagCTCAAACCTTCCAGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGCCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGAAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGCCGATGGTGACACACGCCCACTtc                                                                                                                                                        >  2:227609/5‑251 (MQ=255)
                                                                            tACCGCTCAAACCTGCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACtt                                                                                                                                                         <  1:294314/250‑1 (MQ=255)
                                                                             accgctacaaCCTTCCTTCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGCCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGTCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACCCACGCCTGTTTGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTtc                                                                                                                                                        <  1:81178/242‑1 (MQ=255)
                                                                                gCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCTTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCAc                                                                                                                                                    <  1:141495/251‑1 (MQ=255)
                                                                                      actttccggctatatcGCACAAAATAGCTTAGCGGAACCCAAAATAAAGAGGATCGCCGCAAACAGGACTACATTGACGAACTTCAGGATCTCTTCCCGCGGCGTTCTGGAAAGCACCGTCAGCGGAGACTTTTGCTCGACCGCTCCCACCCCCGCCTGTTGGTTAATCACCACACCAGCGTGCGCCTCCATATGCACAATTTTGTTATGCCCGCCGACGATGGTGACACACGCCAACTTCTCACCGTTcc                                                                                                                                              <  2:401759/236‑1 (MQ=255)
                                                                                         ttCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGCTGGTGACACACGCCCACTTCTCACCGTTCCAgc                                                                                                                                           >  2:96294/1‑250 (MQ=255)
                                                                                            cTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCACCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACtt                                                                                                                                        >  1:138161/1‑251 (MQ=255)
                                                                                            cTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATAGCCGCAACAGGGACTTCATTGACGAACGTCAGAATCTCAGCCCGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACCCGCCTGTTGGGCACACACCACACCATCCTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGCCGATAGTGACACACGCCCACTTCTCACCGTTCCAgcatt                                                                                                                                        >  2:354072/1‑247 (MQ=255)
                                                                                               cGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCGCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTgg                                                                                                                                     <  1:240403/251‑1 (MQ=255)
                                                                                                 aTAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGgcg                                                                                                                                   <  1:269284/251‑1 (MQ=255)
                                                                                                   aaCGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGcgcg                                                                                                                                 >  1:310773/1‑251 (MQ=255)
                                                                                                       cacaATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                               >  2:144173/1‑249 (MQ=255)
                                                                                                          aaTTTAGCTTCGCGGAATCGAGACTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTgaga                                                                                                                            >  2:368079/1‑249 (MQ=255)
                                                                                                              tAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATATCAGCCAGCGTCGTTCTTGAAAGCACCGTCAGCGGAGACTCTTACTCGCCCTCTGCCACACACGCCTGTTGGGTCAACACCACACCATCGTGCGTCTCGATCTGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCCTTTGGCGCGAGAGAAGAggc                                                                                                                      >  1:180919/1‑250 (MQ=255)
                                                                                                                 cTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCCCAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCCTTTGGCGCGTGAGAAGAGGAttt                                                                                                                   >  1:199630/1‑251 (MQ=255)
                                                                                                                   tCGCGGAATCGAGCATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGCACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGAAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCACCGACGATGGGGACACACGCCCACTTCTCACCGTTCCCGACTTTGGCGCGTGAGAAGAGGATTTc                                                                                                                  >  1:339027/1‑250 (MQ=255)
                                                                                                                    cgcgGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCAt                                                                                                                >  1:100237/1‑251 (MQ=255)
                                                                                                                      gggAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCg                                                                                                              <  2:320057/250‑1 (MQ=255)
                                                                                                                        gAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGc                                                                                                             >  2:298127/1‑250 (MQ=255)
                                                                                                                         aaTCGAGAATAAAGCGGCTCGCCGCAAATGGGCCTTCATTGTCGAACTTAAGGATCTCAACCAGCGTCCTTCTGGAAAGCACCGACAGCGGAGACTCTTGCACGCCCTCTGCAACACACGCCTGTTGGGTAAACAACACACCATCGTGCATCTCGATATGCACCATGTTGGAATGCCCGACGACGATGGTGCCACACGCCCACTTCTCACCGTTCCAGACTTTGACGCGTGAGAAGCGGATTTCATCGCa                                                                                                            >  2:338761/1‑250 (MQ=255)
                                                                                                                                 aTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTg                                                                                                   >  2:44513/1‑251 (MQ=255)
                                                                                                                                   aaaGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAATGTTCCTgg                                                                                                  <  1:92573/250‑1 (MQ=255)
                                                                                                                                    aaGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCCGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGAt                                                                                                >  2:203847/1‑251 (MQ=255)
                                                                                                                                            cgccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGATCCTGGATCTTAACgg                                                                                        >  1:324630/1‑251 (MQ=255)
                                                                                                                                                 cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACt                                                                                   >  1:296541/1‑251 (MQ=255)
                                                                                                                                                        gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccg                                                                            >  1:357166/1‑251 (MQ=255)
                                                                                                                                                          cTTCATTGACGAACTTCAAGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCCCAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTAGCGCGTGAGAAGAGGATTTCCTCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgcc                                                                          >  2:143820/1‑251 (MQ=255)
                                                                                                                                                           ttCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgccg                                                                         >  1:37238/1‑251 (MQ=255)
                                                                                                                                                             cATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCgccgcc                                                                       >  2:383847/1‑251 (MQ=255)
                                                                                                                                                               gtGACGAACTCCAGGAGCCAAGCCAGCGTCGTTCTGGAAAGCACCGCCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCGGTGGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcag                                                                     <  1:144173/250‑1 (MQ=255)
                                                                                                                                                                  ccGAACTTCAGGTTCTCAGCCAGCGACGTTCTGGAAAGCTCCGGCAGCGGAGACTCTTGCGCTCCCTCTGCCTCTCACGCCTGTTGGGTAAACCCCACACCAGCGTGCGTCGCGATATGCACAATGTGGGTATGCCCGCCGACGATGGTGTCACACGCCCACTTCTCGCCGTCCCAGTCGTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACGGAGACTTGCCCCGCCGCcagcag                                                                  <  1:257582/250‑1 (MQ=255)
                                                                                                                                                                                   agccagctgctgtCTGGAAACCACCGTCGGCGGAGCCTCTTGCGCGCCCTCTGCCACACACGCCTGTTGGGTAAACCCCACACCATCGGGCGTCTCGATATGCAAAATGTTGGTAGGCCCGCCGACGATGGGGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGTTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCaa                                                 <  2:347838/239‑1 (MQ=255)
                                                                                                                                                                                         cgtcgtCCTGGAAAGCCCCGGCAGCGGGGACTCTTGCTCGCCCTCTGCCACCCACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGTTCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCt                                           <  2:57086/251‑1 (MQ=255)
                                                                                                                                                                                          gtcgtTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCACGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTg                                          >  1:27052/1‑251 (MQ=255)
                                                                                                                                                                                            cgtccTGGAAAGCACCGTCAGCGGAGACTCTTGCCCGCCCTCTGCCCCACCCGCCTGTTGGGTAAACCCCACCCCTTCGTGCGTCCCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTgc                                         <  2:106750/246‑1 (MQ=255)
                                                                                                                                                                                                        ccccttCAGCGGGGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCCAAATGTTGGAATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGTTTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCtt                            <  1:192219/246‑1 (MQ=255)
                                                                                                                                                                                                        ccccGTCAGCGGAGACCCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAAGGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCTTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCt                                      <  1:36945/239‑1 (MQ=255)
                                                                                                                                                                                                        cACCGTCAGCGGAGACTCTTGCTCGCACTCTGCCACACACGCCTGTTGGGTCAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCCGTGCTTTGGCATCGGCAATTGCCTGCGCt                                      >  2:36945/1‑241 (MQ=255)
                                                                                                                                                                                                           cGTCAGCGGAGCCTCTTGCTCGCCCTTTGCCACCCACGCCGGTTGGGTAAACACCACACAATCGTGCGTCTGGATATGCACAATGTTGGAATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGGGAGAAGAGGATTTCATCGCAAGTTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCa                         <  2:143859/251‑1 (MQ=255)
                                                                                                                                                                                                           cGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                          >  1:153533/1‑250 (MQ=255)
                                                                                                                                                                                                              cAGCGGAAACTCTTGCTAGCCCTCTGCCACACCCGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCACGCCGCCAGCAGTGCTTTGGCATCGGCACTTGCCTGCGCTGTTGCGTCTTTCAGc                       >  1:286844/1‑250 (MQ=255)
                                                                                                                                                                                                                   gagtctcttccgcGCCCCCTGCCACACACCCCTTTTGGGGAACCACCACACCAACGTGCGTCTCGATATGCACAATGTGGGGAGGCCCGCCGACGATGGGGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTACCGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                 <  2:132558/239‑1 (MQ=255)
                                                                                                                                                                                                                     gaCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                 >  2:185847/1‑249 (MQ=255)
                                                                                                                                                                                                                     gaCTCTTGCCCTCCCTCTGCCACACACGCCTTTTGGTTAACCCCCACACCACCGTGCGTCTCGATATGCACATGGTGGGTATGCCCGCCGACGATGGTGACACACGCCCGCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGATTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAg               <  1:282961/251‑1 (MQ=255)
                                                                                                                                                                                                                      accacgtctcacccTCTGCCCCACGCCCCTTTTGGGTAAACCCCACCCCATCGGGCTACTCTATTTGCTAAAGGTTGGTATGCCCGGCGACGAGGGGGACACACGCCCCCGTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCatc              <  2:51556/240‑3 (MQ=255)
                                                                                                                                                                                                                         cgtGCTCGCCCTCTGCCACACACTCCTGTTGGGTAAACACCACCCCATCGTGCGGCTCGATATGCAAAATGTTGGTATGCCCGCCGGCGATGGTGACACACGCCCACTTCCCACCGTTCCAGACTGTGGCGCGGGAGAATAGGATTTCATCTCAAGGTTCCTGGATCTTACCGGATACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc            <  1:250270/248‑1 (MQ=255)
                                                                                                                                                                                                                           tGCTCCCCCTCGGCCACACACGCCTGTTGGGTAAACCCCACCCCATCGTGCGTCTAGATAGGCAAAATGTTGGTATGCCCGCCGACGATGGGGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCTCTTCCAGCCCg           <  2:288143/249‑1 (MQ=255)
                                                                                                                                                                                                                            gCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCg        >  1:96746/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCTCCCCCTCTTCCACACTCGCCTGTTGGGTAAACCCCACCCCATCGTTCGTCTCGATATGCACAATTTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAAAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGc         <  2:26256/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ccTCTGCCACACACGCCGTGTGGTTAAACACCACCCCAGCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACGTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACGTCCAGCCCGTCGTTGGCa  <  2:248072/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
GATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCA  >  NC_000913/3255912‑3256388

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