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.8% 266.2 / ‑5.2 82coding (495/1311 nt)cyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (50/31);  minor base . (0/1);  total (50/32)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACG  >  NC_000913/3255919‑3256401
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
agagaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGCTAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAACCGGGACTTCATTGACGACCTTCAGGATCTCAGCCAGCGTCGTTCTGACAAGCACCGTCAGCGGATCATCTTGCTCGCCCTCTACCACCCACGCCTGTTGGTTAAacacc                                                                                                                                                                                                                                          >  1:456844/1‑251 (MQ=255)
agagaCCTCTCGCCAGCAAACCGCGCTCGCCCTGTTTTTTCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCGTCCTGCGAGAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACTCCTGTGGGGTAAacacc                                                                                                                                                                                                                                          <  2:521175/251‑1 (MQ=255)
   gaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAacaccac                                                                                                                                                                                                                                        >  1:50766/1‑250 (MQ=255)
        ttCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAacaccacac                                                                                                                                                                                                                                      <  1:394378/247‑1 (MQ=255)
         acGCCAGAAAACCGCGCCCGCCCTGTCTTTCCAGCGTCGCGCCAATATTCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACGTCAGTGACGAACTTCAGTATCTCAGCCAGCGTCGTTCTGGAAAGCTCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCTCACACGCCTGTTGGGTAAACACCTCACCATCg                                                                                                                                                                                                                                 <  2:34703/250‑1 (MQ=255)
           gCCAGCAAACCGCGCTCTCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAAGCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTg                                                                                                                                                                                                                               <  1:249465/251‑1 (MQ=255)
           gCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATACCGCACACTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCCGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTg                                                                                                                                                                                                                               >  1:516121/1‑251 (MQ=255)
            ccAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGc                                                                                                                                                                                                                              >  1:418095/1‑251 (MQ=255)
                cAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTCCCGCTCACACCTTCCTGCGATAACGCACAATTTAGCTTCCCGGAATCGAGAATAAAGCGGATCGCCGCAACCGGGACTTCATTGACGATCTTCAGGATCTCAGACAGCGTCGGTCTTGAAAGCACCATCAGCGGAAGCTCTTACTATACCTCTACCACCACCGCCTGTTGGGTAAacaccacacactcgatcggtc                                                                                                                                                                                                                          >  2:4885/1‑239 (MQ=255)
                    ccGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGat                                                                                                                                                                                                                      <  1:450816/251‑1 (MQ=255)
                    cagcgcTCGCACTGTTTAACCAGAGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGCATCGAGAATAAAGCGGATCGACGCAACCGGGACTTCATCGACGCACTTCAGGATCTCCGCCAGCGTCGTTCCGGAAAGAACCGTCAGCGGAAACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGAGTCCACCACACCCCATCGTGCGTCTCGat                                                                                                                                                                                                                      >  2:550372/3‑251 (MQ=255)
                           cGCACTGTTTTTCCAGCGTCGCGCCGATATGCAGCCCCCACTTCCCGCTCAACCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAACGCGGATCGCCGCAAACGGGCCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGcac                                                                                                                                                                                                                >  1:193640/1‑250 (MQ=255)
                           cGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGcac                                                                                                                                                                                                                <  1:373377/250‑1 (MQ=255)
                           cGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGcac                                                                                                                                                                                                                >  1:318821/1‑250 (MQ=255)
                                       ccAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATg                                                                                                                                                                                                    <  2:131197/250‑1 (MQ=255)
                                         aGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACTCCAGCGTGCGACTCGATATGCACAATGTTGGTATGccc                                                                                                                                                                                                 <  1:441597/251‑1 (MQ=255)
                                           cGTCGCGCCAATATGCAGCCCCAACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGCAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGACCTTCAGGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGCCTCTTGCTCGCCATCTGCCACACACCCCTGTTGGGGCAACACCCCACCCTCGTGCGCCTCGTTATGCCTAATTGTGGTATGCccgc                                                                                                                                                                                               >  2:267780/1‑251 (MQ=255)
                                            gTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgcc                                                                                                                                                                                              >  1:204069/1‑251 (MQ=255)
                                                cgcCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCAGCGATAACGCACAATTTAGCGTCGCGGAATCGAGAATAAAGCGGAGCGCCGCAAACAGGACTTCTTTGACACATTTCAGGACCTCAGCCAGCGTCGTTATGGAAAGCACCGTCAGCGGAGACTCTTGCTCTCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGGCTCGATATGCACAATGTTGGTATGCGCGCcgacg                                                                                                                                                                                          >  2:143065/1‑251 (MQ=255)
                                                      tatGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGACTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTAGGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAACCACCACACCATCGTGCGTCTCGCTATGCACCATGTTGGTATGCCCGCCGACGATggcg                                                                                                                                                                                    >  2:8036/1‑249 (MQ=255)
                                                         gcaccccccACTTCCCGCTCAAACCTCCCGGCGATAACGCACAATTTAGCTGCGCGGAATCGAGAATAAAGCGGAGCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGaca                                                                                                                                                                                 <  2:482580/247‑1 (MQ=255)
                                                                 cACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTACCACACGccc                                                                                                                                                                          >  1:316673/1‑250 (MQ=255)
                                                                       ccGCTCAAACCTTCCTGCGATAACGCACAATTCAGCTTCGCGGAAGCGAGAATAAAGCGGATGGCCGCAAACGGGACTTCATTGACGAATTTCAGGAACTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCGCTGCCACACACGCCTGTTGGGTAAACTCCACACCATCCTGCGTCTCGATATGCACAATGTTGGTATGCCCGGCGACGATGGTGACACGCGCCCACTtctc                                                                                                                                                                   <  1:415357/251‑1 (MQ=255)
                                                                             aaaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGtt                                                                                                                                                             >  1:420189/1‑251 (MQ=255)
                                                                                ccTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCa                                                                                                                                                          <  1:326350/251‑1 (MQ=255)
                                                                                  ttCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTcc                                                                                                                                                           >  1:260382/1‑248 (MQ=255)
                                                                                        cGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGAGACACACGCCCACTTCTCACCGTTCCAGACTTTgg                                                                                                                                                  >  1:471327/1‑251 (MQ=255)
                                                                                         gATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGACTCTTGCACGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGAGTCTCGCTCTGCACCATGTTGGTATGCCCGCCGACGCTGGTGCCACACGCCCACTTCTCACCGTTCCAGCCTTTGgc                                                                                                                                                 >  2:511353/1‑251 (MQ=255)
                                                                                              cGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCCTCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTACCACCGTTCCAGACTTTGGcgcggg                                                                                                                                            >  2:546911/1‑249 (MQ=255)
                                                                                              cGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                                            >  1:541154/1‑251 (MQ=255)
                                                                                              cGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                                            <  1:503888/251‑1 (MQ=255)
                                                                                              cGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                                            >  1:203865/1‑250 (MQ=255)
                                                                                               gCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGAATGCCCGCCGCCGATGGTGACACACGCCCATTTCTCACCGTTCCaacctctgcacggtcg                                                                                                                                           >  2:87115/1‑236 (MQ=255)
                                                                                               ccacaATTTATCTTCGCTGAATCGAGAATAAAGCGTGTAGCCACAAACGGTACTTCATTGACACACTCCATGATCTCCACCAGCGTGGTTCTTGAAAGCACCGTCACGTGAGACTTTTGCTCGCCCGCTGCCACTCACGCCTGTTGGGTAATCACCACTCCATCGTGCGTTTCGATATGCACAATGGTGGTAGTCCCGCCGACGATGGTGACACACGCCCACTTCTCTCCGTTCCAGACTTTGGCGCGTga                                                                                                                                           <  2:403227/250‑1 (MQ=255)
                                                                                                  caATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACCCACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGagaa                                                                                                                                        <  1:300529/251‑1 (MQ=255)
                                                                                                  caATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATAGACGAACTTCAGGCTATCAGCCACCGTCGTACTGGCAAGCACCGTAGCCGGAGACTCTTGCTAGCCCTCTGCCACACACGCCTGTTGGGTATACTCCACACCATCGGGCGTCTAGATATGCACAATGTTGGCATGACCGCCGACGATGGCACCACACAACCACTTCTCACCGTTCCAGACTTTGGagcgggcgc                                                                                                                                         >  1:68634/1‑245 (MQ=255)
                                                                                                        aGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGACCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCGGTTGGGTAAACACCACACCATCGTGCGTCGCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGAt                                                                                                                                  <  1:76061/251‑1 (MQ=255)
                                                                                                          cTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAACCACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGAttt                                                                                                                                >  2:410581/1‑251 (MQ=255)
                                                                                                            tCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGAttt                                                                                                                                <  1:571398/249‑1 (MQ=255)
                                                                                                                 gaagcGAGAATAAAGCGGATCGCCGCAAACGGGACTTCAATGACGAACTTCAGGATCTCAGCCAGCGTCGTTATGGAAAGCCCCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTGCTCACCGTTCCAGACGTTGGCGCGTGAGAAGAGGATTTCATCGCa                                                                                                                         <  2:98812/247‑1 (MQ=255)
                                                                                                                  aaTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGACACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCa                                                                                                                         >  1:228417/1‑250 (MQ=255)
                                                                                                                   aTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCa                                                                                                                         >  1:499790/1‑249 (MQ=255)
                                                                                                                      gagaAAAAAGCGGATCGCCGCAAACGGTACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCTCACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGtt                                                                                                                    <  1:502363/251‑1 (MQ=255)
                                                                                                                          aTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGACAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGAATGCCCGCCGACGCTGGAGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAAAGGATTTCATCGCAAGGTTCCTg                                                                                                                >  2:409506/1‑251 (MQ=255)
                                                                                                                           tAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGCCGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTgg                                                                                                               >  2:421777/1‑251 (MQ=255)
                                                                                                                                 ggATCGCCGCAAACGGGACTTCATGGACGAACTTCCGGACACCAGCCGGCTGCGTTCTGGAAAGCACCGGCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACCCCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTa                                                                                                         <  2:541154/251‑1 (MQ=255)
                                                                                                                                 ggATCGCCCCAAACGGTACTCCTCTGACGAACGTCAGTACCTCACCCAGCGTCGTTCTGGAAAGCACCTTCAGCGTTGACTCTTGCTCTCCCCTTGCCACACACGCCTGTGGTGTAAACACCACTCCATCGTGCGTCGCGATATGCACAATGGTTGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTGCCATACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTa                                                                                                         <  2:535733/251‑1 (MQ=255)
                                                                                                                                    ccgccacaaactgagcTTCCTTGCCGAACTCCAGGATCGCAGCCAGCTTCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTCCTCTCCCTCTGCCACACACGCCTGTTTGGAAAACACCTCCCCAGCGTGCGTCTCGATATGCACATTGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCCTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACg                                                                                                      <  2:161846/236‑1 (MQ=255)
                                                                                                                                          cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATCTGCACAATGTTGGTATGCCCGCCGACGATGGTGCCACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGGGAGAAGAGGATTTCATCGAAAGGTTCCTGGATCTTAACGGAGACt                                                                                                >  2:358714/1‑251 (MQ=255)
                                                                                                                                             aCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTc                                                                                             <  1:488928/251‑1 (MQ=255)
                                                                                                                                                   cTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgcc                                                                                       >  1:413682/1‑251 (MQ=255)
                                                                                                                                                     tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTgaga                                                                                                                                         <  1:16003/199‑1 (MQ=255)
                                                                                                                                                     tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTgaga                                                                                                                                         >  2:16003/1‑199 (MQ=255)
                                                                                                                                                       aTTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGCTGGTGACACACGCCCACTTCTCACCGTTCCCGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATATTAACGGATACTTTCCCCgccgcc                                                                                    >  2:559501/1‑250 (MQ=255)
                                                                                                                                                       aTTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCACgccgcc                                                                                    >  1:152047/1‑250 (MQ=255)
                                                                                                                                                        ttGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCgccgcc                                                                                    >  1:196103/1‑249 (MQ=255)
                                                                                                                                                         tGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCcagc                                                                                 >  1:156534/1‑251 (MQ=255)
                                                                                                                                                         ggccGAACTTCAGGATCCCAGCCAGCTTCGTTCTGGAAAGCACCTCCAGCGGAGACTCTTGCTCTCCCTCTGCCACACACGCCTGTTGGTTAAACCCCCCACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcag                                                                                  <  1:8036/247‑1 (MQ=255)
                                                                                                                                                          gACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCCAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcagc                                                                                 >  1:322136/1‑250 (MQ=255)
                                                                                                                                                                 gtCAGGATCTCAGCCACCGTGCGTCGGGAACGCCCCGTCAGCGGAGACCCTGGCCCGACCGCTGCCACACACGCCTGTTGGGTAAACCCCACACCATCGTGCTCCACGATATGCACAATGTGGGTATGCCCGCCGACGATGGTGACACCCGCCCTCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCCAGCAGTGCttt                                                                         <  2:155484/250‑1 (MQ=255)
                                                                                                                                                                  tCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGACCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCCGTGCTTTg                                                                        >  1:537976/1‑251 (MQ=255)
                                                                                                                                                                       aTCTCAGACAGCGTCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTTTTGGGTAAACACCACACCATCGTGCGTCTCGATATTCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGACGAGGCTTTCCTCGCCAGGTTCCGGGATCTTAACGGAGACTTTCCCCGCCGCCAGCCGGGCTTTGGcctc                                                                   >  2:339952/1‑248 (MQ=255)
                                                                                                                                                                            agccagcGTCGTTCTGGAAAGCACCGTAAGCGGAGACTCTTGCTAGCCCTCTGCCACACACGCCTGTTGGGTAAACACCCCACCATCGTGCGTCTCGATATGCACAAGGTTGCTATGCCCGCCGCCGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAACCTTTCCCCGCCGCCAGCAGTACTTTGGCATCGGCa                                                               >  1:508065/1‑250 (MQ=255)
                                                                                                                                                                               cagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAAtt                                                            >  1:96265/1‑250 (MQ=255)
                                                                                                                                                                                agcGTCGTTCTGGAAAGCACAGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGCTGCCAAACGCACACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAAGTTCCTGGATGTTAACGGAGACATAACCCGCCGCCAGCCGTGCTTTGGcgccggccatttc                                                          >  2:52905/1‑239 (MQ=255)
                                                                                                                                                                                 gcGTCGTTCTGAAAAGCCCCGTCAGCGGAGCCTTTTGCTCGCCCCCGGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGcc                                                         <  1:562273/251‑1 (MQ=255)
                                                                                                                                                                                           gggAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGGGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGtt                                                <  1:498324/250‑1 (MQ=255)
                                                                                                                                                                                            ggAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGac                                                                                                                                                                                  >  1:74785/1‑120 (MQ=255)
                                                                                                                                                                                                      gTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTACACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCACGGTTCCTGGATCTTAACGGCGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTCCGCTGTTGCGTCttta                                       >  2:61706/1‑248 (MQ=255)
                                                                                                                                                                                                      gTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCac                                     >  2:201710/1‑250 (MQ=255)
                                                                                                                                                                                                         aGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCCAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTTCCTGCCCTGTTGCGTCTTTcaccac                                  >  2:267229/1‑251 (MQ=255)
                                                                                                                                                                                                          gCGGAGACCCTTGCCCGCCCTCTGCCACACCCGCCGGTTGGGAAAACACCACACCATCGTGCGTCTCGATAGGCACAAGGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACt                                 <  1:339952/251‑1 (MQ=255)
                                                                                                                                                                                                             ggAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTc                               >  1:489555/1‑251 (MQ=255)
                                                                                                                                                                                                                    tctTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTcccc                                                                                          >  2:422396/1‑186 (MQ=255)
                                                                                                                                                                                                                    tctTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTcccc                                                                                          <  1:422396/186‑1 (MQ=255)
                                                                                                                                                                                                                    tctTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc                         >  1:371586/1‑251 (MQ=255)
                                                                                                                                                                                                                             cccTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTATTTCAGCACTTCCAGCCCGGCGTTGGc                >  2:393447/1‑251 (MQ=255)
                                                                                                                                                                                                                                tctGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgccg                                                                                      >  1:296415/1‑178 (MQ=255)
                                                                                                                                                                                                                                 cggCCACACACCCCTGTTGGGTAAACACCACACCATCGTGCGTCGCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTGTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAt              <  1:358714/247‑1 (MQ=255)
                                                                                                                                                                                                                                       cacccGCCTGTTGGGTAAACACCACCCCATCGTGCGTCCCTATATGCACAAGGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGGGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGATTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCt      <  1:107656/247‑1 (MQ=255)
                                                                                                                                                                                                                                          acGCCTGTTGGGTAAACACCACGCCATCGTGCGTTCCGATATGCACAATGTGGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCGCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAAc   <  1:93054/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCTCAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACg  >  1:144217/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
AGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACG  >  NC_000913/3255919‑3256401

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