Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,379 T→C 8.1% E13G (GAG→GGG)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,3790TC 8.1% 373.6 / 14.5 120E13G (GAG→GGG) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (53/55);  new base C (8/4);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

Only 100 of 121 total aligned reads displayed.

CGCTCGCGCCAGCGTTTCCACCACATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGG  >  NC_000913/53105‑53652
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                 
cGCTCGCGCCAGCGTTTCCCCCACATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCt                                                                                                                                                                                                                                                                                   <  2:138008/275‑1 (MQ=255)
        ccAGCGTTTCCACCACATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGa                                                                                                                                                                                                                                                                           <  2:276239/275‑1 (MQ=255)
                                            cTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCCATCCCAGCGGGCCCCGCGGGGGGGATCACAACACGTTTGGTTTTAACCATTAGGTg                                                                                                                                                                                                                                        >  4:440464/1‑274 (MQ=255)
                                                                ggTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGgcg                                                                                                                                                                                                                   <  1:106703/275‑1 (MQ=255)
                                                                    cACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTgg                                                                                                                                                                                                               <  2:27303/275‑1 (MQ=255)
                                                                    cACGTAGCGCGACAGGAAGTACGGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGCTGGGGGAATAAGGGCCGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGTGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGTGAGTGATCACAACACGTTGGGTTTTAACCATTATTTGCTCAGGATTTTAACGTAGGCGCTgg                                                                                                                                                                                                               <  1:62563/275‑1 (MQ=255)
                                                                                     aGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCACCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGATAACTAAGTCCGGGCCAATCCCGGCGGGCCCGCCGGGAGTGAACACAACCCGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTgac                                                                                                                                                                                                >  2:18471/1‑271 (MQ=255)
                                                                                                                 gctgttgcggcgggTTGGGGGAATCAGGGTGCTGGTTGAACGGCAAACCCGGGATCGTGTCCCGGTTGGGACGGGGAGTGGCATCGGCCCAAACAACCAGTCCGACCGGCCCCCAACGCTGTGCAAGCTGGACAACTAAGTCCTGGCCAATCCCGGCCGGCTCGCCGGGAGTGTTCACCACACGTTGGGTTTTACCCATTCGTTGCTCAGGATGTTAACGTAGGCGCTTGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCAAGAACTTa                                                                                                                                                                   <  4:380096/261‑1 (MQ=255)
                                                                                                                     tgCAGGGGCGTTGGGGGAATAAGCGCGGCGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTAGGCACGGAGCGTGGCACCGGCCCACACAACCAGTTCGCCCGGCCACTCCCGCTGTGCAAGCTGGCCCACTAAGTCNGGGCCAATCCCGGCCAGCCCGCCGGGAGTCATCCCCACACGTTGCCTTTTCACCATTACTTCCCAACGATTTTGNCGTCGACCCTGGCTCTTTGCTTTTCCAGCCNNCTTTCTTCTTCCTGCgagcattttgcg                                                                                                                                                                >  3:70702/1‑264 (MQ=255)
                                                                                                                               ttGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCa                                                                                                                                                      >  2:307671/1‑273 (MQ=255)
                                                                                                                                 gggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCAAAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCg                                                                                                                                                  <  2:377328/275‑1 (MQ=255)
                                                                                                                                  ggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGc                                                                                                                                                                                                                                                                                >  7:521784/1‑148 (MQ=255)
                                                                                                                                   aggAAGAAGGGCGGAGGGTGAGCGGCAAGCCGAGAATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTGGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCgg                                                                                                                                                 <  2:232256/273‑1 (MQ=255)
                                                                                                                                     gAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                                                                                                                                                                                                                                                              >  5:589112/1‑147 (MQ=255)
                                                                                                                                     gAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                                                                                                                                                                                                                                                              <  6:589112/147‑1 (MQ=255)
                                                                                                                                      aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCggg                                                                                                                                                                                                                                                                            >  6:705782/1‑148 (MQ=255)
                                                                                                                                       aTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGa                                                                                                                                                                                                                                                                           <  8:1044339/148‑1 (MQ=255)
                                                                                                                                              cGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                                                                                                                                                        <  5:504521/144‑1 (MQ=255)
                                                                                                                                              cGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                                                                                                                                                        >  6:504521/1‑144 (MQ=255)
                                                                                                                                                gAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTACCGTAGGCGCTTGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTCCGGTTCATCCGCATGCTGTATCCCCGCTCttt                                                                                                                                   >  4:368639/1‑275 (MQ=255)
                                                                                                                                                      gAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGt                                                                                                                                                                                                                                                            >  5:280837/1‑148 (MQ=255)
                                                                                                                                                      gAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTAGGCACGATCTTTCTgcgc                                                                                                                             >  2:113691/1‑275 (MQ=255)
                                                                                                                                                                   aGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa                                                                                                                                                                                                                                               >  7:28191/1‑148 (MQ=255)
                                                                                                                                                                       tCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa                                                                                                                                                                                                                                               <  6:376648/144‑1 (MQ=255)
                                                                                                                                                                       tCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa                                                                                                                                                                                                                                               >  5:376648/1‑144 (MQ=255)
                                                                                                                                                                         gCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGtt                                                                                                                                                                                                                                         >  5:1010231/1‑148 (MQ=255)
                                                                                                                                                                           tGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                                                                                                                                                                        >  8:622158/1‑147 (MQ=255)
                                                                                                                                                                            gCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCt                                                                                                                                                                                                                                      <  6:1010231/148‑1 (MQ=255)
                                                                                                                                                                             cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATc                                                                                                                                                                                                                                                                     >  7:269118/1‑116 (MQ=255)
                                                                                                                                                                             cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATc                                                                                                                                                                                                                                                                     <  8:269118/116‑1 (MQ=255)
                                                                                                                                                                             cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCCCGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTAc                                                                                                      >  1:9982/1‑275 (MQ=255)
                                                                                                                                                                             cccGGTTGCTACGACGGGTGGCCACGGACCAAACAACCAGTTCGACCGCCCACTCACGCTAGACAAGCTGGACAACTTAGTCCGTGCCAAGCCCGGCGGGCTCGCCGGTGGTGATAACAACACGTTGAGTTTTAACCGTTAAGTGCTCAGGATTTTAACGTAGGCGCTGGCTCGTTGTGCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACGTACGGTTCATCAGCATGCGGTATGCACGATCTTTCGGCGGAGCGTCGGTTTTATCGACATTAc                                                                                                      <  4:18471/275‑1 (MQ=255)
                                                                                                                                                                                    ggTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAttt                                                                                                                                                                                                                              <  8:28191/148‑1 (MQ=255)
                                                                                                                                                                                       aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCACTCCCGGCGAGCCCGCCGGGatt                                                                                                                                                                                                                                                                         >  8:571379/1‑100 (MQ=255)
                                                                                                                                                                                       aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCACTCCCGGCGAGCCCGCCGGGatt                                                                                                                                                                                                                                                                         <  7:571379/102‑3 (MQ=255)
                                                                                                                                                                                       aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                                                                                                                                                                                                 <  6:236154/142‑1 (MQ=255)
                                                                                                                                                                                       aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                                                                                                                                                                                                 >  5:236154/1‑142 (MQ=255)
                                                                                                                                                                                       aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTaa                                                                                                                                                                                                                           <  5:705782/148‑1 (MQ=255)
                                                                                                                                                                                           agagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTa                                                                                                                                                                                                                       <  6:280837/148‑1 (MQ=255)
                                                                                                                                                                                           agagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTa                                                                                                                                                                                                                       >  5:18470/1‑148 (MQ=255)
                                                                                                                                                                                            gagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGt                                                                                        >  4:378803/1‑274 (MQ=255)
                                                                                                                                                                                               tGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatca                                                                                                                                                         >  3:97935/1‑206 (MQ=255)
                                                                                                                                                                                                ggCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATca                                                                                                                                                                                                                                                                    >  8:210770/1‑98 (MQ=255)
                                                                                                                                                                                                  cATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaa                                                                                                                                                                                                                                                                 >  7:1069269/1‑99 (MQ=255)
                                                                                                                                                                                                  cATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaa                                                                                                                                                                                                                                                                 <  8:1069269/99‑1 (MQ=255)
                                                                                                                                                                                                     cGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                                                                                                                                                                                                                                                        <  7:922012/105‑1 (MQ=255)
                                                                                                                                                                                                      ggCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg                                                                                                                                                                                                                                           <  2:433899/117‑1 (MQ=255)
                                                                                                                                                                                                       gCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTa                                                                                                                                                                                                                            <  6:456810/131‑1 (MQ=255)
                                                                                                                                                                                                       gCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCCCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTGGCTGCTTCTTCCGAGCACTTACGGTTCATCAGCATGCGGTTGGCACGCTCTTTCTGCACAGCGTCGGTTTTATCCACATTACGGGTCTCCAGCATCTCGGTTAAAAGc                                                                            >  4:485334/1‑275 (MQ=255)
                                                                                                                                                                                                        cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                                                                                                                                                                                                                                                   <  5:940347/107‑1 (MQ=255)
                                                                                                                                                                                                        cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                                                                                                                                                                                                                                                            >  7:899886/1‑82 (MQ=255)
                                                                                                                                                                                                        cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                                                                                                                                                                                                                                                            <  8:899886/82‑1 (MQ=255)
                                                                                                                                                                                                            aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                                                                                                                                                        >  5:736204/1‑82 (MQ=255)
                                                                                                                                                                                                            aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                                                                                                                                                        <  6:736204/82‑1 (MQ=255)
                                                                                                                                                                                                            aacaacCAGTGCGACCGTCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTGTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTGCGATTAAATGCCAGc                                                                        <  1:70702/274‑1 (MQ=255)
                                                                                                                                                                                                             acaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGcc                                                                                                                                                                                                                                                                               <  5:410660/74‑1 (MQ=255)
                                                                                                                                                                                                              caacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                                                                                                                                                    <  6:991695/148‑1 (MQ=255)
                                                                                                                                                                                                               aacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                                                                                                                                                                                                                                                  >  5:770254/1‑101 (MQ=255)
                                                                                                                                                                                                                  cAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                                                                                                                                                                                                             <  6:104204/103‑1 (MQ=255)
                                                                                                                                                                                                                   aGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTgg                                                                                                                                                                                                               >  7:160927/1‑132 (MQ=255)
                                                                                                                                                                                                                   aGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTgg                                                                                                                                                                                                               <  8:160927/132‑1 (MQ=255)
                                                                                                                                                                                                                           cGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTGAAc                                                        <  2:368640/275‑1 (MQ=255)
                                                                                                                                                                                                                               cACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGTCAATCCCGGCGGGCTCGCCGGGAGTGAACACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTACTGAATCAAGGTTGCTGCTTCCTCCGAGAACTTACGGTTCTTCAGCATGCTGTATGCATGATCTTTTNGCGCAGCGTCGTTTTTATCGACCTTCCGGTTCTTCCGCAGTTCGCTTACATTCCAGTCGGATCAAGAATGGAccc                                                      >  3:118706/1‑272 (MQ=255)
                                                                                                                                                                                                                                 cTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGctt                                                                                                                                                                                 >  5:1071788/1‑148 (MQ=255)
                                                                                                                                                                                                                                    aCGCTGTGCAAGCGGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCa                                                 <  1:44211/273‑1 (MQ=255)
                                                                                                                                                                                                                                        tgtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTAc                                                                                                      >  2:531018/1‑216 (MQ=255)
                                                                                                                                                                                                                                         gtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgtt                                                                                                                                                                                                         <  6:192116/116‑1 (MQ=255)
                                                                                                                                                                                                                                            cAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGt                                                                                                  <  4:531017/216‑1 (MQ=255)
                                                                                                                                                                                                                                                   gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGaa                                                                                                                                                                       <  7:65555/140‑1 (MQ=255)
                                                                                                                                                                                                                                                   gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGaa                                                                                                                                                                       >  8:65555/1‑140 (MQ=255)
                                                                                                                                                                                                                                                    aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCCTTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCa                                                                                                                                                      >  4:556417/1‑156 (MQ=255)
                                                                                                                                                                                                                                                    aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtt                                                                                                                                                                                                                               <  5:920219/83‑1 (MQ=255)
                                                                                                                                                                                                                                                    aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtt                                                                                                                                                                                                                               >  6:920219/1‑83 (MQ=255)
                                                                                                                                                                                                                                                    aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTaa                                                                                                                                                                                                                           <  5:153421/87‑1 (MQ=255)
                                                                                                                                                                                                                                                    aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTaa                                                                                                                                                                                                                           >  6:153421/1‑87 (MQ=255)
                                                                                                                                                                                                                                                    aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTc                               >  1:67664/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGGATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCg                                                                      >  2:209437/1‑236 (MQ=255)
                                                                                                                                                                                                                                                     cAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTAAGAAGAACACGTTGGGTTTGAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGAATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGGGTGAACCGGTGCACTCATTTGACCTTTGTTc                               <  4:307670/274‑1 (MQ=255)
                                                                                                                                                                                                                                                       aCTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcat                                                                                                                                                           >  7:116708/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatc                                                                                                                                                          <  6:1071788/148‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCg                                                                                                                                                  <  2:556418/156‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATg                                                                  <  4:209436/236‑1 (MQ=255)
                                                                                                                                                                                                                                                          aaGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                                                                                                                                                                        <  8:366695/68‑1 (MQ=255)
                                                                                                                                                                                                                                                          aaGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                                                                                                                                                                        >  7:366695/1‑68 (MQ=255)
                                                                                                                                                                                                                                                                    cAATCCCAGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGt               <  2:440465/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCg                                                                                                                                                                           >  5:798148/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                     aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCg                                                                                                                                                                           <  6:798148/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAag                                                               >  2:352870/1‑224 (MQ=255)
                                                                                                                                                                                                                                                                       tCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTGAACCGGTTCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGctcccc            >  4:223415/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                       tCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTACATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTNCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGCCATTACGGGTATCCAGCGGNTCGATTAAATCCCAGCCGACAGAAGAGTGAACCGGTGCACTCATTNGNACTTTGTTCAGGCGCGTCAGGGCGGCAc            >  3:40491/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                       tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                                                                                                                                                                                                             >  6:389783/1‑50 (MQ=37)
                                                                                                                                                                                                                                                                       tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                                                                                                                                                                                                             <  5:389783/50‑1 (MQ=37)
                                                                                                                                                                                                                                                                        cccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                                                                                                                                                    >  6:842748/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                        cccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                                                                                                                                                    <  5:842748/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cccGGCGAGCCCGCCGGGATTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCAc                                                                                                                                          >  5:146580/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                         ccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGgcgc                                                                                                                                                                                                                  <  6:937099/75‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGgcgc                                                                                                                                                                                                                  >  5:937099/1‑75 (MQ=255)
                                                                                                                                                                                                                                                                           ggcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTg                                                           <  4:352869/224‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATc                                                                                                                                      >  6:774906/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                 cTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCTAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGCCCTTTGTTCAGGAGAGTCAGGGCGTCCTGGAATGCCgt  >  2:434029/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                 
CGCTCGCGCCAGCGTTTCCACCACATAATGCCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTCTGCGCAGCGTCGGTTTTATCGACATTACGGGTATCCAGCAGTTCGATTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGG  >  NC_000913/53105‑53652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: