Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,303,981 G→T 7.5% Y43* (TAC→TAA mdtN ← membrane fusion protein of efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,303,9810GT 7.5% 312.0 / 12.2 95Y43* (TAC→TAAmdtNmembrane fusion protein of efflux pump
Reads supporting (aligned to +/- strand):  ref base G (39/48);  new base T (4/3);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 101 total aligned reads displayed.

CTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAGGGCGGTATAAATAATTCCGACAAATGCCAGATTGATATTTGCCCGCTGAATAACTCGCCTCGTGATGAGCGT  >  NC_000913/4303707‑4304244
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                       
cTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGATTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACCATGCGGCCGCTGACTGCCGGTACCACATCCATGGTATCTGCTGGcgcg                                                                                                                                                                                                                                                                         >  1:87541/1‑275 (MQ=255)
     gCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGc                                                                                                                                                                                                                                                                    >  2:271314/1‑275 (MQ=255)
         gTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACCACCATGGTATCTGCTGACGCGTAAGCGTc                                                                                                                                                                                                                                                                 >  2:343434/1‑274 (MQ=255)
               cggcggCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCa                                                                                                                                                                                                                                                                >  2:417572/1‑269 (MQ=255)
                      aCGGGCTTTTTCTACCGTTGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTCCGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGgcg                                                                                                                                                                                                                                                   <  2:188088/275‑1 (MQ=255)
                         ggtttttttTACCCGTGGCATTCCCCGAGTCGGCACAAAAATGTTGCGTGTGAACGCTACGCTGGGTTAGCATAATGTGCTTATCCAGCGCCGCGTGGGAGGTTTCAGCTTTCGCCAGATTGGCTGCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCa                                                                                                                                                                                                                                                 <  4:144265/272‑1 (MQ=255)
                                             ttAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCCTCGTACGGGTCCCGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCTCCAGTTCTACAATGCGTCCGCTGACTTCCGGTCACATACCAATGGACTCTGCTGGCCCATACGCCGCATTATGTGAATGCGCCCTTTCTACGCCCCCGATaa                                                                                                                                                                                                                              >  4:475963/1‑274 (MQ=255)
                                                                        gcgTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCAc                                                                                                                                                                                                                                                <  2:67042/229‑1 (MQ=255)
                                                                        gcgTCAACGCTACGCTGGGTGAGCATAATTTGCTGATCAAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACg                                                                                                                                                                                                   <  2:418608/274‑1 (MQ=255)
                                                                         cgTGAACGGTAGGGTGTGTGAGCAGAGTTTGGTGTTCCGGGCCCGGGAGGGGGGCTGCAGCGTTCGCCAGTTTGGGTGGGTACGGGCGCGGGTCGATGCGGAACAGGAAACCGCCCTGTTGGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGGCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGGAAGCGTCATTAGTGGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCa                                                                                                                                                                                                 <  4:271310/275‑1 (MQ=255)
                                                                         cgTCAACGCTACGCTGGGTGAGCATTATGGGCTTATCCAGCTCCGAGAGGGAGGCTTCAGCTTTCGCCAGTTTGGCTTCGTACGGGCGCGGGTCAATGCGGAGCAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCa                                                                                                                                                                                                 <  2:331737/275‑1 (MQ=255)
                                                                            cAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAAcc                                                                                                                                                                                              >  2:141021/1‑275 (MQ=255)
                                                                                 tgacggtGGGTGAGAAGAATTGGTTTTCCCGGCGCGGCGAGGGAGGGTTCAACTTTTGCCAGATTGGCTTGGTACGGGCGCGGTTCGGTGCGGAACGGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCATTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTaa                                                                                                                                                                                         <  4:141020/273‑1 (MQ=255)
                                                                                 ctacgctGGGTGAGCATAATTTGCTTAGCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTaa                                                                                                                                                                                         <  2:426440/275‑1 (MQ=255)
                                                                                                   aTTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACt                                                                                                                                                                       <  2:475968/275‑1 (MQ=255)
                                                                                                         ttATCCAGCGACGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGAGGGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCgag                                                                                                                                                                  <  2:288252/274‑1 (MQ=255)
                                                                                                                       gagtgaggCTGCAGCTTTTGCCAGATTGGCTTCGTACTGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCGGGTGGTCGGTGACCGGCAGTTCTACAATGCGGCAGCTGACTTCCGGTACCACATCAAGGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAAACAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCg                                                                                                                                                   <  1:43524/275‑1 (MQ=255)
                                                                                                                                  cAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTaa                                                                                                                                                                                                                                                                      <  8:425822/149‑1 (MQ=255)
                                                                                                                                   aGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAg                                                                                                                                                                                                                                                                     <  7:438920/149‑1 (MQ=255)
                                                                                                                                     cTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAgc                                                                                                                                     <  2:180501/275‑1 (MQ=255)
                                                                                                                                     cTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAgc                                                                                                                                     <  1:3653/275‑1 (MQ=255)
                                                                                                                                      tttCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACACTGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGt                                                                                                                                                                                                                                                                  >  5:81984/1‑149 (MQ=255)
                                                                                                                                        tCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCCGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCTCCAGATAACGAAAACAAGGTCAACCAGCGCCAACGCAACCACTACCACAGCCGGGAATTTACTTCCAGGAGCTTTTTTCGGCCTACTTTCCCTGAGcgcg                                                                                                                                  >  4:35698/1‑275 (MQ=255)
                                                                                                                                          gCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAg                                                                                                                                                                                                                                                                     >  7:223521/1‑142 (MQ=255)
                                                                                                                                          gCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAg                                                                                                                                                                                                                                                                     <  8:223521/142‑1 (MQ=255)
                                                                                                                                           ccAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTa                                                                                                                                                                                                                                                             >  5:52115/1‑149 (MQ=255)
                                                                                                                                               aTTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCaaa                                                                                                                           <  2:391278/275‑1 (MQ=255)
                                                                                                                                               aTTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCAGTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCaaa                                                                                                                           <  2:35698/275‑1 (MQ=255)
                                                                                                                                                   gCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATgg                                                                                                                                                                                                                                                     >  7:48938/1‑149 (MQ=255)
                                                                                                                                                     ttCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGgcg                                                                                                                                                                                                                                                   <  7:471129/149‑1 (MQ=255)
                                                                                                                                                           cGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGt                                                                                                                                                                                                                                             <  8:338166/149‑1 (MQ=255)
                                                                                                                                                                      cGATGCGGAACCGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGccagccag                                                                                                    <  2:381369/275‑1 (MQ=255)
                                                                                                                                                                           cGGCACAACAAATGGCCCTGTGTGTCTGCCGGTTTGTCGATTACGGGCGGTTATACGATGAGGCCGCTGACTGTCGGTACCACATCAATGTTATCTGCGGATGCGTAATCGTCAGTAGGTGATGGCGCACTGTCCACGCGCCAGAGAACGAAAACAAGGGCAAGCAGCGCCAACGAAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTGTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGGCAGCCAGAACa                                                                                                <  4:200846/274‑1 (MQ=255)
                                                                                                                                                                           cGGAACAGCAAATCGCCCGGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACa                                                                                                <  2:489327/274‑1 (MQ=255)
                                                                                                                                                                            ggAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCTTACGCGTCagca                                                                                                                                                                                                                                                             <  8:269413/116‑4 (MQ=255)
                                                                                                                                                                            ggAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCTTACGCGTCagca                                                                                                                                                                                                                                                             >  7:269413/1‑113 (MQ=255)
                                                                                                                                                                            ggAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACg                                                                                                                                                                                                                            >  5:73881/1‑149 (MQ=255)
                                                                                                                                                                               aCAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaaa                                                                                                                                                                                                                         >  6:424301/1‑149 (MQ=255)
                                                                                                                                                                                  gCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaaaa                                                                                                                                                                                                                        <  7:54686/147‑1 (MQ=255)
                                                                                                                                                                                  gCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaaaa                                                                                                                                                                                                                        >  8:54686/1‑147 (MQ=255)
                                                                                                                                                                                   cAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCGGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACaa                                                                                                                                                                                                                     <  8:48938/149‑1 (MQ=255)
                                                                                                                                                                                     aaTCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   >  6:318057/1‑149 (MQ=255)
                                                                                                                                                                                        cGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCCTGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTagaga                                                                                  >  1:146778/1‑275 (MQ=255)
                                                                                                                                                                                          cccTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGt                                                                                                                                                                                                                                                                        >  2:335998/1‑91 (MQ=255)
                                                                                                                                                                                           ccTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAAcc                                                                                                                                                                                                             >  5:197689/1‑149 (MQ=255)
                                                                                                                                                                                           ccTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAAcc                                                                                                                                                                                                             >  5:76189/1‑149 (MQ=255)
                                                                                                                                                                                             tGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaa                                                                                                                                                                                                                          <  6:53076/134‑1 (MQ=255)
                                                                                                                                                                                             tGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaa                                                                                                                                                                                                                          >  5:53076/1‑134 (MQ=255)
                                                                                                                                                                                              gTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGc                                                                                                                                                                                                                                                                    <  4:335994/91‑1 (MQ=255)
                                                                                                                                                                                               tttGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATATGCTGACGCGTAAGCGGCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGcc                                                                                                                                                                                                       >  3:150027/1‑151 (MQ=255)
                                                                                                                                                                                                 tggcggtctttttGTCGGTTGCCGCCAGTNNNACAATGCATCCGCTGTATTCAGGGACCACATCAATGGTATGTGCTGAGGGGTAATGGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCATCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCATACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAgg                                                                         <  3:146778/264‑1 (MQ=255)
                                                                                                                                                                                                 ggccTGCTTGGTTAGCGGTGTCCGCCAGTTCTAAAATGCGGCCGCTGACTTCCGGGACCACATCAATGGTATCTGCTGACGAGTAAGAGTCATTAGTTGATGGCGCACGGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGGGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGTAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAgg                                                                         <  4:340439/272‑1 (MQ=255)
                                                                                                                                                                                                   aCTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCaaa                                                                                                                                                                                                     <  7:206555/149‑1 (MQ=255)
                                                                                                                                                                                                   aCTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACg                                                                                                                                                                                                   <  1:150027/151‑1 (MQ=255)
                                                                                                                                                                                                    cTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAAc                                                                                                                                                                                                    >  7:194560/1‑149 (MQ=255)
                                                                                                                                                                                                     tGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGgcg                                                                                                                                                                                                                                                   >  2:35421/1‑101 (MQ=255)
                                                                                                                                                                                                         tGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACt                                                                                                                                                                                                                                               <  4:35421/101‑1 (MQ=255)
                                                                                                                                                                                                         tGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAAc                                                                                                                                                                                               >  8:206638/1‑149 (MQ=255)
                                                                                                                                                                                                         tGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCaaa                                                                             >  2:340443/1‑263 (MQ=255)
                                                                                                                                                                                                           gTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCa                                                                                                                                                                                             >  8:341497/1‑149 (MQ=255)
                                                                                                                                                                                                                  gTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   <  5:222945/120‑1 (MQ=255)
                                                                                                                                                                                                                  gTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   >  6:222945/1‑120 (MQ=255)
                                                                                                                                                                                                                     aCCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGc                                                                                                                                                                                   <  6:76189/149‑1 (MQ=255)
                                                                                                                                                                                                                      ccGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCa                                                                                                                                                                                  <  5:424301/149‑1 (MQ=255)
                                                                                                                                                                                                                       cGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAg                                                                                                                                                                                 >  5:401938/1‑149 (MQ=255)
                                                                                                                                                                                                                           aGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGaa                                                                                                                                                                             <  6:81984/149‑1 (MQ=255)
                                                                                                                                                                                                                             ttCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCAc                                                                                                                                                                                                                                                <  6:308638/80‑1 (MQ=255)
                                                                                                                                                                                                                             ttCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCAGACGCGTAAGCGTCATTAGTTGATGGCGCAc                                                                                                                                                                                                                                                >  5:308638/1‑80 (MQ=255)
                                                                                                                                                                                                                                 aCAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACt                                                                                                                                                                       <  5:71212/149‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGc                                                                                                                                                                     <  5:318057/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATgg                                                                                                                                                                                                                                                     >  3:131924/1‑65 (MQ=255)
                                                                                                                                                                                                                                           cGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCa                                                                                                                                                                                                                                                 <  1:131924/65‑1 (MQ=255)
                                                                                                                                                                                                                                                aCTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCtttttt                                                                                                                                                        <  6:73881/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACaga                                                                                                                                                                                     >  5:17417/1‑114 (MQ=255)
                                                                                                                                                                                                                                                      ggTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACaga                                                                                                                                                                                     <  6:17417/114‑1 (MQ=255)
                                                                                                                                                                                                                                                       gTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTa                                                                                                                                                 >  5:25685/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            acaTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTc                                                                                                                                            <  6:401938/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            acaTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCAt                                                                                                                                         >  1:126800/1‑152 (MQ=255)
                                                                                                                                                                                                                                                              aTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCa                                                                                                                                          >  7:8129/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               tCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCAt                                                                                                                                         >  5:390354/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                cAATTGTATCTGCNGAAGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAgc                                                                                                                                     <  3:126800/152‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGATTTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAgcgc                                                                                                                                   >  8:309096/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     gTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAgcgc                                                                                                                                   <  5:249172/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAGGGCGGTATAAATAATTCCGACAAATGCCAGATTGATATTTGCCCGCTGAATAACCAGCCTCGTGATGAACGt  >  4:16050/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                        tCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAGGGCGGTATAAATAATTCCGACAAATGCCAGATTGATATTTGCCCGCGGAATAACTCGCCTCGGGATGAGCg   >  2:55884/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          taatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGtctc                                                                                                                              >  8:21464/4‑149 (MQ=255)
                                                                                                                                                                                                                                                                           gctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCt                                                                                                                                                             <  7:211313/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCt                                                                                                                                                             >  8:211313/1‑117 (MQ=255)
                                                                                                                                                                                                                                                                           aatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCAt                                                                                                                                         <  8:121031/135‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCAt                                                                                                                                         >  7:121031/3‑137 (MQ=255)
                                                                                                                                                                                                                                                                            atgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   <  5:125250/61‑1 (MQ=255)
                                                                                                                                                                                                                                                                            atgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   >  6:125250/2‑62 (MQ=255)
                                                                                                                                                                                                                                                                             tgACGCGTAAGCGTCATTAGTTGATGTCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCaaa                                                                                                                           >  7:435413/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               aCGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACa                                                                                                                                                                                       >  2:380815/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                               aCGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATa                                                                                                                         <  7:341497/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gcgTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGc                                                                                                                                                              <  5:303961/110‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gcgTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGc                                                                                                                                                              >  6:303961/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                                 gcgTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGa                                                                                                                       <  6:197689/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gcgTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGa                                                                                                                       >  6:454088/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                   gTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGc                                                                                                                                                                                   <  4:380810/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                       
CTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAGGGCGGTATAAATAATTCCGACAAATGCCAGATTGATATTTGCCCGCTGAATAACTCGCCTCGTGATGAGCGT  >  NC_000913/4303707‑4304244

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