Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,303,984 A→C 7.4% A42A (GCT→GCG 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,9840AC 7.4% 316.6 / 17.0 99A42A (GCT→GCGmdtNmembrane fusion protein of efflux pump
Reads supporting (aligned to +/- strand):  ref base A (40/51);  new base C (5/3);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

GCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAGGGCGGTATAAATAATTCCGACAAATGCCAGATTGATATTTGCCCGCTGAATAACTCGCCTCGTGATGAGCGT  >  NC_000913/4303712‑4304244
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                    
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                                                                                                                                                                                                                                                             >  7:269413/1‑113 (MQ=255)
                                                                                                                                                                       ggAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCTTACGCGTCagca                                                                                                                                                                                                                                                             <  8:269413/116‑4 (MQ=255)
                                                                                                                                                                       ggAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACg                                                                                                                                                                                                                            >  5:73881/1‑149 (MQ=255)
                                                                                                                                                                          aCAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaaa                                                                                                                                                                                                                         >  6:424301/1‑149 (MQ=255)
                                                                                                                                                                             gCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaaaa                                                                                                                                                                                                                        >  8:54686/1‑147 (MQ=255)
                                                                                                                                                                             gCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaaaa                                                                                                                                                                                                                        <  7:54686/147‑1 (MQ=255)
                                                                                                                                                                              cAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCGGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACaa                                                                                                                                                                                                                     <  8:48938/149‑1 (MQ=255)
                                                                                                                                                                                aaTCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   >  6:318057/1‑149 (MQ=255)
                                                                                                                                                                                   cGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCCTGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTagaga                                                                                  >  1:146778/1‑275 (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)
                                                                                                                                                                                                                                                                  aTCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGt                                                                                                                                 <  8:8129/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                                                                                                                                                             >  8:211313/1‑117 (MQ=255)
                                                                                                                                                                                                                                                                      gctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCt                                                                                                                                                             <  7:211313/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCAt                                                                                                                                         <  8:121031/135‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCAt                                                                                                                                         >  7:121031/3‑137 (MQ=255)
                                                                                                                                                                                                                                                                       atgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   >  6:125250/2‑62 (MQ=255)
                                                                                                                                                                                                                                                                       atgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAgg                                                                                                                                                                                                                   <  5:125250/61‑1 (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)
                                                                                                                                                                                                                                                                               tAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTa                                                                                                                    <  6:52115/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                    
GCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTAGAAAAATGCCAGCCAGAACAGCATGGCGTAGAGAGCAAAAAGGGCGGTATAAATAATTCCGACAAATGCCAGATTGATATTTGCCCGCTGAATAACTCGCCTCGTGATGAGCGT  >  NC_000913/4303712‑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: