Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,571,280 C→T E99K (GAG→AAG)  yjiR ← fused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,2800CT98.8% 281.8 / NA 82E99K (GAG→AAG) yjiRfused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/54);  minor base G (0/1);  total (26/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTC  >  NC_000913/4571041‑4571505
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
aTGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGAACACCGAGCCACTGTCGAGCATCAGGCGGGGAATCTGCTCTCGCAATACTCGCTGACCAGCTAACCAGGCATAACCGCGAACGGGTTGCAGATTTTTCTGCACCAACCGGCTTAGCTTCCGCCACAGCGGGTTTAGGCTGGGCGCTTCGACATCCGGGGTGCTTTTGCTTAACGGAACAATGGAACAGTCGCTCTGCGGCTTTAGCATATccc                                                                                                                                                                                                                            >  1:734192/1‑250 (MQ=255)
  ggcggCCGCGGGTGATGATGATGTCAGCGGCGGTGCCCACCGAGCCGCTGTCGAGCAGCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCAAAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCGTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAg                                                                                                                                                                                                                           <  1:692047/246‑1 (MQ=255)
  ggCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAAACTGCTCTCGCAATAATCGCTGGACGGCTAACAAGGCATAAACGAGAACGGGTTGCAGATTATGCTGCACCACCCCGCTTAGCTCCCGCCACAACGGGTTCAGGCTGGGCGCTTCGACATCCGGTGTGGTTTTGCTTAACGGAACAATGGGACTGTCGCTTTGCGCCTTCAGCATATcccg                                                                                                                                                                                                                           >  1:580186/1‑248 (MQ=255)
     aGCCGCGGGTGAGGATGATGTCACCGGCGGTGACCACCGAGCGGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCAcc                                                                                                                                                                                                                       <  1:502824/251‑1 (MQ=255)
              ggatgatgatGCTATCGGCGGTGCCCCCCGAGCCTCTGTCGCGCAACAGGCGGGCAATCTGCCCGCGCAATACTCCCTGACCGGCTAACAGGTCATAGCCGAGACCGGTTTGCAGATTATGCGGCCCCACCCGGCTTAGCTCCCGCCACAGCGGTTTCATGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTACAGCATATCCAGCACCTGATCCCAc                                                                                                                                                                                                              <  2:522439/250‑1 (MQ=255)
                          cATCGGCGGTGCCCACCGAGCCGCTTTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTc                                                                                                                                                                                                   <  1:561474/250‑1 (MQ=255)
                            tCGGCGGTGAACACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGCTTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACCTCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCCGCACCTGATCCCACTGGGTAATTTCCAc                                                                                                                                                                                                >  2:485338/1‑251 (MQ=255)
                            gcggcggTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCAc                                                                                                                                                                                                <  1:207246/250‑1 (MQ=255)
                               gcggTGCCCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCgg                                                                                                                                                                                             <  1:574421/251‑1 (MQ=255)
                               gcggTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCgg                                                                                                                                                                                             >  1:597957/1‑251 (MQ=255)
                                cggTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCggg                                                                                                                                                                                            >  1:33259/1‑251 (MQ=255)
                                  gTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGgcg                                                                                                                                                                                          >  1:514393/1‑251 (MQ=255)
                                  gTGACCACCGAGCCGCTGTCGAGCATCAGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGgcgc                                                                                                                                                                                         >  2:108795/1‑251 (MQ=255)
                                    gACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCt                                                                                                                                                                                        <  1:98483/251‑1 (MQ=255)
                                        cccGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCAc                                                                                                                                                                                    <  1:198335/250‑1 (MQ=255)
                                            ggCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTGGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGTCGCTGCACCGGa                                                                                                                                                                                <  1:221865/250‑1 (MQ=255)
                                             gCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCCCCGGAc                                                                                                                                                                               >  2:184984/1‑251 (MQ=255)
                                             gCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAGGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGAc                                                                                                                                                                               <  1:47728/251‑1 (MQ=255)
                                               cGCTGTCGAGCATCCGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCCTATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGAcgc                                                                                                                                                                             >  1:106221/1‑251 (MQ=255)
                                                gCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGAcgcg                                                                                                                                                                            >  1:162418/1‑251 (MQ=255)
                                                 cTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGt                                                                                                                                                                           >  1:360943/1‑251 (MQ=255)
                                                  tGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGATATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTc                                                                                                                                                                          >  1:168721/1‑251 (MQ=255)
                                                      gAGCATCAGGCGGGCAATCTCATCTCGACCTACTCGCTGTCCGGCTAACAAGTCATACCCGAGAAAGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCTGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGTCGCTGCACCGGACGCGTCATcg                                                                                                                                                                      >  2:735190/1‑251 (MQ=255)
                                                        gCATCAGGCGGTCAATCGTCTCTCGCAATACCCGCTGCCCGGCTAAAAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcggc                                                                                                                                                                    <  2:551600/251‑1 (MQ=255)
                                                        gCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGATTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCCCCGGACGCGTCATcggc                                                                                                                                                                    >  1:327995/1‑251 (MQ=255)
                                                         cATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcggcg                                                                                                                                                                   <  1:252103/251‑1 (MQ=255)
                                                           aaaGGCGGTCAATCTGCTCTGCCAACCCTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATGAGGCCGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGt                                                                                                                                                                 <  2:584419/249‑1 (MQ=255)
                                                               gCGGGCAATCGTCTCTCGCAATCCTCGCTGACCGGCAAACAAGTCATACCCGAAAACGGTTTGCAGATTATGCTGCCCCACCCGGCTTAGCGCCCGCCACAGCGGTTTCAGGCGGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTg                                                                                                                                                            <  2:444190/251‑1 (MQ=255)
                                                                  ggCAATATGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAGGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTggcg                                                                                                                                                         <  1:416799/251‑1 (MQ=255)
                                                                   gCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTggacg                                                                                                                                                        >  1:149840/1‑248 (MQ=255)
                                                                         tGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTaacggcgggggggg                                                                                                                                                   >  2:195838/1‑238 (MQ=255)
                                                                             cccGCAATACTCGCTCCCCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTa                                                                                                                                             <  1:92306/249‑1 (MQ=255)
                                                                             cacgcacaccTCGCTGACCGGTAAACAAGTAATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCttt                                                                                                                                              <  2:530845/241‑1 (MQ=255)
                                                                                  aaTACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTtg                                                                                                                                       <  1:73309/251‑1 (MQ=255)
                                                                                      ccggctgccgggtcaaCAAGTCAAACCCGGAAACGGTTGGCGGTTTATGCGCCACCACCGGGCTTGGCCCCCGCCACAGCGGTTTCAGGCGGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCg                                                                                                                                   <  2:97956/237‑1 (MQ=255)
                                                                                           tGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTACTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGCCGCGTCATCGGCGGTACCGGCGGCCgggggtatcattttgcggcaa                                                                                                                                 >  1:640330/1‑232 (MQ=255)
                                                                                              ccGGCAAACAAGTCTAAACCGAGAACGGTTGGCGGATTATGCGGCCCCACCCGGTTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATa                                                                                                                           <  2:602294/251‑1 (MQ=255)
                                                                                                   tAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGc                                                                                                                                                                                                                                   >  1:175004/1‑145 (MQ=255)
                                                                                                   tAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGc                                                                                                                                                                                                                                   <  2:175004/145‑1 (MQ=255)
                                                                                                     aaaaGTCATAACCGAGACCGTTTTCAAGATGTTGCTGCCCCACCCGGCTGAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGaa                                                                                                                   <  2:594607/249‑1 (MQ=255)
                                                                                                     aCAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAACTGCCGGTGGGGTTTTCGTtgttggccaaaaaaaaaccg                                                                                                                     >  1:108880/1‑233 (MQ=255)
                                                                                                      aaaGTCATACCCGAGAACGGTTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAc                                                                                                                  <  2:557950/250‑1 (MQ=255)
                                                                                                       aaGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGCTTACGTTTGTGGGCCAAAAAAAcccaacc                                                                                                                  >  1:610099/1‑247 (MQ=255)
                                                                                                       aaGTCATAACCGGGACGGTTTGGCAGATTATGCTGCCCCCCCCGGTTTGGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAcg                                                                                                                 <  2:508689/251‑1 (MQ=255)
                                                                                                               cccGAAAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCg                                                                                                         <  2:153868/250‑1 (MQ=255)
                                                                                                                     accGGTTTGCATATTATGCGGCACCGCCCGGCTTAGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtga                                                                                                   <  2:618558/249‑1 (MQ=255)
                                                                                                                          gttggtagtttaggttgccccccccGGCTTAGCCCCCGCCGCAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTACTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGCTCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGGCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAg                                                                                               <  2:640330/229‑1 (MQ=255)
                                                                                                                                  aTTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGCTTACGTTTTGCGGCAAAATAACCCCCACGCGGCCGCGGGGGGGTGAGCTTTATcc                                                                                       >  1:635474/1‑250 (MQ=255)
                                                                                                                                     aTGCTGCACCACCCGGCTTAGCTCCCCCCACAGCGGCTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGGGGGTTTTCCTTGTGCGGCAAAAAAACCCCAACGCGGCGGGGGGGGGGTGAGGTTCATCGtct                                                                                    >  1:514220/1‑251 (MQ=255)
                                                                                                                                              caccCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCt                                                                           <  2:162418/251‑1 (MQ=255)
                                                                                                                                                cccGTCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGa                                                                         <  2:2220/251‑1 (MQ=255)
                                                                                                                                                   ggCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATAc                                                                      <  1:246742/251‑1 (MQ=255)
                                                                                                                                                    ccttcgCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCTAGCGTCTGATACg                                                                     <  1:566257/246‑1 (MQ=255)
                                                                                                                                                        aGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctg                                                                 <  1:515454/251‑1 (MQ=255)
                                                                                                                                                            cccGCCACAGGGGTGTCAGGGTGGGCGCAGCGACATCCGGGGTGCTTGTGCGTAACGGAACAATGGAACTGGCGCTATGCGCCTGCAGCAGAGCCAGCACCGGATCCGACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCGTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGAGGAGCTTCAGCGGCTCCAGCGTCTGATACGCctgctt                                                              <  2:266830/250‑2 (MQ=255)
                                                                                                                                                               gCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGGGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACg                                                          <  1:259880/251‑1 (MQ=255)
                                                                                                                                                                             aGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGAc                                            <  2:33259/251‑1 (MQ=255)
                                                                                                                                                                             aGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGAc                                            <  1:158657/251‑1 (MQ=255)
                                                                                                                                                                                   ggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTg                                      <  1:549402/251‑1 (MQ=255)
                                                                                                                                                                                    gcgcTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGc                                     <  1:167928/251‑1 (MQ=255)
                                                                                                                                                                                      gcTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGGCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACGGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTTTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCGGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTc                                   <  2:64256/251‑1 (MQ=255)
                                                                                                                                                                                      ccTTCGACAGCCGGTGTGCTGTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGTGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATATTGACGCCGTGCTc                                   <  2:106221/250‑1 (MQ=255)
                                                                                                                                                                                         tCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAAGTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCGCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCt                                 <  2:734192/250‑1 (MQ=255)
                                                                                                                                                                                           gACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAATGGGGGGTGGGGCTTTCGTTGGGCGGCAAAAAAAACCCGAACCGGCTTCGGGGTGGTGAGGTTCATCGCCTCCCCCGGCCTATACGCCCGCTGCCCGGGGGTTATACTCACGCCGGGCTccctac                              >  1:192293/1‑247 (MQ=255)
                                                                                                                                                                                            aCATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGCGTTTCGTTTTGCGGCAAATAAACCCCACGCGGCCGCGGAGGGGTGGGCATCATCCTCTCCCCCGTCATTTACGCCCGCTGCCCGGGGCTGATAcagtcaccccgctcccctatt                            >  1:573053/1‑231 (MQ=255)
                                                                                                                                                                                                ccGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTaa                          <  1:646484/250‑1 (MQ=255)
                                                                                                                                                                                                    tgtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTCGTTTTGGGGCAAAATAACCCCCAAGCGGCCGCGGGGGGATGAAGTTCCTCGTCTTCAGCGGCTGGTAAGCCTGCTGCACCGTGCTGGTTCTGGCGCCGGGGTCCTTACTTTAGccgc                     >  1:276846/1‑248 (MQ=255)
                                                                                                                                                                                                     gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTTCGTTGTGGGACAAAATAACCCCAACGCGGCCGCGGGGGGATGGGCTTCCTCGTCTCCCGCGGCTGATACGCCCGCTGCAc                                                           >  1:106160/1‑212 (MQ=255)
                                                                                                                                                                                                     gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCAc                                                           <  2:106160/212‑1 (MQ=255)
                                                                                                                                                                                                     gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTgc                     <  1:75314/250‑1 (MQ=255)
                                                                                                                                                                                                     gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTGAAGCTgc                     <  1:195838/250‑1 (MQ=255)
                                                                                                                                                                                                       gCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGGTACTGACGCCGTGCTCCTGACTTAAGCTgcgc                   <  1:640066/250‑1 (MQ=255)
                                                                                                                                                                                                           ttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCg               <  1:305882/250‑1 (MQ=255)
                                                                                                                                                                                                           ttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATAGCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCg               <  2:189140/250‑1 (MQ=255)
                                                                                                                                                                                                             gCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACg            <  2:583201/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaCGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCaa        <  2:115886/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaCGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCa         <  1:615384/250‑1 (MQ=255)
                                                                                                                                                                                                                   cGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGTAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAtt      <  2:304657/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaCAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCCCCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCGGCGCACCGACGGCAATTtct   <  2:360943/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTtctc  <  2:406836/251‑1 (MQ=255)
                                                                                                                                                                                                                         aaTGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACt                             <  2:699247/222‑1 (MQ=255)
                                                                                                                                                                                                                            ggAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCg                                                                                       <  1:558669/161‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
ATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTC  >  NC_000913/4571041‑4571505

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