Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,061,818 G→C 6.3% W323C (TGG→TGC) ‡ scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,061,8180GC 6.3% 648.8 / 26.3 188W323C (TGG→TGC) ‡scpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (80/96);  new base C (3/9);  total (83/105)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

Only 100 of 189 total aligned reads displayed.

ACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGC  >  NC_000913/3061547‑3062081
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
acaTGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAg                                                                                                                                                                                                                                                                       <  4:545807/274‑1 (MQ=255)
               gTGCCAATTGCGTGGAGCAGGTAGAATTTACGCTTGCTGATGGGATTGAGTATATAAAAGCAGCAATCTCTGCCGGATTGAAAATTGATGACTTCGCTCCTAGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                                                                                                                                       <  2:629135/275‑1 (MQ=255)
               gTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTGTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                                                                                                                                       <  2:234441/275‑1 (MQ=255)
                                          ttACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcacca                                                                                                                                                                                                                              >  2:193400/1‑273 (MQ=255)
                                            aCGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCACCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAGGCGGTCAGTGGATTTGGCGCACAGGCCCCGAAATCACTGGCGTTGCGTCCCCACTGCCAGAACTCAGGCTGGCGACTGCCTGAACAGGGTCCGTAAAACGACGTTATCCCAATCACCAt                                                                                                                                                                                                                          >  4:249718/1‑275 (MQ=255)
                                                 cgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      >  4:246365/1‑274 (MQ=255)
                                                                               aaTCTCTGCCGGCCTGAAAATTGTTGCCTTCGCTCCTGGCCTGTCGTTCTTCTTCGGCATCGGAATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGt                                                                                                                                                                                       <  4:193400/275‑1 (MQ=255)
                                                                                     tGCCGGACTGAAAATTGAGTACTTCGCTCCTCGCCGGTCGTTCTTCTTCGGCATCGGCATGGATCGGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTGGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGc                                                                                                                                                                                 <  1:143253/275‑1 (MQ=255)
                                                                                          gACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCCGTCACTGCATAcc                                                                                                                                                                            >  2:101240/1‑275 (MQ=255)
                                                                                                      attccTTCGCCCCTCGCATGTCGTTCTTCTTCCGCTTCGTCATGGATCTGTTCAGGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGCGGATTTGGCGAACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACCACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGAc                                                                                                                                                                <  1:24958/271‑1 (MQ=255)
                                                                                                           ttCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCCGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTAGAACAACAGTAACCGCCACACCATTGAAGCGCCTGCTGCCACGCTGGGCGTTACCCAGTCACTCCATACCACACCCTTTACCGAAgc                                                                                                                                                           >  4:489844/1‑275 (MQ=255)
                                                                                                                  ctcGCCTGTCGTTTTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAATTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGt                                                                                                                                                    <  2:56241/275‑1 (MQ=255)
                                                                                                                            gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGa                                                                                                                                                                                                                                                                        <  6:573466/149‑1 (MQ=255)
                                                                                                                            gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGa                                                                                                                                                                                                                                                                        <  7:262068/149‑1 (MQ=255)
                                                                                                                                 tcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAACCGCTTGGTTGTCCTACCGAttt                                                                                                                                      >  4:417021/1‑274 (MQ=255)
                                                                                                                                   ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgac                                                                                                                                                                                                                                                                 <  7:920672/149‑1 (MQ=255)
                                                                                                                                   ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                                                                                                                                                 <  6:266349/149‑2 (MQ=255)
                                                                                                                                   ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                                                                                                                                                 <  5:605768/149‑2 (MQ=255)
                                                                                                                                   ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                                                                                                                                                 <  7:119763/149‑2 (MQ=255)
                                                                                                                                    tcggcatcggcaTGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCTAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgagg                                                                                                                                                                                                                                                                <  6:943851/149‑3 (MQ=255)
                                                                                                                                        catcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAAc                                                                                                                                                                                                                                                            <  5:280652/149‑1 (MQ=255)
                                                                                                                                            ggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGAACCCGAAATCACTGGCGATGCGTACCCACTGCCAGACCTCGGGCTGGAGCCTGACTGAACAGGa                                                                                                                                                                                                                                                        >  8:643937/1‑149 (MQ=255)
                                                                                                                                             gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                                                                                                                                       >  7:831902/1‑149 (MQ=255)
                                                                                                                                               atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATcc                                                                                                                                                                                                                                                     >  7:218629/1‑149 (MQ=255)
                                                                                                                                               atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATcc                                                                                                                                                                                                                                                     <  7:547801/149‑1 (MQ=255)
                                                                                                                                               atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                                                                                                                                                 >  8:106095/1‑136 (MQ=255)
                                                                                                                                                tGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCg                                                                                                                                                                                                                                                    >  5:929220/1‑149 (MQ=255)
                                                                                                                                                tGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCg                                                                                                                                                                                                                                                    <  8:273426/149‑1 (MQ=255)
                                                                                                                                                  gATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGta                                                                                                                                                                                                                                                  >  8:210858/1‑149 (MQ=255)
                                                                                                                                                   aTCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtat                                                                                                                                                                                                                                                 >  7:481916/1‑149 (MQ=255)
                                                                                                                                                       gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCCCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                                                                                                                                             >  8:13939/1‑149 (MQ=255)
                                                                                                                                                        tttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGTCCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCa                                                                                                                                                                                         >  4:45488/1‑200 (MQ=255)
                                                                                                                                                         ttATGAACGTCTGCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACAcc                                                                                                             <  2:593067/275‑1 (MQ=255)
                                                                                                                                                               aCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATc                                                                                                                                                                                                                                     >  7:690034/1‑149 (MQ=255)
                                                                                                                                                               aCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGAACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatc                                                                                                       <  2:249718/275‑1 (MQ=255)
                                                                                                                                                                cGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATcc                                                                                                                                                                                                                                    >  6:837635/1‑149 (MQ=255)
                                                                                                                                                                cGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCAt                                                                                                                                                                               >  2:517080/1‑202 (MQ=255)
                                                                                                                                                                 ggcGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAGTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGGCCTCAGGCTGGAGCCTGACTGAACAGAATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatca                                                                                                      <  1:190615/272‑1 (MQ=255)
                                                                                                                                                                  tCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTTCCAGACCTCAGCCTGGAGCCTGACTGAACAGGATCCGtat                                                                                                                                                                                                                                                 >  1:62483/1‑134 (MQ=255)
                                                                                                                                                                   cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGca                                                                                                                                                                                                                                 <  5:712398/149‑1 (MQ=255)
                                                                                                                                                                     ccATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcacc                                                                                                                                                                                                                               <  7:577896/149‑1 (MQ=255)
                                                                                                                                                                     ccATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatcatc                                                                                                 <  2:417021/275‑1 (MQ=255)
                                                                                                                                                                      cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCAc                                                                                                                                                                                    <  4:57946/191‑1 (MQ=255)
                                                                                                                                                                      cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTAAGTGGATTTGGAGGACAGGACCCGAAATCACTGGCGCAGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGAACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCa                                                                                                           <  2:230829/264‑1 (MQ=255)
                                                                                                                                                                           tgcgtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAtt                                                                                                                                                                                                                         <  5:559433/149‑1 (MQ=255)
                                                                                                                                                                              gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      <  8:791838/149‑1 (MQ=255)
                                                                                                                                                                              gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      <  7:202136/149‑1 (MQ=255)
                                                                                                                                                                               tgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACa                                                                                                                                                                                                                                                           >  2:589167/1‑111 (MQ=255)
                                                                                                                                                                                  ggCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAgcgc                                                                                                                                                                                                                  >  5:399204/1‑149 (MQ=255)
                                                                                                                                                                                    cACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGc                                                                                                                                                                                                      >  2:342931/1‑159 (MQ=255)
                                                                                                                                                                                       gTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGaa                                                                                                                                                                                                                                                             >  3:56841/1‑101 (MQ=255)
                                                                                                                                                                                       gTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctggct                                                                                                                                                                                                             >  8:1040257/1‑149 (MQ=255)
                                                                                                                                                                                        ttatttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTggg                                                                                                                                                                                                  <  4:342931/159‑1 (MQ=255)
                                                                                                                                                                                                ggAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTg                                                                                                                                                                                                    >  6:526574/1‑149 (MQ=255)
                                                                                                                                                                                                   gCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGc                                                                                                                                                                                                 >  6:201415/1‑149 (MQ=255)
                                                                                                                                                                                                     gAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCgg                                                                                                                                                                                               <  5:526574/149‑1 (MQ=255)
                                                                                                                                                                                                     gAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCgg                                                                                                                                                                                               <  7:402631/149‑1 (MQ=255)
                                                                                                                                                                                                      aaGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACt                                                                                                                                                                                           >  4:163824/1‑152 (MQ=255)
                                                                                                                                                                                                          ggTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGt                                                                                                                                                                                       <  2:163824/152‑1 (MQ=255)
                                                                                                                                                                                                           gTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCa                                                                                                                                                                                         <  5:165449/149‑1 (MQ=255)
                                                                                                                                                                                                               gTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCa                                                                                                                                                                                     >  5:99021/1‑149 (MQ=255)
                                                                                                                                                                                                                tGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAg                                                                                                                                                                                        <  8:359848/145‑1 (MQ=255)
                                                                                                                                                                                                                gggttttGGCGAACAGAACCCGAAAGCACTGGCGCTGCGTCCCCCCTGCCATACCCAAGGCTGGGGCCTGACTGAACAGGACCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGcc                                                       <  4:321646/270‑1 (MQ=255)
                                                                                                                                                                                                                 ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                                                                                                                                   <  5:837635/149‑1 (MQ=255)
                                                                                                                                                                                                                 ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                                                                                                                                   >  8:466789/1‑149 (MQ=255)
                                                                                                                                                                                                                 ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                                                                                                                                   <  7:5004/149‑1 (MQ=255)
                                                                                                                                                                                                                   aTTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGa                                                    <  2:659804/274‑1 (MQ=255)
                                                                                                                                                                                                                       ggCCCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTggg                                                                                                                                                                                                  <  6:854948/128‑1 (MQ=255)
                                                                                                                                                                                                                        gcgcCCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTa                                               <  2:32344/274‑1 (MQ=255)
                                                                                                                                                                                                                          gcACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAg                                                                                    <  4:604507/235‑1 (MQ=255)
                                                                                                                                                                                                                             caGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACTCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCTATCAACTGGCTGGATCCTCTTACa                                          >  1:208802/1‑274 (MQ=255)
                                                                                                                                                                                                                              aGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgct                                                                                                                                                                                                                 >  8:517970/1‑106 (MQ=255)
                                                                                                                                                                                                                                gACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCtt                                                                                                                                                                    <  8:23666/149‑1 (MQ=255)
                                                                                                                                                                                                                                 aCCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCGGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCttt                                                                                                                                                                   <  7:643937/149‑1 (MQ=255)
                                                                                                                                                                                                                                 aCCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCttt                                                                                                                                                                   <  7:442425/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTg                                                                                                                                                                  <  8:739593/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTAAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                                                                                                                                                 >  5:661687/1‑149 (MQ=255)
                                                                                                                                                                                                                                      aaaTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTa                                                                                                                                             >  1:164483/1‑166 (MQ=255)
                                                                                                                                                                                                                                                   tGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTg                                                                                                                                                 <  7:966088/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGc                                                                                                                                                <  8:546746/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        aCCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTAc                                                                                                                                            >  7:1041826/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                                                                                                                                          >  7:1036357/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                                                                                                                                          >  7:229425/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           cACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGa                                                                                                                                         <  8:837233/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAt                                                                                                                                        >  8:594012/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              tGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                                                                                                                                      >  7:996674/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGt                                                                                                                                                                                                                                                   >  8:181098/1‑38 (MQ=255)
                                                                                                                                                                                                                                                                 cAGACCTCAGGCTGGTGCCTGACTGATCAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGTCGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCTTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCCAATCGTCAAACAAGCCag     <  2:78478/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ttcagtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGaa                                                                                                                                                             >  6:660937/6‑122 (MQ=255)
                                                                                                                                                                                                                                                                  aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa                                                                                                                                  <  8:1041826/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa                                                                                                                                  <  8:416983/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa                                                                                                                                  >  5:1072701/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     ccTCAGGCTGGAGCCTGACTGAAAAGGATCCGTATAACAACGTTATCCGGACCACCATTGAAGCGCTGGCTCCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCa                                                                                                                               >  7:217188/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     acTCAGGCTGGAGAAGATCTGGAGAGGATCCGTGTAGCAACGTTATGCGCGCCGCAAAGGAAGGGCGGGTCGCGGGGCTGGGCGGTAGTCAGACACTGCATCCCAACGCGTGTGACGAAGCGAGTGGTTGGACTACCGATGTCTCAGCACGCATTGCCCGAAACACCCAAATCATCATCCAGGAAGAGTCACAATTCTGTCGCACCGTCGAACCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAGATCGTCAAACAAGCCAGAGc  <  4:101240/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tCAGTCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTAACTCCATACCAACGCCTTTGACGAATCGCTTGGTTTGCCTACCGATTTCTCAGCACg                                                                                                                             >  8:1036343/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACg                                                                                                                             >  7:245440/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        cAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGNGTTGCCTACCGATTTCTCAGCACNCATTNcc                                                                                                                      >  4:657229/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                          ggCTGGAGCCTGACTGAACAGGATCCGTATAACAACg                                                                                                                                                                                                                                          <  2:291111/37‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCAtt                                                                                                                         >  7:537042/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               gAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGccc                                                                                                                     >  6:628623/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
ACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGC  >  NC_000913/3061547‑3062081

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