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

GATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCA  >  NC_000913/3255912‑3256388
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
gATGAAGAGAGATCTTTCGCCAGCAACCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCCCTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGt                                                                                                                                                                                                                                    <  2:15527/251‑1 (MQ=255)
  tGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACGATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTa                                                                                                                                                                                                                                   >  1:329504/1‑250 (MQ=255)
                    cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCg                                                                                                                                                                                                                <  1:192716/251‑1 (MQ=255)
                      gCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGCTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTCACCACCACACCATCGTGCGtc                                                                                                                                                                                                              >  1:204658/1‑251 (MQ=255)
                         aaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCg                                                                                                                                                                                                           <  1:344762/251‑1 (MQ=255)
                              cgcTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGCTAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATAGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTACACACCACACCATCGTGCGTCTCGATATg                                                                                                                                                                                                      >  1:373470/1‑251 (MQ=255)
                                  cGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGcac                                                                                                                                                                                                   <  1:258242/250‑1 (MQ=255)
                                       tGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                                                                                                                                                                                             >  1:132558/1‑251 (MQ=255)
                                           tttCCAGCGGCGCGCCAATATGCAGCCCCCCCTTACCGCTCAAACCTTCCTGCGATAACGCCCAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGGCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGt                                                                                                                                                                                          <  1:172287/250‑1 (MQ=255)
                                            ttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                                                                                                                                                                                         <  1:21308/250‑1 (MQ=255)
                                            ttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTGCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                                                                                                                                                                                         <  1:75888/250‑1 (MQ=255)
                                                 gCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccg                                                                                                                                                                                   <  1:365633/251‑1 (MQ=255)
                                                     cgcgCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgccg                                                                                                                                                                                >  1:149922/1‑250 (MQ=255)
                                                      gcgcCAATATGCAGCCCCCACTTACCGCTCAAACCCTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgac                                                                                                                                                                              >  1:106750/1‑251 (MQ=255)
                                                       cgcCAATCTGCCGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGCATAAAGCGGATCGCCGAAAACGGGACTTCATTGTCGACCTTCAGCATCTCAGCCAGCGTCGCTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGACCTCTGCCACACACCCC‑TGTTGGGTAAACACCACACCATCGTGAGTCTCGATATGCACAATGTTGGTCTGCACGCcgacg                                                                                                                                                                             >  2:165983/1‑251 (MQ=255)
                                                       cgcCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacg                                                                                                                                                                             <  1:80353/251‑1 (MQ=255)
                                                             tatGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTg                                                                                                                                                                       >  1:165172/1‑251 (MQ=255)
                                                             tatGCACCCCCCACTTACCGCTCAACCCTTCATGCGCTAACGCACAATTTACCTTCGCGGAATCGAGAATAAAGCGCATCGCCGCAAACGGCACTTCATTGACGCACTTCAGGATCTGAGCCAGCGTCGTTCTGGATAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCGCAAACGCC‑TGTTGCGCAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACTATGGTg                                                                                                                                                                       >  2:107119/1‑251 (MQ=255)
                                                                        cACTTACCGCTCAAACCTTCCTGAGATAACGCACAACTTAGCTTCGCGGAATCGAGAATAAATCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACAATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGGCACACGccc                                                                                                                                                             >  2:238947/1‑250 (MQ=255)
                                                                         aCTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCACAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTCTGCCCGCCGCCGATGGCGCCCCACGCCCAc                                                                                                                                                           >  2:43602/1‑251 (MQ=255)
                                                                            tacagCTCAAACCTTCCAGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGCCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGAAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGCCGATGGTGACACACGCCCACTtc                                                                                                                                                        >  2:227609/5‑251 (MQ=255)
                                                                            tACCGCTCAAACCTGCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACtt                                                                                                                                                         <  1:294314/250‑1 (MQ=255)
                                                                             accgctacaaCCTTCCTTCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGCCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGTCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACCCACGCC‑TGTTTGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTtc                                                                                                                                                        <  1:81178/242‑1 (MQ=255)
                                                                                gCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCTTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCAc                                                                                                                                                    <  1:141495/251‑1 (MQ=255)
                                                                                      actttccggctatatcGCACAAAATAGCTTAGCGGAACCCAAAATAAAGAGGATCGCCGCAAACAGGACTACATTGACGAACTTCAGGATCTCTTCCCGCGGCGTTCTGGAAAGCACCGTCAGCGGAGACTTTTGCTCGACCGCTCCCACCCCCGCC‑TGTTGGTTAATCACCACACCAGCGTGCGCCTCCATATGCACAATTTTGTTATGCCCGCCGACGATGGTGACACACGCCAACTTCTCACCGTTcc                                                                                                                                              <  2:401759/236‑1 (MQ=255)
                                                                                         ttCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGCTGGTGACACACGCCCACTTCTCACCGTTCCAgc                                                                                                                                           >  2:96294/1‑250 (MQ=255)
                                                                                            cTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCACCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACtt                                                                                                                                        >  1:138161/1‑251 (MQ=255)
                                                                                            cTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATAGCCGCAACAGGGACTTCATTGACGAACGTCAGAATCTCAGCCCGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACCCGCC‑TGTTGGGCACACACCACACCATCCTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGCCGATAGTGACACACGCCCACTTCTCACCGTTCCAgcatt                                                                                                                                        >  2:354072/1‑247 (MQ=255)
                                                                                               cGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCGCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTgg                                                                                                                                     <  1:240403/251‑1 (MQ=255)
                                                                                                 aTAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGgcg                                                                                                                                   <  1:269284/251‑1 (MQ=255)
                                                                                                   aaCGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGcgcg                                                                                                                                 >  1:310773/1‑251 (MQ=255)
                                                                                                       cacaATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                               >  2:144173/1‑249 (MQ=255)
                                                                                                          aaTTTAGCTTCGCGGAATCGAGACTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTgaga                                                                                                                            >  2:368079/1‑249 (MQ=255)
                                                                                                              tAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATATCAGCCAGCGTCGTTCTTGAAAGCACCGTCAGCGGAGACTCTTACTCGCCCTCTGCCACACACGCC‑TGTTGGGTCAACACCACACCATCGTGCGTCTCGATCTGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCCTTTGGCGCGAGAGAAGAggc                                                                                                                      >  1:180919/1‑250 (MQ=255)
                                                                                                                 cTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCCCAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCCTTTGGCGCGTGAGAAGAGGAttt                                                                                                                   >  1:199630/1‑251 (MQ=255)
                                                                                                                   tCGCGGAATCGAGCATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGCACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGAAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCACCGACGATGGGGACACACGCCCACTTCTCACCGTTCCCGACTTTGGCGCGTGAGAAGAGGATTTc                                                                                                                  >  1:339027/1‑250 (MQ=255)
                                                                                                                    cgcgGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCAt                                                                                                                >  1:100237/1‑251 (MQ=255)
                                                                                                                      gggAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCg                                                                                                              <  2:320057/250‑1 (MQ=255)
                                                                                                                        gAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGc                                                                                                             >  2:298127/1‑250 (MQ=255)
                                                                                                                         aaTCGAGAATAAAGCGGCTCGCCGCAAATGGGCCTTCATTGTCGAACTTAAGGATCTCAACCAGCGTCCTTCTGGAAAGCACCGACAGCGGAGACTCTTGCACGCCCTCTGCAACACACGCC‑TGTTGGGTAAACAACACACCATCGTGCATCTCGATATGCACCATGTTGGAATGCCCGACGACGATGGTGCCACACGCCCACTTCTCACCGTTCCAGACTTTGACGCGTGAGAAGCGGATTTCATCGCa                                                                                                            >  2:338761/1‑250 (MQ=255)
                                                                                                                                 aTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTg                                                                                                   >  2:44513/1‑251 (MQ=255)
                                                                                                                                   aaaGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAATGTTCCTgg                                                                                                  <  1:92573/250‑1 (MQ=255)
                                                                                                                                    aaGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCCGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGAt                                                                                                >  2:203847/1‑251 (MQ=255)
                                                                                                                                            cgccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGATCCTGGATCTTAACgg                                                                                        >  1:324630/1‑251 (MQ=255)
                                                                                                                                                 cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACt                                                                                   >  1:296541/1‑251 (MQ=255)
                                                                                                                                                        gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccg                                                                            >  1:357166/1‑251 (MQ=255)
                                                                                                                                                          cTTCATTGACGAACTTCAAGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCCCAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTAGCGCGTGAGAAGAGGATTTCCTCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgcc                                                                          >  2:143820/1‑251 (MQ=255)
                                                                                                                                                           ttCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgccg                                                                         >  1:37238/1‑251 (MQ=255)
                                                                                                                                                             cATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCgccgcc                                                                       >  2:383847/1‑251 (MQ=255)
                                                                                                                                                               gtGACGAACTCCAGGAGCCAAGCCAGCGTCGTTCTGGAAAGCACCGCCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑GGTGGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcag                                                                     <  1:144173/250‑1 (MQ=255)
                                                                                                                                                                  ccGAACTTCAGGTTCTCAGCCAGCGACGTTCTGGAAAGCTCCGGCAGCGGAGACTCTTGCGCTCCCTCTGCCTCTCACGCC‑TGTTGGGTAAACCCCACACCAGCGTGCGTCGCGATATGCACAATGTGGGTATGCCCGCCGACGATGGTGTCACACGCCCACTTCTCGCCGTCCCAGTCGTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACGGAGACTTGCCCCGCCGCcagcag                                                                  <  1:257582/250‑1 (MQ=255)
                                                                                                                                                                                   agccagctgctgtCTGGAAACCACCGTCGGCGGAGCCTCTTGCGCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACCCCACACCATCGGGCGTCTCGATATGCAAAATGTTGGTAGGCCCGCCGACGATGGGGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGTTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCaa                                                 <  2:347838/239‑1 (MQ=255)
                                                                                                                                                                                         cgtcgtCCTGGAAAGCCCCGGCAGCGGGGACTCTTGCTCGCCCTCTGCCACCCACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGTTCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCt                                           <  2:57086/251‑1 (MQ=255)
                                                                                                                                                                                          gtcgtTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCACGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTg                                          >  1:27052/1‑251 (MQ=255)
                                                                                                                                                                                            cgtccTGGAAAGCACCGTCAGCGGAGACTCTTGCCCGCCCTCTGCCCCACCCGCC‑TGTTGGGTAAACCCCACCCCTTCGTGCGTCCCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTgc                                         <  2:106750/246‑1 (MQ=255)
                                                                                                                                                                                                        ccccttCAGCGGGGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCCAAATGTTGGAATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGTTTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCtt                            <  1:192219/246‑1 (MQ=255)
                                                                                                                                                                                                        ccccGTCAGCGGAGACCCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAAGGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCTTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCt                                      <  1:36945/239‑1 (MQ=255)
                                                                                                                                                                                                        cACCGTCAGCGGAGACTCTTGCTCGCACTCTGCCACACACGCC‑TGTTGGGTCAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCCGTGCTTTGGCATCGGCAATTGCCTGCGCt                                      >  2:36945/1‑241 (MQ=255)
                                                                                                                                                                                                           cGTCAGCGGAGCCTCTTGCTCGCCCTTTGCCACCCACGCC‑GGTTGGGTAAACACCACACAATCGTGCGTCTGGATATGCACAATGTTGGAATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGGGAGAAGAGGATTTCATCGCAAGTTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCa                         <  2:143859/251‑1 (MQ=255)
                                                                                                                                                                                                           cGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                          >  1:153533/1‑250 (MQ=255)
                                                                                                                                                                                                              cAGCGGAAACTCTTGCTAGCCCTCTGCCACACCCGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCACGCCGCCAGCAGTGCTTTGGCATCGGCACTTGCCTGCGCTGTTGCGTCTTTCAGc                       >  1:286844/1‑250 (MQ=255)
                                                                                                                                                                                                                   gagtctcttccgcGCCCCCTGCCACACACCCC‑TTTTGGGGAACCACCACACCAACGTGCGTCTCGATATGCACAATGTGGGGAGGCCCGCCGACGATGGGGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTACCGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                 <  2:132558/239‑1 (MQ=255)
                                                                                                                                                                                                                     gaCTCTTGCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                 >  2:185847/1‑249 (MQ=255)
                                                                                                                                                                                                                     gaCTCTTGCCCTCCCTCTGCCACACACGCC‑TTTTGGTTAACCCCCACACCACCGTGCGTCTCGATATGCACATGGTGGGTATGCCCGCCGACGATGGTGACACACGCCCGCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGATTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAg               <  1:282961/251‑1 (MQ=255)
                                                                                                                                                                                                                      accacgtctcacccTCTGCCCCACGCCCC‑TTTTGGGTAAACCCCACCCCATCGGGCTACTCTATTTGCTAAAGGTTGGTATGCCCGGCGACGAGGGGGACACACGCCCCCGTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCatc              <  2:51556/240‑3 (MQ=255)
                                                                                                                                                                                                                         cgtGCTCGCCCTCTGCCACACACTCC‑TGTTGGGTAAACACCACCCCATCGTGCGGCTCGATATGCAAAATGTTGGTATGCCCGCCGGCGATGGTGACACACGCCCACTTCCCACCGTTCCAGACTGTGGCGCGGGAGAATAGGATTTCATCTCAAGGTTCCTGGATCTTACCGGATACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc            <  1:250270/248‑1 (MQ=255)
                                                                                                                                                                                                                           tGCTCCCCCTCGGCCACACACGCC‑TGTTGGGTAAACCCCACCCCATCGTGCGTCTAGATAGGCAAAATGTTGGTATGCCCGCCGACGATGGGGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCTCTTCCAGCCCg           <  2:288143/249‑1 (MQ=255)
                                                                                                                                                                                                                            gCTCGCCCTCTGCCACACACGCC‑TGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCg        >  1:96746/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCTCCCCCTCTTCCACACTCGCC‑TGTTGGGTAAACCCCACCCCATCGTTCGTCTCGATATGCACAATTTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAAAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGc         <  2:26256/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ccTCTGCCACACACGCCGTG‑TGGTTAAACACCACCCCAGCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACGTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACGTCCAGCCCGTCGTTGGCa  <  2:248072/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
GATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCC‑TGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCA  >  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: ‑