Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,535 A→G D177D (GAT→GAC yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,5350AG98.7% 219.7 / ‑4.7 75D177D (GAT→GACyjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (35/38);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

CGGCAATGGCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCG  >  NC_000913/4616291‑4616775
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
cggtaaTGGCTTTCTGGTTAGTGCGTACCAGCGTTCTCCCACCGCCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAACTTTATCAGGCACCGGACCGGTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTGCGGGGAAACCTGCTTCAACtcggcgtctg                                                                                                                                                                                                                                               <  2:656823/245‑1 (MQ=255)
cGGCAATGGCTTACAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACtctgcgtctgc                                                                                                                                                                                                                                              >  1:173264/1‑250 (MQ=255)
 ggCAATGGCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCAc                                                                                                                                                                                                                                            >  1:452209/1‑251 (MQ=255)
  gCAATGGCTTTCAGTATAAGGCTTACCAGCGTGCTGCCTCCCCCCGCCAGCGGCAGTAAGTGTTAACAGCTCGCAATCCCGTAGTCCCCTTTTGTGAAAATTTTATCATGCTCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGAGGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGAGCGTAGATAAAGGGCAGGTTCGTGGAAACCTGCTTCAACTCGGCGTCTGCACg                                                                                                                                                                                                                                           <  2:64576/251‑1 (MQ=255)
  gCAATGGCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACg                                                                                                                                                                                                                                           >  1:17418/1‑251 (MQ=255)
        gCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGAt                                                                                                                                                                                                                                     <  1:591484/251‑1 (MQ=255)
           ttCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCAAGTAAGCTcgc                                                                                                                                                                                                                                  >  2:565428/1‑251 (MQ=255)
               ggTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCgg                                                                                                                                                                                                                              >  1:224137/1‑251 (MQ=255)
                 ttAAGGCGTACCAGCGTGCTGCCACCCCCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACtctgcgtctg                                                                                                                                                                                                                                               <  1:562713/232‑1 (MQ=255)
                 ttAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACtctgcgtctg                                                                                                                                                                                                                                               >  2:562713/1‑232 (MQ=255)
                    aGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGATCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCACGCGGGTAAt                                                                                                                                                                                                                         >  2:5178/1‑251 (MQ=255)
                     ggCGTACCAGCGGGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATg                                                                                                                                                                                                                        <  1:420296/251‑1 (MQ=255)
                         tACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCg                                                                                                                                                                                                                    >  1:21732/1‑251 (MQ=255)
                              gCGTGCTGCCCCCCCCCGCCAGCGGCACTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCTATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCa                                                                                                                                                                                                               <  2:460611/251‑1 (MQ=255)
                                gTGCTGGCACCACCCGCCAGCGGCAGAGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCACCTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAAATAAAGGTCAGGTTCGGGGAAACCTGCTGCAACTCTGCGTCTGCAGGTAAGATCGCGCGGGTCAAGGGCAAATc                                                                                                                                                                                                                >  2:438375/1‑248 (MQ=255)
                                gTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTAATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc                                                                                                                                                                                                                >  1:496679/1‑248 (MQ=255)
                                gTGCTGCCACCACACGCCAGCGGCAGTGAGTTGTAACCGCTCACAATACCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc                                                                                                                                                                                                                >  1:581171/1‑248 (MQ=255)
                                  gccgcccccccccGCCAGCGGCAGTGAGTTGTAACAGCTCACAACTCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGAc                                                                                                                                                                                                            <  2:428348/240‑1 (MQ=255)
                                     gcCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACggg                                                                                                                                                                                                         <  1:568712/250‑1 (MQ=255)
                                         ccaccCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGat                                                                                                                                                                                                    >  1:394710/1‑251 (MQ=255)
                                              cGCCATCGGCAGTGATTTGTTACAGCTCCCAATCCCGTAGCCCCCTTTGGTGAAAATTTTATCATGCACCGTACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTCATGGGCGAATCTGCCGTGCCGTTATTgg                                                                                                                                                                                               <  2:224590/251‑1 (MQ=255)
                                                ccAGCGGAAGCGAGTTGTAACAGCTCCCAAACCCGTAGACCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGAGATGCGCGGTTTGCTACATTAACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGCTTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCACCTCTGCGTCTGAACGTAAGATCGAGCGGGTAATGGGCGAATCAGACGGGCCGATACTGGCg                                                                                                                                                                                             >  2:286556/1‑251 (MQ=255)
                                                       cagggAGTTGTAACAGCTCAAAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGGGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGCTTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGATTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAAGGGGCAACTCAGACGGGACGGTATTTGCGTGAATa                                                                                                                                                                                       >  2:660502/5‑250 (MQ=255)
                                                       cAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATaa                                                                                                                                                                                      <  1:441110/251‑1 (MQ=255)
                                                       cAGTGAGCTGTAACAGCACACAATACCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACAGTTGGCGATGAGCGGTTTGCTACATTCACAGATGTTCTTCGACACTTCCAGAAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTACTAGGAGAATCAGACCGGCCGATATTGGAGTGCATa                                                                                                                                                                                       >  2:138898/1‑250 (MQ=255)
                                                           gAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGc                                                                                                                                                                                  >  1:193642/1‑251 (MQ=255)
                                                            ggaggTAACAGCTCAAAACCCCGTGGCCCCCGTTTGGGAAAATTTTATCATGCCCCGGACCGTTGGCGATGTGCTGGTTGCTACATTCACGGATGTGCTTCGCCACTTCCAGCAGCAGGGCATCAGGGGTGATTGCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGGCTGCACGTAAGATCGCGCGGGTAAGGGGCGAAGCAGACGGGCCGATATTGGCGTGCATAAAGGCg                                                                                                                                                                                 <  2:301759/247‑1 (MQ=255)
                                                                        tCACAAACCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGGTGTGCGGTTTGCTACATTCACAGATGTTCTTCACCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAATCTGCTTCAACTCTGCGTCTGAACGTAATATCGCGAGGGTACTGGGAGAATCAGTCGCGCAGATAGTGGCGTGCATAAAGGCGTCTGGCAgggat                                                                                                                                                                     >  2:262755/1‑249 (MQ=255)
                                                                               cccGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCgg                                                                                                                                                              >  1:145416/1‑250 (MQ=255)
                                                                                 cGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAgg                                                                                                                                                            <  1:461265/251‑1 (MQ=255)
                                                                                  gTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGt                                                                                                                                                           >  1:256220/1‑251 (MQ=255)
                                                                                    tgCCCCCTTTTGTGAAAATTTTATCATGCTCCGGACCGTTGGCGATGTGCGTTTTGCTACATTCACAGATGTTCTTCGCCTCTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTa                                                                                                                                                          <  1:655079/249‑1 (MQ=255)
                                                                                    aGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGAGTCTGGCAGGGGTCTATCGAGGTa                                                                                                                                                          >  2:563451/1‑250 (MQ=255)
                                                                                              tgtgAAAATTTTATCATGCACAGGACCGTGGGCGATGTGCGGTTTGTTACAGTTACAGATTTTCTTCGCCACTTCCAGCAGCAGGTCAGAAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCGGCGTCTGCACGTAAGAGCGCGCGGGTAAGGGGCGAATCAGTCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGaaa                                                                                                                                               <  2:98312/250‑1 (MQ=255)
                                                                                                tgAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGt                                                                                                                                            >  1:617212/1‑251 (MQ=255)
                                                                                                  aaaaTTTAATCATGCACCGGACCGTTGGCGATGTGCGTTTTGCAACATTCACAGATGTTCTTCGCCACTCCAAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAAGTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGttt                                                                                                                                          <  2:262278/251‑1 (MQ=255)
                                                                                                         taccaggCACCGGTCCGGTGGCGAGGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCg                                                                                                                                    <  2:8783/244‑1 (MQ=255)
                                                                                                              tgctccGGACCGTTGGCGATGTGCGGTTTGCTCCATGCACAGAGGTTCTGGGCCACTTCCAGCTGCAGGTCATCAGGGGTGATTGCAGGATCGTAGATAAAGGGCAGGTGCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGaagat                                                                                                                               <  2:599238/246‑4 (MQ=255)
                                                                                                               gCACCGGATCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATc                                                                                                                              <  1:232227/250‑1 (MQ=255)
                                                                                                                 aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGAt                                                                                                                           >  1:339681/1‑251 (MQ=255)
                                                                                                                      aCCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCCGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGACGGGTTCTGTCGGGGTAACGCAAGAAACGTTTTATTCGAATATCGctctttt                                                                                                                      >  2:421853/1‑248 (MQ=255)
                                                                                                                        cGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGaaa                                                                                                                                               >  1:631979/1‑224 (MQ=255)
                                                                                                                            ggCGATGTGCGGTTTGCTCCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGGGATTTCAGGATCGTAGATAAAGTTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTtgtg                                                                                                                <  2:21732/251‑1 (MQ=255)
                                                                                                                                       gTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCt                                                                                                      >  1:65698/1‑250 (MQ=255)
                                                                                                                                          tGCTAAATTCACAGCGGTTCTTGGCCACTGCCAGACGCAGTTAATCAGGGGTGATTACAGGATCTGAGATAAAGGTCAGGTTAGGGGAAACCTGCTTCAACTCTGCGTTTGCACGGAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATGGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCGTGTGTTAGAATTCTAAc                                                                                                   <  2:102972/250‑1 (MQ=255)
                                                                                                                                              aCATTCACATATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGAGCGTAGAGAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAAGGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCATAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATa                                                                                                <  2:256220/249‑1 (MQ=255)
                                                                                                                                               cagtCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGCCAGGTTCGGGGAATCCTGCTTCAACTCTGCGTCTGCACTTAATATCGCGCGGGTAATGGGCGAATCAGACTGGCCGATATTGGCGTGCATAAAGGCGTCTTGCAGGGTTCTGTCGAGGTAACTCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATatgg                                                                                             <  1:326775/248‑4 (MQ=255)
                                                                                                                                                 ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGCTTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAACCGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGtt                                                                                           >  2:267569/1‑251 (MQ=255)
                                                                                                                                                  tCACAGATGTACTTAGCCACTTCCAGCAGCAGGTCATCAGGGTCGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGCAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGCGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGCTATGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGCTCTCGTCTTGTGGTAGAATTCTAACATATGGTTg                                                                                          >  2:92063/1‑251 (MQ=255)
                                                                                                                                                    tcaGATGTTCTTAGACAATTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACATGCTTCAACTATGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGACCGATATTGGCGTGAATAAAGGCGTATGGCCGGGGTCTGCCGAGGAACAGCACGAAACGTTTTATTCGAAAATCGATCTCATCATGTGTTAGAATTATAAAATagtgttgc                                                                                         >  2:209478/2‑243 (MQ=255)
                                                                                                                                                      agttGTTCTTCGCAATTGCCAGCAGCAGGTCATCAGGGTTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCGGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGTCAGGGTTCTGTCGAGGTAACGCCAGAAACTTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGcaac                                                                                      <  2:259790/248‑1 (MQ=255)
                                                                                                                                                        aTGTTCTGCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGTTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGcaacaa                                                                                    <  1:421853/251‑1 (MQ=255)
                                                                                                                                                              ttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTaa                                                                                                                                                         >  2:326518/1‑177 (MQ=255)
                                                                                                                                                              ttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTaa                                                                                                                                                         <  1:326518/177‑1 (MQ=255)
                                                                                                                                                               tcGCCCCTCCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCAt                                                                               <  1:526363/249‑1 (MQ=255)
                                                                                                                                                                     cTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACATGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGACTCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATATCGTCTTGTGTTAGAATTCTAACCTATGGTTGCAACAACGCATCCAGTTGc                                                                       >  2:644094/1‑251 (MQ=255)
                                                                                                                                                                       tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGAGCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGcc                                                                      <  1:644094/250‑1 (MQ=255)
                                                                                                                                                                        ccAGCATCAGGTGATCAGGGGGGATTTCAGGATCGTAGATAAAGGTCATGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGGCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGccc                                                                     <  1:295162/250‑1 (MQ=255)
                                                                                                                                                                            cagcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                                >  1:191345/1‑251 (MQ=255)
                                                                                                                                                                                 ggTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACg                                                                                                                                                      <  1:555117/160‑1 (MQ=255)
                                                                                                                                                                                 ggTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTAAGGTTCGGGGCAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAGAGGCGTCTGGCAGGGTTCTGTCGAGGGAACg                                                                                                                                                      >  2:555117/1‑160 (MQ=255)
                                                                                                                                                                                    caccaGGGATGATTGCAGGCTCGTAGATAAAGGGCAGGTGCGTGGAAATCTGCTTCAATACCGCGTCTGCACGTAAGCTCGCGCGGGTAATGGCCGAATCAGACGGGCCGATTTTGGCGTGCATAAAGGCGTCTGGCATGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATGCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                                        <  2:193642/248‑1 (MQ=255)
                                                                                                                                                                                       caGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAATCCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGGCTGGCAGGGTGCTGTCGAGGTAACGCCAGAAGCGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAAGTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                                                     <  1:286556/251‑1 (MQ=255)
                                                                                                                                                                                          gggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                                                  <  1:132524/251‑1 (MQ=255)
                                                                                                                                                                                                 ttctggaccGTAGATAAAGGTCAGGTTCGGGGAAACCGGCTGCAGCTCTGCGGCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc                                            <  2:191345/242‑1 (MQ=255)
                                                                                                                                                                                                 ttCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc                                            <  1:598407/250‑1 (MQ=255)
                                                                                                                                                                                                   cAGGATCGTAGATAACGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATCTTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGGAGACCGGCAtgtatatgacctac                                         >  2:41857/1‑238 (MQ=255)
                                                                                                                                                                                                   cAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCTTCTGCGCGATAGAGCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc                                         <  2:145416/251‑1 (MQ=255)
                                                                                                                                                                                                   cAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGATCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATTTGACCGAc                                         >  2:42081/1‑251 (MQ=255)
                                                                                                                                                                                                     ggATCGTAGATAAAGGTCAGGTGCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACgg                                       <  1:228447/251‑1 (MQ=255)
                                                                                                                                                                                                      gATCGTACATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGGTCTGTCGAGGTAACGACAGAAACGGTTTATTAGACTATCGATCTCGCTTTGTGTTAGAATTCTAACGTATGTATGCCACAAGGCATCCCATTGCCCCAGGTAGACCGGCGACGATGGTACCGACGGt                                      >  2:268538/1‑251 (MQ=255)
                                                                                                                                                                                                         cGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc                                    <  1:565428/250‑1 (MQ=255)
                                                                                                                                                                                                                  aGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGTCCGGCATCGATGTGACCGACGGTACGTGGTGGTaa                          <  1:520898/251‑1 (MQ=255)
                                                                                                                                                                                                                  aGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaa                          <  1:267569/251‑1 (MQ=255)
                                                                                                                                                                                                                           tCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGAGCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCATAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGGCGGTACGTGGTGGTAAAGAATGGt                  <  2:339681/250‑1 (MQ=255)
                                                                                                                                                                                                                             ggggaacccttcTTCAACTCTGCGTCTGCACGTAACATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcag               <  1:544939/240‑1 (MQ=255)
                                                                                                                                                                                                                                      tGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCagagag        >  1:209633/1‑249 (MQ=255)
                                                                                                                                                                                                                                        cTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCCTCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTg    >  1:303039/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGACTCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGAAGACCGGCATCGATGTGACCGAAGGTACGTGGTGGTAAAGAATGGGCAGAagagaggttgcg  >  2:511526/1‑245 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
CGGCAATGGCTTTCAGGTTAAGGCGTACCAGCGTGCTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCG  >  NC_000913/4616291‑4616775

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