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.

CGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGC  >  NC_000913/3608418‑3608950
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                      
cgaccacgcaTTACCCCATACCCCCCTGCGCTGCCTGGTGCACATGCTTGCACCTGCCCCCGCCACTCACGCCAATATCCGCCAGTACATCACGATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGGGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGcc                                                                                                                                                                                                                                                                      <  4:478579/267‑1 (MQ=255)
cgACGCCGCACTCACCCATTGCCACATGCGCGGCCTGGGTCCAATGATTGACCTGGCACGCCCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTTGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGccc                                                                                                                                                                                                                                                                     <  1:81213/275‑1 (MQ=255)
             aCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑Agg                                                                                                                                                                                                                                                       >  2:236254/1‑275 (MQ=255)
                  tAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCC‑GATCGCAACACTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAg                                                                                                                                                                                                                                                  >  1:210826/1‑275 (MQ=255)
                      caccTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCg                                                                                                                                                                                                                                                                        <  2:278822/250‑1 (MQ=255)
                                                tGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGCTTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCCCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCa                                                                                                                                                                                                                     >  1:79776/1‑274 (MQ=255)
                                                                                                       gggCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                                                                                                                                 >  2:506526/1‑191 (MQ=255)
                                                                                                         gCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                                                                                                                                  >  1:56305/1‑173 (MQ=255)
                                                                                                           tGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGa                                                                                                                                                                                                                                             <  4:506524/191‑1 (MQ=255)
                                                                                                             aanGGGATCTTCCTCGTCACCACGCTGTTAGGGANGACCGGGTNNTGGCTGGCAGTGCTGGCAGATACNGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                                                                                                                              <  3:56304/173‑1 (MQ=255)
                                                                                                                                 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCct                                                                                                                                                                                                                                                                   >  6:588417/1‑148 (MQ=255)
                                                                                                                                 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAACAGTTCGTGTTcc                                                                                                                                    >  2:353331/1‑274 (MQ=255)
                                                                                                                                  acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                                                                                                                                 >  7:139364/1‑149 (MQ=255)
                                                                                                                                  acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                                                                                                                                  >  6:432853/1‑148 (MQ=255)
                                                                                                                                     cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                                                                                                                              <  8:816949/149‑1 (MQ=255)
                                                                                                                                       gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                                                                                                                             >  5:754217/1‑148 (MQ=255)
                                                                                                                                       gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                                                                            <  7:265978/149‑1 (MQ=255)
                                                                                                                                        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                                                                               <  8:130433/145‑1 (MQ=255)
                                                                                                                                        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                                                                               >  7:130433/1‑145 (MQ=255)
                                                                                                                                          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGag‑ag                                                                                                                                                                                                                                                        <  7:580578/149‑1 (MQ=255)
                                                                                                                                           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCTCCTCTGGGGAG‑Agg                                                                                                                                                                                                                                                       <  7:67490/149‑1 (MQ=255)
                                                                                                                                           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGag‑ag                                                                                                                                                                                                                                                        <  6:226194/148‑1 (MQ=255)
                                                                                                                                           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGag‑ag                                                                                                                                                                                                                                                        <  5:319121/148‑1 (MQ=255)
                                                                                                                                             gATGACCGGGTTGTGGTTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑Aggg                                                                                                                                                                                                                                                      >  5:580439/1‑148 (MQ=255)
                                                                                                                                             gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGt                                                                                                                                                                                                                                                     >  7:30833/1‑149 (MQ=255)
                                                                                                                                               tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGtt                                                                                                                                                                                                                                                    <  6:510813/148‑1 (MQ=255)
                                                                                                                                                 aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                                                                                                                                 >  8:112758/1‑149 (MQ=255)
                                                                                                                                                 aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                                                                                                                                 <  8:470791/149‑1 (MQ=255)
                                                                                                                                                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCATCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGAAACCAGTTCGTGTTACATAAGGGTACAAAAcaaa                                                                                                                  >  1:204875/1‑272 (MQ=255)
                                                                                                                                                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAggg                                                                                                                                                                                                                                                <  6:148691/148‑1 (MQ=255)
                                                                                                                                                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGt                                                                                                                                                                                                                                               <  7:644412/149‑1 (MQ=255)
                                                                                                                                                    gggTTGTGGCTGGCAGTGCTGGCAGATGCGGGGGCGACGGTGCTGGTGAAAACGAATGCGATAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GAGCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGACAggga                                                                                                                                                                                                                                               >  6:689923/1‑147 (MQ=255)
                                                                                                                                                    gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGt                                                                                                                                                                                                                                               >  6:197909/1‑148 (MQ=255)
                                                                                                                                                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAg                                                                                                                                                                                                                                            <  7:238879/149‑1 (MQ=255)
                                                                                                                                                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAg                                                                                                                                                                                                                                            <  7:781725/149‑1 (MQ=255)
                                                                                                                                                      gTTGTGGCTGGCAGTGCTGGCACATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGCTAGGGTGa                                                                                                                                                                                                                                             <  5:783568/148‑1 (MQ=255)
                                                                                                                                                       ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAgg                                                                                                                                                                                                                                           <  7:309811/149‑1 (MQ=255)
                                                                                                                                                       ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAgg                                                                                                                                                                                                                                           <  8:785094/149‑1 (MQ=255)
                                                                                                                                                             ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGaaa                                                                                                                                                                                                                                      <  5:805179/148‑1 (MQ=255)
                                                                                                                                                               ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                                                                                                                                 <  7:296398/135‑1 (MQ=255)
                                                                                                                                                               ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                                                                                                                                 >  8:296398/1‑135 (MQ=255)
                                                                                                                                                                      cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                                                                                                                            >  7:365480/1‑149 (MQ=255)
                                                                                                                                                                          cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                                                                                                                                         >  6:737947/1‑148 (MQ=255)
                                                                                                                                                                              taggggggCGAGGGTGCGGGTAACAGCGAAGGCGTAAAGATTGTTGCCACGGAGAAAAGGCAACCG‑GGCAGCAACATTCGGCGCGACCGGCCCCCCCGCCCCTCTGGGGGG‑GGGGTTCGGGTGAGGGAAAAAGGCGGCATCGAAGCCaa                                                                                                                                                                                                                    <  8:365480/146‑1 (MQ=255)
                                                                                                                                                                                   gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGCCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAGCGCGGGTAGTGGCCGGATCGTCGCCg                                                                                 >  1:228888/1‑275 (MQ=255)
                                                                                                                                                                                     gggagggtggtGGGGACAGCGAAGGGGTTAAGGTTGTTGCGAGGGAAATAAGGAAAACC‑AACCCAAACATGGAGGGCATACGTTCCCCTCGCCCCGCTGGGGAG‑AGGGTGGGGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                                                                                                             <  8:895378/144‑1 (MQ=255)
                                                                                                                                                                                      cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                                                                                                             >  5:366333/1‑148 (MQ=255)
                                                                                                                                                                                       gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                                                                                                           >  8:667718/1‑149 (MQ=255)
                                                                                                                                                                                         cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                                                                                                          >  5:648394/1‑148 (MQ=255)
                                                                                                                                                                                                ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTaagatcgg                                                                                                                                                                                                   >  5:664310/1‑141 (MQ=255)
                                                                                                                                                                                                 gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                                                                                                  >  5:564041/1‑148 (MQ=255)
                                                                                                                                                                                                     cggaaaaTGGGTTAAGTTCTTTGCGGGGGGGA‑AAGGAAAACCAGATCACAAAATTGAGCGTGATCGGTCCCCTCGCCCCGC‑GGGGAGAAGGGTTGGGCGGGGAGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                                                                                                                             <  7:112758/144‑1 (MQ=255)
                                                                                                                                                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCC‑GATCGGAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGCATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAAAGtttag                                                               >  1:72956/1‑272 (MQ=255)
                                                                                                                                                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCCATCAGCCCCTATCAACCGCCTTTACGAATCAACTCACGATAAGGCAGTCCATCCGTCTCGTCAGCAACCAGTTCGTGTTCCAGAAAGGTACAAAACCCAGGAATATCTCGGGTGGTGGCCGGATCGACGGGGATACACGGCACAGTTTCg                                                               >  4:76196/1‑274 (MQ=255)
                                                                                                                                                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTACCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCaaa                                                                                                                                                                                                                    >  2:395310/1‑125 (MQ=255)
                                                                                                                                                                                                        caaaTGCTTTAAGATTGTGGCGAAGGAGAAAAGGAAAACC‑GATCGCAACATTGGGCGCGACCGCCCCCCCCCCCCTTGGGGGGAG‑GGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                                                                                                           <  5:737947/146‑1 (MQ=255)
                                                                                                                                                                                                        cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGcc                                                             >  2:709386/1‑274 (MQ=255)
                                                                                                                                                                                                        cGAACGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCCACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGCCAGGCCATCCGTCTCTTTAGCAACCAGTTCTTGTTCCATAACGTTACAAACCCCAGGAATATAGCGGGTAGTGGCCGGATCGTCGGCGATAATCAAAAAAGttttccc                                                             >  1:80825/1‑270 (MQ=255)
                                                                                                                                                                                                         gAATGCTTTAAGTTTGTGGCGCAGGAGAAAAGGAAACCC‑GACCGAAAAATGGAGCGCGATCGGCCCCCCCGCCCCTCGGGGGAG‑AGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                                                                                                          <  5:432853/148‑1 (MQ=255)
                                                                                                                                                                                                          aaTGCGTTAAGTTTGTTGCGAAGGAAATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGCCCCCCCGCCCCTCTGGGGAG‑GGGGTTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                                                                                                                                        <  8:823348/149‑1 (MQ=255)
                                                                                                                                                                                                                ttAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                                                                                                   <  6:651202/100‑1 (MQ=255)
                                                                                                                                                                                                                     aTTTTTGCGAAGAAAATAAGGCAAACC‑GATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑GGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                                                                                                                                             <  7:884661/149‑1 (MQ=38)
                                                                                                                                                                                                                         ttgCGCAGGAGATAAGGAAAACG‑GAGCGAAACATTGCGCGCGATCGGCCCCCCCGCCCCTCGGGGGAG‑GGGGTTGGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                                                                                                                         <  8:139364/149‑1 (MQ=255)
                                                                                                                                                                                                                         tcgggaaGGAGAAAAGAAAAACC‑AAACAAAATATGGCGCGCGAGCGCCCCCCCCGGCCCTCGGGGGGG‑GGGGTGAGGGTGAGGAGAAAAGGCGGCCTCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                                                                                                                          <  5:831349/142‑1 (MQ=255)
                                                                                                                                                                                                                             gcAGGAAAGAAGGAAAACC‑GATCGAAACATGGAGCGCGATCGGCCCCCTCGCCCCTCGGGGGAG‑AGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                                                                                                                     <  8:30833/149‑1 (MQ=255)
                                                                                                                                                                                                                                     tAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCATTTACGCATCTAATAACGATAAGGCAGTCCATCCGt                                                                                                                                                             >  8:560097/1‑149 (MQ=39)
                                                                                                                                                                                                                                      aaGGCAAACC‑GATCGCAACATTGGGCGCGATCGGCCCCCTCGCCCTTCGGGGGAG‑AGGTTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                                                                                                                             <  5:60639/148‑1 (MQ=255)
                                                                                                                                                                                                                                          cAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                                                                                                         <  6:754217/148‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaCG‑GATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGCG‑AGGGTTAGGGGGAGGTGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACGGCCTTACCGAATCAAATAAAGATAAAGCAGTCCANCCGGCTCTTTCGAAAGCAGTTCGTGTTCAACACAGATACAGAACCCAGGAAAAGCGAGGGTAGAGGACGGAAGANCGGCGACAAACAGCAGCATTTCCGCACGATNNATATTGCGGCCGGGGTTGGGCAcccacac                         >  3:101161/1‑270 (MQ=255)
                                                                                                                                                                                                                                            aaCC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                                                                                                       >  5:458674/1‑148 (MQ=255)
                                                                                                                                                                                                                                            aaCC‑GATCGAAACATTGAGCGGGATCGGCCCCCCCGCCCCTCGGGGGGG‑GGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                                                                                                       <  6:366333/148‑1 (MQ=255)
                                                                                                                                                                                                                                             ccc‑GATCAAAATATTGGGGGCGAGCGGCCCCCCGGCCCCTGTGGGGAG‑GGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                                                                                                      <  6:315070/147‑1 (MQ=255)
                                                                                                                                                                                                                                                  aTCGCACAATGGAGCGCGATCGGCCCCCCCGCCCCTCGGGGGAG‑AGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                                                                                                 <  7:560097/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                                                                                                                             >  7:192878/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        actttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                                                                                                            <  6:564041/145‑1 (MQ=255)
                                                                                                                                                                                                                                                          caatgttgcgATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                                                                                                                                          >  5:845521/8‑148 (MQ=255)
                                                                                                                                                                                                                                                            tGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAg         >  2:309407/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTCTCAACCGCCTTTACGAATCAAATCACGATAATGCAGTCCTTCTTTCTCTCTATCAACCAGTTCGTGTTCCCGTCATGCACAAAACCCAGGAATCTCGCGGGTAGCGCCCGGTTCGGCGACAAAAACCCGCCACGTTTCCCCAGGCTCTAAGTTGCCCACCGTTTTGAGGACAACCACCACCTGTTCCGGGCAgc        >  4:358789/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGaaa                                                                                                                                                                                                                                      >  6:910059/1‑52 (MQ=255)
                                                                                                                                                                                                                                                              aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGaaa                                                                                                                                                                                                                                      <  5:910059/52‑1 (MQ=255)
                                                                                                                                                                                                                                                               gcgcGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                                                                                                    >  7:538243/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 gcgATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAACGCGGGATGGGGCCGGATCGTCGGCGATAAACAGCAACGTTTCGCCAGGCGGCATATTGCGCACGGTCTTGCGCCCCATCTTCCCCGTTTCCcggaccct       >  4:243145/1‑267 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGGG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                                                                                                >  8:479188/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                                                                                                >  7:659032/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGc  >  2:183559/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGAATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCCGCGCAGGc  >  2:584008/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                    aTCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                                                                                                >  6:280899/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    aTCGGTCCCCTCGCCCCTCAGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAGTAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCACTTCGTGTTCCATaa                                                                                                                                >  6:81827/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     gcGGTCCCCTCGCCCCCATGGGGAA‑GGGGTTAGGGATAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                                                                                                               <  5:689923/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ctaacccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                                                                                             >  5:381634/5‑148 (MQ=255)
                                                                                                                                                                                                                                                                        taacccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                                                                                                            >  6:379286/4‑148 (MQ=255)
                                                                                                                                                                                                                                                                        taacccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                                                                                                                           >  7:722052/4‑149 (MQ=255)
                                                                                                                                                                                                                                                                         tCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                                                                                                                           >  5:908689/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          acccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                                                                                                          >  6:587681/2‑148 (MQ=255)
                                                                                                                                                                                                                                                                           cccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                                                                                                               <  8:821723/142‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                                                                                                       >  7:75240/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                            ccTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                                                                                                       <  8:75240/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                                                                                                                        >  5:907364/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             cTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa                                                                                                                      <  7:667718/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                                                                                                     >  5:387282/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                      
CGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGC  >  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: ‑