Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,223,855 G→A 5.4% A10T (GCG→ACG)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,8550GA 5.4% 433.8 / 10.1 130A10T (GCG→ACG) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (60/63);  new base A (3/4);  total (63/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 138 total aligned reads displayed.

CAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGG  >  NC_000913/4223583‑4224128
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
cAACGGCGGGTTTTCTGCCGCGTTTCGCGTGAATTGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTAGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgc                                                                                                                                                                                                                                                                                  <  2:25686/274‑1 (MQ=255)
                    cgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAACTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCt                                                                                                                                                                                                                                                             >  4:287914/1‑275 (MQ=255)
                             ggAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGa                                                                                                                                                                                                                                                    <  2:257593/275‑1 (MQ=255)
                                                      gttttttCTGTATCTTGGAAATCATTTTCATTTTTGTTTGTAGCGAATACAATAGTTAGTGAACTGCTACACTGAGTTAATGTTGAACACATCTCATGTTTCGTGGTGGTTGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGGATTTTTATTGAATCTGGGTTGAGCCGGTCGGGAGGAAGTGTGTGCAGCAAAGTGGAACAACTGCGTGCGCAGTTGAATGAACGTAGTCTGGTGCTGGACGGCGGTATGGGCACCATGGTCCa                                                                                                                                                                                                                            <  4:274142/273‑1 (MQ=255)
                                                      ctTTTTGCTGTATCTGTGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTTAACAAATCTTATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCa                                                                                                                                                                                                                            <  1:57539/274‑1 (MQ=255)
                                                                                  cattttTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGGTGAGTTAATGTTGAACAAATCTCATGTTTCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTc                                                                                                                                                                                               <  4:328497/275‑1 (MQ=255)
                                                                                                        aaTAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTg                                                                                                                                                                         <  2:4668/275‑1 (MQ=255)
                                                                                                                ctgaactgaTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCg                                                                                                                                                                 <  4:134150/275‑1 (MQ=255)
                                                                                                                       gaTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc                                                                                                                                                                                                                                                        >  2:543114/1‑181 (MQ=255)
                                                                                                                        aTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAACTGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTTAACGCTTTGCCGACTGGCCATGCTCCCTc                                                                                                                                                            >  1:112106/1‑272 (MQ=255)
                                                                                                                           cGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGa                                                                                                                                                                                                                                                    <  4:543113/181‑1 (MQ=255)
                                                                                                                            gATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTGGAAGTTTTATTAAAGCTGGGGTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAgg                                                                                                                                                       <  2:114374/273‑1 (MQ=255)
                                                                                                                                  ttAATGTTGAACAAATTTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGACAGTATGGTCTGCTGAATTTTTATTACACCTGGGTTGAGCGTGTCGGGAGCAGGTGTGAGCAGCAAAGTGGAACAAGTGCGTGCGCTGTTAAATGAATGTATTCTGGTGCTGGACGACGGTATGGTCACCATGAGCCAGAGTTACCGATCGAACGAAACAGATTTCCGCGCGGCAACATTCGCAGACTGGCCACTCAACACCAAAGGAACAACGa                                                                                                                                              >  4:323672/1‑275 (MQ=255)
                                                                                                                                     aTGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttg                                                                                                                                                                                                                                                                           <  7:613873/148‑2 (MQ=255)
                                                                                                                                       gTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa                                                                                                                                                                                                                                                                         <  7:796494/148‑1 (MQ=255)
                                                                                                                                       gTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa                                                                                                                                                                                                                                                                         >  8:796494/1‑148 (MQ=255)
                                                                                                                                        ttGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                                                                                                                                                                                                        >  7:915064/1‑148 (MQ=255)
                                                                                                                                           aaCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                                                                                                                                                                                                                                                                     <  8:750707/148‑1 (MQ=255)
                                                                                                                                           aaCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAGGGCAACAACGACCTGCTgg                                                                                                                                      >  2:127843/1‑275 (MQ=255)
                                                                                                                                             cAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa                                                                                                                                                                                                                                                                          <  8:949747/141‑1 (MQ=255)
                                                                                                                                               aaTCTCATTTGGGTTGGTGGTAGCTTTTACCACAGGTGCGTTTTTGCCAGTATGGTTTGTTGAATTTGTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACt                                                                                                                                  <  4:16887/275‑1 (MQ=255)
                                                                                                                                                    cATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgcg                                                                                                                                                                                                                                                                                 >  5:934439/1‑127 (MQ=255)
                                                                                                                                                    cATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                                                                                                                                                                                                            >  5:2178/1‑148 (MQ=255)
                                                                                                                                                    cATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                                                                                                                                                                                                            <  8:137990/148‑1 (MQ=255)
                                                                                                                                                    cATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg                                                                                                                                                                                                                                                            <  8:915064/148‑1 (MQ=255)
                                                                                                                                                     aTGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTa                                                                                                                             >  2:480759/1‑274 (MQ=255)
                                                                                                                                                          gCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGCGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTCAACGCTTTGCAGACTGGCCATGCGAACTCAAAGGCCACCACGACCTGCCGGTACTCAGTAAGCCgg                                                                                                                       >  4:448077/1‑275 (MQ=255)
                                                                                                                                                               gtggtCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcgg                                                                                                                                                                                                                                                 <  6:2178/148‑1 (MQ=255)
                                                                                                                                                                tggtCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATTCGACCTCAAAGGCAACAACGACCTGCTGGTACCCAGTAAACCGGacttg                                                                                                                  >  2:205822/1‑270 (MQ=255)
                                                                                                                                                                 ggtCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGAt                                                                                                                <  1:44849/275‑1 (MQ=255)
                                                                                                                                                                    cGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcgg                                                                                                                                                                                                                                                 <  6:341822/143‑1 (MQ=255)
                                                                                                                                                                    cGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTa                                                                                                                                                                                                                                            >  7:873822/1‑148 (MQ=255)
                                                                                                                                                                    cGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTAATAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcgg                                                                                                                                                                                                                                                 >  5:341822/1‑143 (MQ=255)
                                                                                                                                                                        tttACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttg                                                                                                                                                                                                                                                                           >  7:267923/1‑112 (MQ=255)
                                                                                                                                                                              acaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa                                                                                                                                                                                                                                                                         >  7:673094/1‑109 (MQ=255)
                                                                                                                                                                              acaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCAt                                                                                                                                                                                                                                  >  7:613885/1‑148 (MQ=255)
                                                                                                                                                                                 gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGc                                                                                                                                                                                                                                       >  4:508862/1‑140 (MQ=255)
                                                                                                                                                                                 gATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTTAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTTCGGGTACTCAGAAAACCGGAAGTGATCGCCGCTATCCACAAc                                                                                                >  4:493588/1‑275 (MQ=255)
                                                                                                                                                                                   tGCGTGTATGCCAGTATGGTTTTTTGAAGTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCGTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGc                                                                                              <  2:13991/275‑1 (MQ=255)
                                                                                                                                                                                    gCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGa                                                                                                                                                                                                                                                                      <  6:1045206/106‑1 (MQ=255)
                                                                                                                                                                                    gCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGa                                                                                                                                                                                                                                                                      >  5:1045206/1‑106 (MQ=255)
                                                                                                                                                                                     cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCa                                                                                                                                                                                                                                   <  2:508863/140‑1 (MQ=255)
                                                                                                                                                                                        ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt                                                                                                                                                                                                                                                       >  8:366833/1‑117 (MQ=255)
                                                                                                                                                                                        ttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt                                                                                                                                                                                                                                                       <  7:366833/117‑1 (MQ=255)
                                                                                                                                                                                                gTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATc                                                                                                               >  2:155260/1‑245 (MQ=255)
                                                                                                                                                                                                 tATGGTTTGTTGAATTTTTATTGAATCTGGGTTGAGCGGGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGc                                                                                  <  2:416210/273‑1 (MQ=255)
                                                                                                                                                                                                 tATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTc                                                                                                                                                                                                                                                              <  6:156602/101‑1 (MQ=255)
                                                                                                                                                                                                 tATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCg                                                                                 <  4:480758/274‑1 (MQ=255)
                                                                                                                                                                                                 tATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACGGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGAAAACCGGAAGTGATCGCCGATATCAACAACGCCCACTTTGAAGCg                                                                                 >  4:56002/1‑274 (MQ=255)
                                                                                                                                                                                                 tATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGGGAACGCTTTGCCGACTGGCCATGCGACCCCAAAGGCCACAACGACCTGCGGGTACCCCGGAAACCGGAAGTACTCCCCGGCATCCAAAAACCCTACTTTCAACCgg                                                                                >  4:158762/1‑275 (MQ=255)
                                                                                                                                                                                                      tttGTTTAATTTTTTTTGAGTCTGGGTGTGGCGTGTAGGGAGCAAGGTTGGGCAGCAAAGTGGAAAAACTGCGTGCGAAGTTAAATGAACGTATGCTGGTGCTGNACGGCGGTATGGGCACCATGAGCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATcgccg                                                                                                           <  3:112106/243‑1 (MQ=255)
                                                                                                                                                                                                         tttgAATTTTTATTAAATTGGGGTTGGGCGGGTGGGGGGGAAGGGTCGACAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGGATGGGCACCATGTGCCAGGGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGat                                                                        <  4:513865/274‑1 (MQ=255)
                                                                                                                                                                                                         gttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGAt                                                                                                                                                                                                                               <  6:312336/124‑1 (MQ=255)
                                                                                                                                                                                                             aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATg                                                                                                                                                                                                                                 <  8:196865/118‑1 (MQ=255)
                                                                                                                                                                                                             aaTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATg                                                                                                                                                                                                                                 >  7:196865/1‑118 (MQ=255)
                                                                                                                                                                                                               tttttATTAAATCTGGGTTGAGCGTGTAGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGCGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGACTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACACCGACCTGCTGGTACTCAGTAAACTGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATaccatc                                                                  >  2:385592/1‑275 (MQ=255)
                                                                                                                                                                                                                    attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTAt                                                                                                                                                                                                                     >  8:453343/1‑123 (MQ=255)
                                                                                                                                                                                                                      taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATg                                                                                                                                                                                                                                 <  7:1013757/109‑1 (MQ=255)
                                                                                                                                                                                                                      taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                                                                                                                                                                                                          <  7:536913/116‑1 (MQ=255)
                                                                                                                                                                                                                      taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                                                                                                                                                                                                          >  8:536913/1‑116 (MQ=255)
                                                                                                                                                                                                                      taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCag                                                                                                                                                                                                                           <  8:244538/115‑1 (MQ=255)
                                                                                                                                                                                                                      taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCg                                                                                                                                                                                              <  6:350907/144‑1 (MQ=255)
                                                                                                                                                                                                                      taAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCg                                                                                                                                                                                              >  5:350907/1‑144 (MQ=255)
                                                                                                                                                                                                                       aaaTCTGGGGTGGGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGGAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCaa                                                          <  2:3060/275‑1 (MQ=255)
                                                                                                                                                                                                                         aTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt                                                                                                                                                                                                                                             >  6:1019026/1‑94 (MQ=255)
                                                                                                                                                                                                                         aTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt                                                                                                                                                                                                                                             <  5:1019026/94‑1 (MQ=255)
                                                                                                                                                                                                                           cTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa                                                                                                                                                                                                                                                                         <  7:976564/64‑1 (MQ=255)
                                                                                                                                                                                                                           cTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACg                                                                                                                                                                                     >  5:410101/1‑148 (MQ=255)
                                                                                                                                                                                                                             gggTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAtt                                                                                                                                                                                                  <  7:1115020/133‑1 (MQ=255)
                                                                                                                                                                                                                              ggTTGAGCGTGTCGGGGGAAAGGGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTc                                                   <  4:425646/275‑1 (MQ=255)
                                                                                                                                                                                                                              ggTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcgg                                                                                                                                                                                                                                                 >  8:1063748/1‑85 (MQ=255)
                                                                                                                                                                                                                              ggTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcgg                                                                                                                                                                                                                                                 <  7:1063748/85‑1 (MQ=255)
                                                                                                                                                                                                                                   aGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAc                                                                                                                                                                                                                                                   >  6:656653/1‑78 (MQ=255)
                                                                                                                                                                                                                                   aGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAc                                                                                                                                                                                                                                                   <  5:656653/78‑1 (MQ=255)
                                                                                                                                                                                                                                    gCGTGTCGGGAGCAAGGGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTcc                                             <  1:128696/275‑1 (MQ=255)
                                                                                                                                                                                                                                                cAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGaa                                                                                                                                                                                                             <  6:1017649/103‑1 (MQ=255)
                                                                                                                                                                                                                                                cAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAg                                                                                                                                                                                                        <  8:1009971/108‑1 (MQ=255)
                                                                                                                                                                                                                                                cAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAg                                                                                                                                                                                                        >  7:1009971/1‑108 (MQ=255)
                                                                                                                                                                                                                                                  aGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACCACGCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCACCACCTTCAACTCCACGACCATTGCGAt                               >  2:318092/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    tgtgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAtt                                                                                                                                                                                                  <  5:1001881/110‑1 (MQ=255)
                                                                                                                                                                                                                                                         gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCag                                                                                                                                                                                                                           >  8:64329/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                aGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTGATCGACTGAACGAGGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaacaa                                                                                                                                                 <  5:299495/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGa                                                                                                                                                                >  7:654065/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                  tGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGa                                                                                                                                                                <  8:654065/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTgg                                                                                                                                                                        >  5:115858/4‑119 (MQ=255)
                                                                                                                                                                                                                                                                     aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATg                                                                                                                                                                                                                                 >  6:424106/1‑62 (MQ=255)
                                                                                                                                                                                                                                                                     aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAg                                                                                                                                                                                                        >  8:1044751/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                     aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTc                                                                                                                                                            >  7:1009845/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                     aacaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTc                                                                                                                                                            <  8:1009845/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttaacTGCGCACGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCaa                                                                                                                                                          >  8:763770/3‑132 (MQ=255)
                                                                                                                                                                                                                                                                       caacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATcc                                                                                                                                                                                                                             <  5:696782/64‑1 (MQ=255)
                                                                                                                                                                                                                                                                       caacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATcc                                                                                                                                                                                                                             >  5:542434/1‑64 (MQ=255)
                                                                                                                                                                                                                                                                       caacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATcc                                                                                                                                                                                                                             <  6:542434/64‑1 (MQ=255)
                                                                                                                                                                                                                                                                       caacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATcc                                                                                                                                                                                                                             >  6:696782/1‑64 (MQ=255)
                                                                                                                                                                                                                                                                          cTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCtt                                                                                                                                                                                  >  5:897022/1‑104 (MQ=255)
                                                                                                                                                                                                                                                                          cTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCtt                                                                                                                                                                                  <  6:897022/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACTTTGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCANCAACGACCTGCTGGCACCCAGTAAACCGGAAGTGATCGCCGCTATCCACACCGCCAACTTTCAAGCGGGNGCGGACACCATCGAAACAAACACAGTAAACTCCACCACCATTGCGNANGGCGATTACCAGAAGGAATaccatgcgg  >  3:45375/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                               tgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCGGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATATCAACGAAACCAACACCTTCAACTTCACGACCATTGCTATGGCGGATTCCCAGATGGAATCCCTGTcag  >  2:535000/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                               tgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACCACGCCTACTTTGAAACGGGCACGGATATCATCAAAACCAACACCTTCCACCCCAACACCATTGCGCTGCCGGATTACCAGATGGAATCCCttcccg  >  4:225910/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                               tgcgCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTcgg  >  2:307804/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
CAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGG  >  NC_000913/4223583‑4224128

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