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.

TACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTC  >  NC_000913/4040397‑4040932
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                
                       ccTTACAACGCCTAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGtt‑gtttgt                                                                                                                                                                                                                                                <  2:316537/275‑1 (MQ=25)
                                gCCGAAGGGGTTTTGGCGGATTGAGAGAAGAGTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGTATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAGCTGTT‑GTTTGTCGGTGAACg                                                                                                                                                                                                                                       <  2:98327/275‑1 (MQ=25)
                                 ccGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAGGGCCCAGGCGAAAGAAAGGGCCTTGCGTTTTAACTGTT‑GTTTGTCGGTGAACg                                                                                                                                                                                                                                       >  4:263141/1‑274 (MQ=255)
                                 ccGAAGGTGTTTTGGCGGATTGAGAGAAGATTGTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACg                                                                                                                                                                                                                                       <  2:197340/274‑1 (MQ=25)
                                 ccGAAGGTGGTTTGGCGGAATGAGAGAAGAGTTTCAGCCTGATACCGATTAAATCATAAAGCAGAAGCGGTCTGATAAAACAGAATTTGCCTTGCGGAAGGAGCGCGGTGGTCCCACCTTACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACg                                                                                                                                                                                                                                       <  4:69312/274‑1 (MQ=11)
                                                               tgttCAGCATGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGTTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGAATAAATCCGCCGGGAg                                                                                                                                                                                                         <  2:82129/272‑1 (MQ=12)
                                                                                     aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTGTATCtg                                                                                                                                                                                                                                                         >  1:9110/1‑205 (MQ=17)
                                                                                         gAACGCAGGAACGGTCTGATAAAACAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgtt‑gt                                                                                                                                                                                                                                                    <  3:9110/205‑1 (MQ=17)
                                                                                               aGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGCGTCGGCCAAATCCGCCGGTAACGGATTGTAACGTTTCGAAGCAACGCGCCggg                                                                                                                                                                         >  4:329407/1‑273 (MQ=11)
                                                                                                      gTCTGAGAAAAAAGCATAGGCGCGGGGNCAGGAGGGGGGTGGTCCCACCTGTCNCCGTGACGAACTCAGAAGTGAAACGCCGGAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGNCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                                                                                                                                                                   <  3:57124/273‑1 (MQ=11)
                                                                                                                     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGa                                                                                                                                                                                                                                                                                 <  5:3951/149‑1 (MQ=34)
                                                                                                                     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTCTCTGTT‑GTTTGTCGGTGAACGCTCTCAGGAGTAGGACAACTCCGCCGGGAGCGGTTTTGAACGGTGCCACGCAACGGCCCGGGGGAGGGCGGGCCGGAcgcctgac                                                                                                                                                  >  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)
                                                                                                                                                 cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGGTTTATCtgtt‑gt                                                                                                                                                                                                                                                    >  7:1263877/1‑149 (MQ=21)
                                                                                                                                                   cACCTGACCCCATGCCGAACTCAGAAGTGATACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGtt‑gttt                                                                                                                                                                                                                                                  <  5:958268/149‑1 (MQ=25)
                                                                                                                                                            ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGaa                                                                                                                                                                                                                                         >  5:1107478/1‑149 (MQ=34)
                                                                                                                                                             cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAAc                                                                                                                                                                                                                                        <  8:622050/149‑1 (MQ=35)
                                                                                                                                                                gCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGGGGGTGGAGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAACGCCCTCCTTACg                                                                                                       >  4:266705/1‑275 (MQ=12)
                                                                                                                                                                 ccGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGctc                                                                                                                                                                                                                                    >  7:611668/1‑149 (MQ=35)
                                                                                                                                                                  cGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTgg                                                                                                                                                                                                                                                                          >  5:526701/1‑111 (MQ=38)
                                                                                                                                                                    aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTcc                                                                                                                                                                                                                                 >  5:725942/1‑149 (MQ=35)
                                                                                                                                                                       tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGa                                                                                                                                                                                                                              <  6:234801/149‑1 (MQ=34)
                                                                                                                                                                       tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGa                                                                                                                                                                                                                              >  8:745365/1‑149 (MQ=35)
                                                                                                                                                                        cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAg                                                                                                                                                                                                                             <  7:1342956/149‑1 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           <  7:94696/149‑1 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  5:686428/1‑149 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  6:568172/1‑149 (MQ=35)
                                                                                                                                                                          gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  5:300517/1‑149 (MQ=34)
                                                                                                                                                                            aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                                                                         <  6:295441/149‑1 (MQ=34)
                                                                                                                                                                            aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                                                                         >  7:696135/1‑149 (MQ=35)
                                                                                                                                                                             gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGa                                                                                                                                                                                                                        <  5:1016816/149‑1 (MQ=34)
                                                                                                                                                                             gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGa                                                                                                                                                                                                                        <  5:1363963/149‑1 (MQ=34)
                                                                                                                                                                             gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGa                                                                                                                                                                                                                        <  6:686428/149‑1 (MQ=35)
                                                                                                                                                                                  aCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAAt                                                                                                                                                                                                                   <  8:329553/149‑1 (MQ=35)
                                                                                                                                                                                       gTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                                                                                                                                                              <  8:351383/149‑1 (MQ=34)
                                                                                                                                                                                          gcgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGa                                                                                                                                                                                                                              <  8:1252577/130‑1 (MQ=37)
                                                                                                                                                                                           cgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGa                                                                                                                                                                                                          <  6:1333491/149‑1 (MQ=35)
                                                                                                                                                                                              cGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCg                                                                                                                                                                                                       <  7:289230/149‑1 (MQ=35)
                                                                                                                                                                                                  ggTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                                                                                                                                                   <  5:711185/149‑1 (MQ=34)
                                                                                                                                                                                                  ggTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                                                                                                                                                   >  7:341208/1‑149 (MQ=35)
                                                                                                                                                                                                      gtgtGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGaa                                                                                                                                                                                               >  6:502479/1‑149 (MQ=34)
                                                                                                                                                                                                       tgtgGGGTCTCCTCATGCGGGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                                                                              <  6:1158699/149‑1 (MQ=21)
                                                                                                                                                                                                         tgGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCATGCATCAAATTAAGCAGAACGCCAGCCTGACGGATGGGCTTTTTACATTTGCCCAGAAAAAAATCCCggct                                                                >  4:236006/1‑271 (MQ=255)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           <  7:1446356/113‑1 (MQ=37)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                                                                                                                                                           >  8:1446356/1‑113 (MQ=35)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGa                                                                                                                                                                                       >  5:77764/1‑149 (MQ=34)
                                                                                                                                                                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCACAAACTGCCAGGCATCAAATTACGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACg                                                          >  2:266174/1‑274 (MQ=255)
                                                                                                                                                                                                                tcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg                                                                                                                                                                                     >  6:248037/1‑149 (MQ=35)
                                                                                                                                                                                                                tcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAATGCCAGGCATACAATTAAGCAGAATGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAGGCCGGATGCGACGcc                                                        >  1:29257/1‑273 (MQ=255)
                                                                                                                                                                                                                  ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                                                                                                                                                   >  5:649417/1‑149 (MQ=34)
                                                                                                                                                                                                                   tcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCaa                                                                                                                                                                                  <  6:708006/149‑1 (MQ=32)
                                                                                                                                                                                                                        cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCTTCGCAACgttc                                                                                                                                                                             <  6:1352552/149‑4 (MQ=255)
                                                                                                                                                                                                                         gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                                                            >  5:407339/1‑149 (MQ=32)
                                                                                                                                                                                                                            agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGa                                                                                                                                                                         <  6:649417/149‑1 (MQ=32)
                                                                                                                                                                                                                              tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg                                                                                                                                                                       <  6:277695/149‑1 (MQ=32)
                                                                                                                                                                                                                              tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg                                                                                                                                                                       >  8:583193/1‑149 (MQ=32)
                                                                                                                                                                                                                               aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAggg                                                                                                                                                                      >  8:220877/1‑149 (MQ=32)
                                                                                                                                                                                                                                 ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                                                                                    >  7:447382/1‑149 (MQ=32)
                                                                                                                                                                                                                                 ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                                                                                    >  6:606152/1‑149 (MQ=32)
                                                                                                                                                                                                                                 ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                                                                                    <  7:220877/149‑1 (MQ=32)
                                                                                                                                                                                                                                   aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGc                                                                                                                                                                  <  6:109709/149‑1 (MQ=32)
                                                                                                                                                                                                                                        ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCACAGAAAAAACTGCCGGATTCAACGCTGGCGCGTCTTATCCAGCCTActct                                >  2:49633/1‑274 (MQ=255)
                                                                                                                                                                                                                                        ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAAGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCCGAGAAAAAAATGCCGAATGCCGAGCCGGCGCCGCCTATCCCGCCtcccat                                >  4:253418/1‑269 (MQ=14)
                                                                                                                                                                                                                                          aGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                                                           >  7:797034/1‑149 (MQ=32)
                                                                                                                                                                                                                                           ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                                                                                                                                                          <  6:1309782/149‑1 (MQ=32)
                                                                                                                                                                                                                                            gCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAc                                                                                                                                                         >  8:186247/1‑149 (MQ=32)
                                                                                                                                                                                                                                              ctCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCGCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAGTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTa                          <  1:29260/273‑1 (MQ=255)
                                                                                                                                                                                                                                               tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcc                                                                                                                                                      >  7:1245063/1‑149 (MQ=32)
                                                                                                                                                                                                                                                   aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                                                                              <  8:113009/105‑1 (MQ=33)
                                                                                                                                                                                                                                                   aTAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccgcacactaccatcggcgctacggcgtttcacttctgagttcggcaggggggcaggtgggaccaccg                                                                                                                                                   >  5:1086843‑M1/1‑44 (MQ=255)
                                                                                                                                                                                                                                                     aaaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt                                                                                                                                                >  6:833287/1‑149 (MQ=32)
                                                                                                                                                                                                                                                      aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                                                                              <  5:274536/102‑1 (MQ=32)
                                                                                                                                                                                                                                                      aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATa                                                                                                                                               >  6:535432/1‑149 (MQ=32)
                                                                                                                                                                                                                                                       aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa                                                                                                                                              >  5:142902/1‑149 (MQ=32)
                                                                                                                                                                                                                                                       aaCGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgct                                                                                                                                               >  6:750440‑M1/1‑40 (MQ=255)
                                                                                                                                                                                                                                                        aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                                             <  7:745365/149‑1 (MQ=32)
                                                                                                                                                                                                                                                        aCGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATctgatgcctggcagttccctactctcgcatggggagaccccac                                                                                                                                                                                                                 >  7:708761‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                                                          gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctact                                                                                                                                            >  8:846620‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                                                                          gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                                                           >  5:1234077/1‑149 (MQ=32)
                                                                                                                                                                                                                                                            aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgc                                                                                                                                          >  6:1445534‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                                                            aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                                                           >  4:18297/1‑131 (MQ=32)
                                                                                                                                                                                                                                                             aGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgcc                                                                                                                                         >  5:818280‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                                                             aGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgcc                                                                                                                                         >  6:446108‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                                                             aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                                                        >  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)
                                                                                                                                                                                                                                                                cTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgc                                                                                                                                                       <  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)
                                                                                                                                                                                                                                                                 tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCtctc                                                                                                                                                                                                                                  <  5:1313079/55‑1 (MQ=35)
                                                                                                                                                                                                                                                                     ctggAAGCTTGGGCTTTTCGTTTTAT‑TGTTGGTTGGCCGGTGAACCCTCTCCTGAGAAGGACAAACCCGCCGGGGGCGGTTTTGAACGTTGCGAAGCAACGGCCCGGAGGTGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccactc  <  2:253418/273‑1 (MQ=21)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                
TACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTT‑GTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTC  >  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: ‑