Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,216 A→G 6.0% D122G (GAT→GGT)  rluB → 23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,327,2160AG 6.0% 733.7 / 22.8 221D122G (GAT→GGT) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base A (103/104);  new base G (6/8);  total (109/112)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-02

Only 100 of 227 total aligned reads displayed.

GTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTGGTGTGCAGGTGAGCCGCCTGATCCGT  >  NC_000913/1326942‑1327487
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                 
           aaTTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTg                                                                                                                                                                                                                                                                         >  2:262531/1‑274 (MQ=255)
               gcaaaaCGCGGCGTTCGTGTTTAAGTTACCCCTGGGTGGAAAATCCGTATCGATGGTCAACTGATTGCGGTAGGTGGGTCCGCTGAGCAAATTTGTCGCGTTTTGGCCGATTACAATCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGGCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCaa                                                                                                                                                                                                                                                                   <  4:16391/271‑1 (MQ=255)
               gCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCaaa                                                                                                                                                                                                                                                                  <  4:262530/275‑1 (MQ=255)
                         ccGATCTTATTAAAGTTACCCCTGGTATAAAATACCGTATCTAGCGTCACCTGATTTCGGTACGGGAGTCCGCGGACATGATTGTTCGCGTTCTGGACGATTGCAAACCGGAAGGTGAGTTGTGTACCCGTAACGAGCCGGAAGGACGTCCAACGGTGTTTGCCCGTCTGCCTAAACTACGTGGCGCACGCTGGATTGCCGTGGGTAGCTTGGACGTTAATACCTCTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCTTTTAATg                                                                                                                                                                                                                                                        <  4:210035/274‑1 (MQ=255)
                                            cccTGGATTGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGTGTCGCCTGGACGGTAATACCTGTTGCATGCTTCTGTTCCCCACGGATGGTGAACTG‑‑GCAAACGGTTTAAGTCACCCAAGCCGTGACGTTg                                                                                                                                                                                                                                     >  4:148192/1‑275 (MQ=255)
                                                    tGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCCAGCCGTGAAGTGGAAAGTTaa                                                                                                                                                                                                                             >  4:116800/1‑275 (MQ=255)
                                                      aaaaTCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGACCGTGAata                                                                                                                                                                                                                           >  2:74200/1‑275 (MQ=255)
                                                                tcaaTGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                                                                                                                                                  <  2:319312/271‑1 (MQ=255)
                                                                tagATGGTTACCTGATATCGGTACGTGAGTCCGCTTAACAAATTTTGCGCGTTCTGGCCTAGTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGGCCAACGGTGTTTGACCGTCTGCCAAAACTGCGGGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                                                                                                                                                  <  2:240451/272‑1 (MQ=255)
                                                                                  cGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAg                                                                                                                                                                                                                                                <  4:490155/226‑1 (MQ=255)
                                                                                             tCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAGGGTGAGCCG‑‑GCAAACCGTTTAATGAACCAAAGCCGTGCAGTTGAACGTGAATATGCCGTGNNNGTGTGTGGGTAGGGTGACGCTGCGAAAcg                                                                                                                                                                                    >  3:147159/1‑274 (MQ=255)
                                                                                                gcggAACAAATTTGTAGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGAACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCg                                                                                                                                                                                 <  2:148193/272‑1 (MQ=255)
                                                                                                    aaCAAATTGGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                                                                                                                                             <  2:116801/275‑1 (MQ=255)
                                                                                                                                   aCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCggtggt                                                                                                                                                                                                                                                                               >  7:586097/1‑148 (MQ=255)
                                                                                                                                     cGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGa                                                                                                                                                                                                                                                                             <  8:586097/148‑1 (MQ=255)
                                                                                                                                        aaGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCTATGGTGa                                                                                                                                                                                                                                                                             <  7:676546/145‑1 (MQ=255)
                                                                                                                                        aaGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGa                                                                                                                                                                                                                                                                             >  8:676546/1‑145 (MQ=255)
                                                                                                                                              aGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCCGGTTGACGACGCGAAACTGCGTGATTTCAGTCGTGGCGTGCAGTTGGAAGGTGGTCCGGCAGCttt                                                                                                                                    >  2:347745/1‑274 (MQ=255)
                                                                                                                                                 tgtgTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                                                                                                                                                               <  5:820055/148‑1 (MQ=255)
                                                                                                                                                   tgtACCCGTGACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGtt                                                                                                                                                                                                                                                             <  7:113591/148‑1 (MQ=255)
                                                                                                                                                      aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTg‑‑g                                                                                                                                                                                                                                                                      >  5:682430/1‑136 (MQ=255)
                                                                                                                                                      aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                                                                                                                                                               <  7:517480/143‑1 (MQ=255)
                                                                                                                                                          gTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCa                                                                                                                                                                                                                                                      <  5:477395/148‑1 (MQ=255)
                                                                                                                                                          gTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGCACGTg                                                                                                                                                                                                                               >  4:438681/1‑171 (MQ=255)
                                                                                                                                                               gATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAg                                                                                                                  <  2:437053/275‑1 (MQ=255)
                                                                                                                                                                  ccGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGcc                                                                                                                                                                                                                                              >  5:905499/1‑148 (MQ=255)
                                                                                                                                                                   cGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTaa                                                                                                                                                                                                                                                          <  8:244527/135‑1 (MQ=255)
                                                                                                                                                                    ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                                                                                                                                                              >  8:413689/1‑130 (MQ=255)
                                                                                                                                                                        ggagggccAACGGTGTTTGACGGTTGGGAAAAAAGGCGGGGGGCACGATGGATGTCCGTGGGTCGTGTGGACGTTAAGCACTGTGGACTGCTGCTGTTCACGACGGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGGAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGa                                                                                                         <  4:227006/269‑1 (MQ=255)
                                                                                                                                                                         gACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAt                                                                                                                                                                                                                                                                                  >  5:513758/1‑107 (MQ=255)
                                                                                                                                                                           cGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTga                                                                                                                                                                                                <  2:217361/185‑1 (MQ=255)
                                                                                                                                                                             tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCg                                                                                                                                                                                                                                                               <  5:59703/120‑1 (MQ=255)
                                                                                                                                                                             tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaa                                                                                                                                                                                                                                   <  8:338756/148‑1 (MQ=255)
                                                                                                                                                                                  gggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCGAGCTGGATTGCCGTGGGTCGCGTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAAc                                                                                               <  1:148910/274‑1 (MQ=255)
                                                                                                                                                                                  cGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAc                                                                                                                                                                                                                                                                           >  7:739052/1‑105 (MQ=255)
                                                                                                                                                                                     tgtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                                                                                                                                                          >  5:429509/1‑133 (MQ=255)
                                                                                                                                                                                        ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                                                                                                                                                  <  7:85070/122‑1 (MQ=255)
                                                                                                                                                                                        ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                                                                                                                                                  >  8:85070/1‑122 (MQ=255)
                                                                                                                                                                                          gACCGTCTGCCAAAACTGCGTGGCGCCCGCTGGATTGCCGTGGGTCGCTTGGACTTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTGAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTa                                                                                       <  1:147159/275‑1 (MQ=255)
                                                                                                                                                                                                     aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAttt                                                                                                                                                                                                           >  7:554786/1‑148 (MQ=255)
                                                                                                                                                                                                       aaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGt                                                                                                                                                                                                                                >  7:315634/1‑125 (MQ=255)
                                                                                                                                                                                                           gCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                                                                                                                                                     >  5:946321/1‑148 (MQ=255)
                                                                                                                                                                                                            cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCCGTGGGAAGATGTTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCCAAGGGACACACCCGTGCTACACCGTGACTCTGACCg                                                                      >  4:292142/1‑274 (MQ=255)
                                                                                                                                                                                                             ttGGCGCCCGCTGGATTGCCGTGGGTCTCTTCGACGTTAATTCCTGTGGGCGGCTGCTGTTCACCACCGATTGTGACCTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGaa                                                                    <  1:55750/274‑1 (MQ=255)
                                                                                                                                                                                                             gTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGaa                                                                    >  2:95757/1‑275 (MQ=255)
                                                                                                                                                                                                             gTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGaa                                                                    >  2:189090/1‑275 (MQ=255)
                                                                                                                                                                                                              tGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                                                                                                                                                         >  8:1080442/1‑125 (MQ=255)
                                                                                                                                                                                                              tGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGGTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                                                                                                                                                         <  7:1080442/125‑1 (MQ=255)
                                                                                                                                                                                                                gcgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAtt                                                                                                                                                                                                            >  8:881497/1‑136 (MQ=255)
                                                                                                                                                                                                                     cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg                                                                                                                                                                                           >  7:405990/1‑148 (MQ=255)
                                                                                                                                                                                                                       cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCa                                                                                                                                                                                                                                                  <  8:710000/91‑1 (MQ=255)
                                                                                                                                                                                                                       cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGa                                                                                                                                                                                                                              >  8:225575/1‑111 (MQ=255)
                                                                                                                                                                                                                             tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                                                                                                                                                   <  8:120264/148‑1 (MQ=255)
                                                                                                                                                                                                                             tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                                                                                                                                                   >  8:81630/1‑148 (MQ=255)
                                                                                                                                                                                                                             tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                                                                                                                                                   <  8:521708/148‑1 (MQ=255)
                                                                                                                                                                                                                                cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg                                                                                                                                                                                         <  6:207600/139‑1 (MQ=255)
                                                                                                                                                                                                                                cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg                                                                                                                                                                                         >  5:207600/1‑139 (MQ=255)
                                                                                                                                                                                                                                cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt                                                                                                                                                                                >  6:946308/1‑148 (MQ=255)
                                                                                                                                                                                                                                  tGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCAAAACTGCGTGa                                                                                                                                                                              <  7:81630/148‑1 (MQ=255)
                                                                                                                                                                                                                                       cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTga                                                                                                                                                                                                >  7:252674/1‑125 (MQ=255)
                                                                                                                                                                                                                                         cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCtgt                                        >  2:462417/1‑275 (MQ=255)
                                                                                                                                                                                                                                               cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                                                                                                                                                 >  6:745318/1‑84 (MQ=255)
                                                                                                                                                                                                                                               cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg                                                                                                                                                                                   >  5:13188/1‑130 (MQ=255)
                                                                                                                                                                                                                                                 ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGt                                                                                                                                                                                                        <  7:1034222/107‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGt                                                                                                                                                                                                        >  8:1034222/1‑107 (MQ=255)
                                                                                                                                                                                                                                                 ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                                                                                                                                                     >  6:42123/1‑110 (MQ=255)
                                                                                                                                                                                                                                                         tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTa                                                                                                                                                                                                                                                           <  6:675221/48‑1 (MQ=255)
                                                                                                                                                                                                                                                         tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTgg                                                                                                                                                       >  6:800855/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCg                                                                                                                                                                     >  6:905347/1‑133 (MQ=255)
                                                                                                                                                                                                                                                            ggTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAg                                                                                                                                                    <  6:505280/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAggg                                                                                                     >  4:85732/1‑195 (MQ=255)
                                                                                                                                                                                                                                                               ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                                                                                                                                                              <  6:664493/39‑1 (MQ=255)
                                                                                                                                                                                                                                                               ctgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGt                                                                                                                                                                                                                                                              >  5:664493/1‑39 (MQ=255)
                                                                                                                                                                                                                                                                tgctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGACCAACGTGACTATGACCGAAGGCCGCAACCGGGAAGTTCGTCGTCTGTGGGAAGCGGTTGGTGTGCAGGTg                 >  2:536817/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 gctgctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACATGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGt                                                                                                                                                          <  7:430835/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gctgctgTTCACCACCGATGGTGAAATG‑‑GCAAACCGTTTAATGCACCCTAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                                                                                                                                                    <  6:271723/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCa                                                                                                                                                                                                      >  8:1135176/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                    gctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCg                                                                                                                                                                                 >  7:605967/1‑111 (MQ=255)
                                                                                                                                                                                                                                                                    gctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCg                                                                                                                                                                                 <  8:605967/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gctgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCg                                                                                                                                            >  5:418635/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      tgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                                                                                                                                                           <  6:224592/51‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tgTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                                                                                                                                                           >  5:224592/1‑51 (MQ=255)
                                                                                                                                                                                                                                                                        ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                                                                                                                                                         >  5:89308/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                        ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                                                                                                                                                                                                         <  6:89308/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ttCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                                                                                                                                             >  5:174077/1‑111 (MQ=255)
                                                                                                                                                                                                                                                                         tCACCATCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                                                                                                                                                  >  7:232952/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                         tCACCATCGGTGGTGAACAGCAGCAGACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGAt                                                                                                                                                  <  5:759337/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAat                                                                                                                                                                                                                            <  7:952208/63‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                                                                                                                                                                                                     >  6:565157/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                         tCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGc                                                                                                                                       >  5:1030414/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          caccatCGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCt                                                                                                                                      <  8:299229/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                           accaccGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                                                                                                                                                 <  5:377013/56‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccatcGGTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                                                                                                                                             <  6:437980/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccaccGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTGGTGTGCAGGTGAGCCGCCTGATc     >  2:387244/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                              accGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                                                                                                                                                          <  7:200121/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ccGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGCCTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAGGCGGTTGGTGTGCAGGTGAGCCGCCTGATCCgt  >  1:5170/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                 gATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTgg                                                                                                                                                                  <  6:134817/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTgg                                                                                                                                                                  <  8:510999/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaaa                                                                                                                               >  7:1126646/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  aTGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAAc                                                                                                                              >  6:11623/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                 
GTTGATGGCAAAATTGCTAAACTCGGCGATCGTGTTGAAGTTACCCCTGGATTGAAAATCCGTATCGATGGTCACCTGATTTCGGTACGTGAGTCCGCTGAACAAATTTGTCGCGTTCTGGCCTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTG‑‑GCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTGGTGTGCAGGTGAGCCGCCTGATCCGT  >  NC_000913/1326942‑1327487

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