Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,040,661 G→A 5.6% intergenic (+25/+245) rrfA → / ← mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,040,6610GA 5.6% 560.0 / 16.0 187intergenic (+25/+245)rrfA/mobB5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (86/90);  new base A (8/3);  total (94/93)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 189 total aligned reads displayed.

TACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTC  >  NC_000913/4040397‑4040932
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                
                       ccTTACAACGCCTAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttgt                                                                                                                                                                                                                                                <  2:316537/275‑1 (MQ=25)
                                gCCGAAGGGGTTTTGGCGGATTGAGAGAAGAGTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGTATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAGCTGTTGTTTGTCGGTGAACg                                                                                                                                                                                                                                       <  2:98327/275‑1 (MQ=25)
                                 ccGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAGGGCCCAGGCGAAAGAAAGGGCCTTGCGTTTTAACTGTTGTTTGTCGGTGAACg                                                                                                                                                                                                                                       >  4:263141/1‑274 (MQ=255)
                                 ccGAAGGTGTTTTGGCGGATTGAGAGAAGATTGTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                                                                                                                                                       <  2:197340/274‑1 (MQ=25)
                                 ccGAAGGTGGTTTGGCGGAATGAGAGAAGAGTTTCAGCCTGATACCGATTAAATCATAAAGCAGAAGCGGTCTGATAAAACAGAATTTGCCTTGCGGAAGGAGCGCGGTGGTCCCACCTTACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                                                                                                                                                       <  4:69312/274‑1 (MQ=11)
                                                               tgttCAGCATGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGTTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGAATAAATCCGCCGGGAg                                                                                                                                                                                                         <  2:82129/272‑1 (MQ=12)
                                                                                     aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTGTATCtg                                                                                                                                                                                                                                                         >  1:9110/1‑205 (MQ=17)
                                                                                         gAACGCAGGAACGGTCTGATAAAACAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgttgt                                                                                                                                                                                                                                                    <  3:9110/205‑1 (MQ=17)
                                                                                               aGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGCGTCGGCCAAATCCGCCGGTAACGGATTGTAACGTTTCGAAGCAACGCGCCggg                                                                                                                                                                         >  4:329407/1‑273 (MQ=11)
                                                                                                      gTCTGAGAAAAAAGCATAGGCGCGGGGNCAGGAGGGGGGTGGTCCCACCTGTCNCCGTGACGAACTCAGAAGTGAAACGCCGGAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGNCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                                                                                                                                                                   <  3:57124/273‑1 (MQ=11)
                                                                                                                     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGa                                                                                                                                                                                                                                                                                 <  5:3951/149‑1 (MQ=34)
                                                                                                                     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTCTCTGTTGTTTGTCGGTGAACGCTCTCAGGAGTAGGACAACTCCGCCGGGAGCGGTTTTGAACGGTGCCACGCAACGGCCCGGGGGAGGGCGGGCCGGAcgcctgac                                                                                                                                                  >  4:306444/1‑271 (MQ=11)
                                                                                                                      tttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGaa                                                                                                                                                                                                                                                                                <  6:980201/149‑1 (MQ=34)
                                                                                                                      tttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGaa                                                                                                                                                                                                                                                                                <  8:1307696/149‑1 (MQ=35)
                                                                                                                            tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTgg                                                                                                                                                                                                                                                                          <  8:516681/149‑1 (MQ=35)
                                                                                                                             ggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTggg                                                                                                                                                                                                                                                                         <  6:927712/149‑1 (MQ=35)
                                                                                                                                      aGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGtt                                                                                                                                                                                                                                                                <  8:1406685/149‑1 (MQ=35)
                                                                                                                                         gcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTA                                                                                                                                                                                                                                                             >  7:1374355‑M1/1‑149 (MQ=255)
                                                                                                                                         gcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTA                                                                                                                                                                                                                                                             >  5:903966‑M1/1‑149 (MQ=255)
                                                                                                                                          cgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT                                                                                                                                                                                                                                                            <  6:1405849‑M1/149‑1 (MQ=255)
                                                                                                                                                 cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGGTTTATCtgttgt                                                                                                                                                                                                                                                    >  7:1263877/1‑149 (MQ=21)
                                                                                                                                                   cACCTGACCCCATGCCGAACTCAGAAGTGATACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgttt                                                                                                                                                                                                                                                  <  5:958268/149‑1 (MQ=25)
                                                                                                                                                            ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGaa                                                                                                                                                                                                                                         >  5:1107478/1‑149 (MQ=34)
                                                                                                                                                             cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAc                                                                                                                                                                                                                                        <  8:622050/149‑1 (MQ=35)
                                                                                                                                                                gCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGGGGGTGGAGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAACGCCCTCCTTACg                                                                                                       >  4:266705/1‑275 (MQ=12)
                                                                                                                                                                 ccGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGctc                                                                                                                                                                                                                                    >  7:611668/1‑149 (MQ=35)
                                                                                                                                                                  cGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTgg                                                                                                                                                                                                                                                                          >  5:526701/1‑111 (MQ=38)
                                                                                                                                                                    aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTcc                                                                                                                                                                                                                                 >  5:725942/1‑149 (MQ=35)
                                                                                                                                                                       tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGa                                                                                                                                                                                                                              <  6:234801/149‑1 (MQ=34)
                                                                                                                                                                       tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGa                                                                                                                                                                                                                              >  8:745365/1‑149 (MQ=35)
                                                                                                                                                                        cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                                                                                                                                                             <  7:1342956/149‑1 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           <  7:94696/149‑1 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  5:686428/1‑149 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  6:568172/1‑149 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  5:300517/1‑149 (MQ=34)
                                                                                                                                                                            aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                                                                         <  6:295441/149‑1 (MQ=34)
                                                                                                                                                                            aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                                                                         >  7:696135/1‑149 (MQ=35)
                                                                                                                                                                             gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGa                                                                                                                                                                                                                        <  5:1016816/149‑1 (MQ=34)
                                                                                                                                                                             gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGa                                                                                                                                                                                                                        <  5:1363963/149‑1 (MQ=34)
                                                                                                                                                                             gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGa                                                                                                                                                                                                                        <  6:686428/149‑1 (MQ=35)
                                                                                                                                                                                  aCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAAt                                                                                                                                                                                                                   <  8:329553/149‑1 (MQ=35)
                                                                                                                                                                                       gTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                                                                                                                                                              <  8:351383/149‑1 (MQ=34)
                                                                                                                                                                                          gcgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGa                                                                                                                                                                                                                              <  8:1252577/130‑1 (MQ=37)
                                                                                                                                                                                           cgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGa                                                                                                                                                                                                          <  6:1333491/149‑1 (MQ=35)
                                                                                                                                                                                              cGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCg                                                                                                                                                                                                       <  7:289230/149‑1 (MQ=35)
                                                                                                                                                                                                  ggTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                                                                                                                                                   <  5:711185/149‑1 (MQ=34)
                                                                                                                                                                                                  ggTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                                                                                                                                                   >  7:341208/1‑149 (MQ=35)
                                                                                                                                                                                                      gtgtGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGaa                                                                                                                                                                                               >  6:502479/1‑149 (MQ=34)
                                                                                                                                                                                                       tgtgGGGTCTCCTCATGCGGGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                                                                              <  6:1158699/149‑1 (MQ=21)
                                                                                                                                                                                                         tgGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCATGCATCAAATTAAGCAGAACGCCAGCCTGACGGATGGGCTTTTTACATTTGCCCAGAAAAAAATCCCggct                                                                >  4:236006/1‑271 (MQ=255)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           <  7:1446356/113‑1 (MQ=37)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  8:1446356/1‑113 (MQ=35)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGa                                                                                                                                                                                       >  5:77764/1‑149 (MQ=34)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCACAAACTGCCAGGCATCAAATTACGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACg                                                          >  2:266174/1‑274 (MQ=255)
                                                                                                                                                                                                                tcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg                                                                                                                                                                                     >  6:248037/1‑149 (MQ=35)
                                                                                                                                                                                                                tcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAATGCCAGGCATACAATTAAGCAGAATGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAGGCCGGATGCGACGcc                                                        >  1:29257/1‑273 (MQ=255)
                                                                                                                                                                                                                  ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                                                                                                                                                   >  5:649417/1‑149 (MQ=34)
                                                                                                                                                                                                                   tcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCaa                                                                                                                                                                                  <  6:708006/149‑1 (MQ=32)
                                                                                                                                                                                                                        cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCTTCGCAACgttc                                                                                                                                                                             <  6:1352552/149‑4 (MQ=255)
                                                                                                                                                                                                                         gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                                                            >  5:407339/1‑149 (MQ=32)
                                                                                                                                                                                                                            agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGa                                                                                                                                                                         <  6:649417/149‑1 (MQ=32)
                                                                                                                                                                                                                              tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg                                                                                                                                                                       <  6:277695/149‑1 (MQ=32)
                                                                                                                                                                                                                              tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg                                                                                                                                                                       >  8:583193/1‑149 (MQ=32)
                                                                                                                                                                                                                               aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAggg                                                                                                                                                                      >  8:220877/1‑149 (MQ=32)
                                                                                                                                                                                                                                 ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                                                                                    >  7:447382/1‑149 (MQ=32)
                                                                                                                                                                                                                                 ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                                                                                    >  6:606152/1‑149 (MQ=32)
                                                                                                                                                                                                                                 ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                                                                                    <  7:220877/149‑1 (MQ=32)
                                                                                                                                                                                                                                   aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGc                                                                                                                                                                  <  6:109709/149‑1 (MQ=32)
                                                                                                                                                                                                                                        ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCACAGAAAAAACTGCCGGATTCAACGCTGGCGCGTCTTATCCAGCCTActct                                >  2:49633/1‑274 (MQ=255)
                                                                                                                                                                                                                                        ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAAGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCCGAGAAAAAAATGCCGAATGCCGAGCCGGCGCCGCCTATCCCGCCtcccat                                >  4:253418/1‑269 (MQ=14)
                                                                                                                                                                                                                                          aGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                                                           >  7:797034/1‑149 (MQ=32)
                                                                                                                                                                                                                                           ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                                                                                                                                                          <  6:1309782/149‑1 (MQ=32)
                                                                                                                                                                                                                                            gCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAc                                                                                                                                                         >  8:186247/1‑149 (MQ=32)
                                                                                                                                                                                                                                              ctCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCGCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAGTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTa                          <  1:29260/273‑1 (MQ=255)
                                                                                                                                                                                                                                               tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcc                                                                                                                                                      >  7:1245063/1‑149 (MQ=32)
                                                                                                                                                                                                                                                   aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                                                                              <  8:113009/105‑1 (MQ=33)
                                                                                                                                                                                                                                                   aTAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccgcacactaccatcggcgctacggcgtttcacttctgagttcggcaggggggcaggtgggaccaccg                                                                                                                                                   >  5:1086843‑M1/1‑44 (MQ=255)
                                                                                                                                                                                                                                                     aaaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt                                                                                                                                                >  6:833287/1‑149 (MQ=32)
                                                                                                                                                                                                                                                      aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                                                                              <  5:274536/102‑1 (MQ=32)
                                                                                                                                                                                                                                                      aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATa                                                                                                                                               >  6:535432/1‑149 (MQ=32)
                                                                                                                                                                                                                                                       aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa                                                                                                                                              >  5:142902/1‑149 (MQ=32)
                                                                                                                                                                                                                                                       aaCGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgct                                                                                                                                               >  6:750440‑M1/1‑40 (MQ=255)
                                                                                                                                                                                                                                                        aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                                             <  7:745365/149‑1 (MQ=32)
                                                                                                                                                                                                                                                        aCGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATctgatgcctggcagttccctactctcgcatggggagaccccac                                                                                                                                                                                                                 >  7:708761‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                                                          gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctact                                                                                                                                            >  8:846620‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                                                                          gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                                                           >  5:1234077/1‑149 (MQ=32)
                                                                                                                                                                                                                                                            aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgc                                                                                                                                          >  6:1445534‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                                                            aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                                                           >  4:18297/1‑131 (MQ=32)
                                                                                                                                                                                                                                                             aGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgcc                                                                                                                                         >  5:818280‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                                                             aGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgcc                                                                                                                                         >  6:446108‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                                                             aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                                                        >  6:1012345/1‑149 (MQ=32)
                                                                                                                                                                                                                                                              ggCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccg                                                                                                                                        >  5:1167571‑M1/1‑33 (MQ=255)
                                                                                                                                                                                                                                                              ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATctgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttc                                                                                                                                                                                        <  7:1291704‑M1/101‑69 (MQ=255)
                                                                                                                                                                                                                                                               gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccgc                                                                                                                                       >  5:1363542‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                                                cTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgc                                                                                                                                                       <  2:18297/131‑1 (MQ=32)
                                                                                                                                                                                                                                                                cTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtt                                                                                                                                                                                          <  5:337263‑M1/97‑67 (MQ=255)
                                                                                                                                                                                                                                                                cTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtt                                                                                                                                                                                          >  6:337263‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                                                                cTCAGTCGAAAAACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccgcc                                                                                                                                      >  7:1146878‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                                                                 tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCtctc                                                                                                                                                                                                                                  <  5:1313079/55‑1 (MQ=35)
                                                                                                                                                                                                                                                                     ctggAAGCTTGGGCTTTTCGTTTTATTGTTGGTTGGCCGGTGAACCCTCTCCTGAGAAGGACAAACCCGCCGGGGGCGGTTTTGAACGTTGCGAAGCAACGGCCCGGAGGTGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccactc  <  2:253418/273‑1 (MQ=21)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                
TACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTC  >  NC_000913/4040397‑4040932

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