Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,689 G→T 5.3% intergenic (+40/+62) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6890GT 5.3% 331.7 / 4.5 105intergenic (+40/+62)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (55/43);  new base T (4/2);  total (59/46)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 106 total aligned reads displayed.

CGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGC  >  NC_000913/3608418‑3608950
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                      
cgaccacgcaTTACCCCATACCCCCCTGCGCTGCCTGGTGCACATGCTTGCACCTGCCCCCGCCACTCACGCCAATATCCGCCAGTACATCACGATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGGGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGcc                                                                                                                                                                                                                                                                      <  4:478579/267‑1 (MQ=255)
cgACGCCGCACTCACCCATTGCCACATGCGCGGCCTGGGTCCAATGATTGACCTGGCACGCCCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTTGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGccc                                                                                                                                                                                                                                                                     <  1:81213/275‑1 (MQ=255)
             aCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                                                                                                                       >  2:236254/1‑275 (MQ=255)
                  tAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCCGATCGCAACACTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                                                                                                                                  >  1:210826/1‑275 (MQ=255)
                      caccTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCg                                                                                                                                                                                                                                                                        <  2:278822/250‑1 (MQ=255)
                                                tGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGCTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCCCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCa                                                                                                                                                                                                                     >  1:79776/1‑274 (MQ=255)
                                                                                                       gggCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                                                                                                                                 >  2:506526/1‑191 (MQ=255)
                                                                                                         gCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                                                                                                                                  >  1:56305/1‑173 (MQ=255)
                                                                                                           tGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                                                                                                                             <  4:506524/191‑1 (MQ=255)
                                                                                                             aanGGGATCTTCCTCGTCACCACGCTGTTAGGGANGACCGGGTNNTGGCTGGCAGTGCTGGCAGATACNGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                                                                                                                              <  3:56304/173‑1 (MQ=255)
                                                                                                                                 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCct                                                                                                                                                                                                                                                                   >  6:588417/1‑148 (MQ=255)
                                                                                                                                 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAACAGTTCGTGTTcc                                                                                                                                    >  2:353331/1‑274 (MQ=255)
                                                                                                                                  acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                                                                                                                                 >  7:139364/1‑149 (MQ=255)
                                                                                                                                  acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                                                                                                                                  >  6:432853/1‑148 (MQ=255)
                                                                                                                                     cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                                                                                                                              <  8:816949/149‑1 (MQ=255)
                                                                                                                                       gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                                                                                                                             >  5:754217/1‑148 (MQ=255)
                                                                                                                                       gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                                                                            <  7:265978/149‑1 (MQ=255)
                                                                                                                                        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                                                                               <  8:130433/145‑1 (MQ=255)
                                                                                                                                        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                                                                               >  7:130433/1‑145 (MQ=255)
                                                                                                                                          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                                                                        <  7:580578/149‑1 (MQ=255)
                                                                                                                                           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCTCCTCTGGGGAGAgg                                                                                                                                                                                                                                                       <  7:67490/149‑1 (MQ=255)
                                                                                                                                           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                                                                        <  6:226194/148‑1 (MQ=255)
                                                                                                                                           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                                                                        <  5:319121/148‑1 (MQ=255)
                                                                                                                                             gATGACCGGGTTGTGGTTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                                                                                                                      >  5:580439/1‑148 (MQ=255)
                                                                                                                                             gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                                                                                                                                     >  7:30833/1‑149 (MQ=255)
                                                                                                                                               tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                                                                                                                    <  6:510813/148‑1 (MQ=255)
                                                                                                                                                 aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                                                                                                                                 >  8:112758/1‑149 (MQ=255)
                                                                                                                                                 aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                                                                                                                                 <  8:470791/149‑1 (MQ=255)
                                                                                                                                                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCATCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGAAACCAGTTCGTGTTACATAAGGGTACAAAAcaaa                                                                                                                  >  1:204875/1‑272 (MQ=255)
                                                                                                                                                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                                                                                                <  6:148691/148‑1 (MQ=255)
                                                                                                                                                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                                                                                                                               <  7:644412/149‑1 (MQ=255)
                                                                                                                                                    gggTTGTGGCTGGCAGTGCTGGCAGATGCGGGGGCGACGGTGCTGGTGAAAACGAATGCGATAAGATTGTTGCGCAGGAGATAAGGCAAACCGAGCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGACAggga                                                                                                                                                                                                                                               >  6:689923/1‑147 (MQ=255)
                                                                                                                                                    gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                                                                                                                               >  6:197909/1‑148 (MQ=255)
                                                                                                                                                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                                                                                                                            <  7:238879/149‑1 (MQ=255)
                                                                                                                                                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                                                                                                                            <  7:781725/149‑1 (MQ=255)
                                                                                                                                                      gTTGTGGCTGGCAGTGCTGGCACATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGCTAGGGTGa                                                                                                                                                                                                                                             <  5:783568/148‑1 (MQ=255)
                                                                                                                                                       ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                                                                                                                           <  7:309811/149‑1 (MQ=255)
                                                                                                                                                       ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                                                                                                                           <  8:785094/149‑1 (MQ=255)
                                                                                                                                                             ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                                                                                                                      <  5:805179/148‑1 (MQ=255)
                                                                                                                                                               ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                                                                                                                                 <  7:296398/135‑1 (MQ=255)
                                                                                                                                                               ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                                                                                                                                 >  8:296398/1‑135 (MQ=255)
                                                                                                                                                                      cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                                                                                                                            >  7:365480/1‑149 (MQ=255)
                                                                                                                                                                          cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                                                                                                                                         >  6:737947/1‑148 (MQ=255)
                                                                                                                                                                              taggggggCGAGGGTGCGGGTAACAGCGAAGGCGTAAAGATTGTTGCCACGGAGAAAAGGCAACCGGGCAGCAACATTCGGCGCGACCGGCCCCCCCGCCCCTCTGGGGGGGGGGTTCGGGTGAGGGAAAAAGGCGGCATCGAAGCCaa                                                                                                                                                                                                                    <  8:365480/146‑1 (MQ=255)
                                                                                                                                                                                   gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGCCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAGCGCGGGTAGTGGCCGGATCGTCGCCg                                                                                 >  1:228888/1‑275 (MQ=255)
                                                                                                                                                                                     gggagggtggtGGGGACAGCGAAGGGGTTAAGGTTGTTGCGAGGGAAATAAGGAAAACCAACCCAAACATGGAGGGCATACGTTCCCCTCGCCCCGCTGGGGAGAGGGTGGGGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                                                                                                             <  8:895378/144‑1 (MQ=255)
                                                                                                                                                                                      cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                                                                                                             >  5:366333/1‑148 (MQ=255)
                                                                                                                                                                                       gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                                                                                                           >  8:667718/1‑149 (MQ=255)
                                                                                                                                                                                         cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                                                                                                          >  5:648394/1‑148 (MQ=255)
                                                                                                                                                                                                ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTaagatcgg                                                                                                                                                                                                   >  5:664310/1‑141 (MQ=255)
                                                                                                                                                                                                 gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                                                                                                  >  5:564041/1‑148 (MQ=255)
                                                                                                                                                                                                     cggaaaaTGGGTTAAGTTCTTTGCGGGGGGGAAAGGAAAACCAGATCACAAAATTGAGCGTGATCGGTCCCCTCGCCCCGCGGGGAGAAGGGTTGGGCGGGGAGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                                                                                                                             <  7:112758/144‑1 (MQ=255)
                                                                                                                                                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCCGATCGGAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGCATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAAAGtttag                                                               >  1:72956/1‑272 (MQ=255)
                                                                                                                                                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCCATCAGCCCCTATCAACCGCCTTTACGAATCAACTCACGATAAGGCAGTCCATCCGTCTCGTCAGCAACCAGTTCGTGTTCCAGAAAGGTACAAAACCCAGGAATATCTCGGGTGGTGGCCGGATCGACGGGGATACACGGCACAGTTTCg                                                               >  4:76196/1‑274 (MQ=255)
                                                                                                                                                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTACCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCaaa                                                                                                                                                                                                                    >  2:395310/1‑125 (MQ=255)
                                                                                                                                                                                                        caaaTGCTTTAAGATTGTGGCGAAGGAGAAAAGGAAAACCGATCGCAACATTGGGCGCGACCGCCCCCCCCCCCCTTGGGGGGAGGGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                                                                                                           <  5:737947/146‑1 (MQ=255)
                                                                                                                                                                                                        cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGcc                                                             >  2:709386/1‑274 (MQ=255)
                                                                                                                                                                                                        cGAACGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCCACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGCCAGGCCATCCGTCTCTTTAGCAACCAGTTCTTGTTCCATAACGTTACAAACCCCAGGAATATAGCGGGTAGTGGCCGGATCGTCGGCGATAATCAAAAAAGttttccc                                                             >  1:80825/1‑270 (MQ=255)
                                                                                                                                                                                                         gAATGCTTTAAGTTTGTGGCGCAGGAGAAAAGGAAACCCGACCGAAAAATGGAGCGCGATCGGCCCCCCCGCCCCTCGGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                                                                                                          <  5:432853/148‑1 (MQ=255)
                                                                                                                                                                                                          aaTGCGTTAAGTTTGTTGCGAAGGAAATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCCCGCCCCTCTGGGGAGGGGGTTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                                                                                                                                        <  8:823348/149‑1 (MQ=255)
                                                                                                                                                                                                                ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                                                                                                   <  6:651202/100‑1 (MQ=255)
                                                                                                                                                                                                                     aTTTTTGCGAAGAAAATAAGGCAAACCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                                                                                                                                             <  7:884661/149‑1 (MQ=38)
                                                                                                                                                                                                                         ttgCGCAGGAGATAAGGAAAACGGAGCGAAACATTGCGCGCGATCGGCCCCCCCGCCCCTCGGGGGAGGGGGTTGGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                                                                                                                         <  8:139364/149‑1 (MQ=255)
                                                                                                                                                                                                                         tcgggaaGGAGAAAAGAAAAACCAAACAAAATATGGCGCGCGAGCGCCCCCCCCGGCCCTCGGGGGGGGGGGTGAGGGTGAGGAGAAAAGGCGGCCTCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                                                                                                                          <  5:831349/142‑1 (MQ=255)
                                                                                                                                                                                                                             gcAGGAAAGAAGGAAAACCGATCGAAACATGGAGCGCGATCGGCCCCCTCGCCCCTCGGGGGAGAGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                                                                                                                     <  8:30833/149‑1 (MQ=255)
                                                                                                                                                                                                                                     tAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCATTTACGCATCTAATAACGATAAGGCAGTCCATCCGt                                                                                                                                                             >  8:560097/1‑149 (MQ=39)
                                                                                                                                                                                                                                      aaGGCAAACCGATCGCAACATTGGGCGCGATCGGCCCCCTCGCCCTTCGGGGGAGAGGTTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                                                                                                                             <  5:60639/148‑1 (MQ=255)
                                                                                                                                                                                                                                          cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                                                                                                         <  6:754217/148‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaCGGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGCGAGGGTTAGGGGGAGGTGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACGGCCTTACCGAATCAAATAAAGATAAAGCAGTCCANCCGGCTCTTTCGAAAGCAGTTCGTGTTCAACACAGATACAGAACCCAGGAAAAGCGAGGGTAGAGGACGGAAGANCGGCGACAAACAGCAGCATTTCCGCACGATNNATATTGCGGCCGGGGTTGGGCAcccacac                         >  3:101161/1‑270 (MQ=255)
                                                                                                                                                                                                                                            aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                                                                                                       >  5:458674/1‑148 (MQ=255)
                                                                                                                                                                                                                                            aaCCGATCGAAACATTGAGCGGGATCGGCCCCCCCGCCCCTCGGGGGGGGGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                                                                                                       <  6:366333/148‑1 (MQ=255)
                                                                                                                                                                                                                                             cccGATCAAAATATTGGGGGCGAGCGGCCCCCCGGCCCCTGTGGGGAGGGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                                                                                                      <  6:315070/147‑1 (MQ=255)
                                                                                                                                                                                                                                                  aTCGCACAATGGAGCGCGATCGGCCCCCCCGCCCCTCGGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                                                                                                 <  7:560097/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                                                                                                                             >  7:192878/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        actttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                                                                                                            <  6:564041/145‑1 (MQ=255)
                                                                                                                                                                                                                                                          caatgttgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                                                                                                                                          >  5:845521/8‑148 (MQ=255)
                                                                                                                                                                                                                                                            tGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAg         >  2:309407/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTCTCAACCGCCTTTACGAATCAAATCACGATAATGCAGTCCTTCTTTCTCTCTATCAACCAGTTCGTGTTCCCGTCATGCACAAAACCCAGGAATCTCGCGGGTAGCGCCCGGTTCGGCGACAAAAACCCGCCACGTTTCCCCAGGCTCTAAGTTGCCCACCGTTTTGAGGACAACCACCACCTGTTCCGGGCAgc        >  4:358789/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                                                                                                                      >  6:910059/1‑52 (MQ=255)
                                                                                                                                                                                                                                                              aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                                                                                                                      <  5:910059/52‑1 (MQ=255)
                                                                                                                                                                                                                                                               gcgcGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                                                                                                    >  7:538243/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 gcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAACGCGGGATGGGGCCGGATCGTCGGCGATAAACAGCAACGTTTCGCCAGGCGGCATATTGCGCACGGTCTTGCGCCCCATCTTCCCCGTTTCCcggaccct       >  4:243145/1‑267 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGGGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                                                                                                >  8:479188/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                                                                                                >  7:659032/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGc  >  2:183559/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGAATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCCGCGCAGGc  >  2:584008/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                    aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                                                                                                >  6:280899/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    aTCGGTCCCCTCGCCCCTCAGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAGTAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCACTTCGTGTTCCATaa                                                                                                                                >  6:81827/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     gcGGTCCCCTCGCCCCCATGGGGAAGGGGTTAGGGATAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                                                                                                               <  5:689923/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                                                                                             >  5:381634/5‑148 (MQ=255)
                                                                                                                                                                                                                                                                        taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                                                                                                            >  6:379286/4‑148 (MQ=255)
                                                                                                                                                                                                                                                                        taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                                                                                                                           >  7:722052/4‑149 (MQ=255)
                                                                                                                                                                                                                                                                         tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                                                                                                                           >  5:908689/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          acccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                                                                                                          >  6:587681/2‑148 (MQ=255)
                                                                                                                                                                                                                                                                           cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                                                                                                               <  8:821723/142‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                                                                                                       >  7:75240/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                            ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                                                                                                       <  8:75240/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                                                                                                                        >  5:907364/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa                                                                                                                      <  7:667718/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                                                                                                     >  5:387282/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                      
CGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGC  >  NC_000913/3608418‑3608950

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