Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,707 C→T 7.6% D106N (GAT→AAT)  yneE ← bestrophin family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,608,7070CT 7.6% 424.3 / 13.7 140D106N (GAT→AAT) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (58/70);  new base T (6/5);  total (64/76)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

Only 100 of 141 total aligned reads displayed.

AGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAG  >  NC_000913/1608439‑1608969
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
             tatCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTagc                                                                                                                                                                                                                                                     <  2:222273/275‑3 (MQ=255)
             tatCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                                                                                                                                                                                                     <  2:243334/275‑1 (MQ=255)
                         gATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTCTTAACCTCGCCCAGTAAc                                                                                                                                                                                                                                         >  1:60334/1‑275 (MQ=255)
                                                 tatCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTc                                                                                                                                                                                                                                                              >  2:657160/1‑230 (MQ=255)
                                                     tGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttt                                                                                                                                                                                                                                                          <  4:657158/230‑1 (MQ=255)
                                                             gTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATc                                                                                                                                                                                                                                                                        <  2:291282/208‑1 (MQ=255)
                                                              tCCATTGCGGGGCTGAACCGCAAACCATTCCCCCATGATTAACAAGATACGGTTAGCTGGAAAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGGGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                                                                                                                                                                                     <  4:522236/274‑1 (MQ=255)
                                                              tCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTGTTAACAAGATACGGTTAGCTGGCGAGTTCGAAGCCACTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGAACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGCAAACAGTGGGCTAAAGCGATTTTCAGCCGGGCAAACTCCCTTACACTTGCCGACTCCGCCCATGTCGTTTTAACCTCGCCCAATCACCACCATTAGGCAATCATCAAATCCACCCACAGcttcc                                                                                                                                                                                                    >  4:101760/1‑273 (MQ=255)
                                                              tCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGACAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTTCCGAATCCGGCAATGTCGTTTTTACCTc                                                                                                                                                                                                                                                   >  4:590510/1‑228 (MQ=255)
                                                                  ttGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                                                                                                                                                                                               <  2:590512/228‑1 (MQ=255)
                                                                  ttGCGGCGCTGAACCGCCAACCATGCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                                                                                                                                                                                               <  2:377507/244‑1 (MQ=255)
                                                                     cggcgCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGCGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTACTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGCCGTTTTTACCTAGCCACATAACTACCGTGAGGCCATCCTACAATGCCCCCAACGTTTTCccgcac                                                                                                                                                                                              >  4:128400/1‑272 (MQ=255)
                                                                      ggcgCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAACGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCa                                                                                                                                                                                            <  2:128401/275‑1 (MQ=255)
                                                                       gcgcTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                                                                                                                                                                                           <  2:607667/275‑1 (MQ=255)
                                                                                     aCCATTCTCCCAGTATTAACAAGAGACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAATATCTTCAGTGTTGAGATAATGAGCCAACACTTCCGCCTGTGGCTGTTGGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTa                                                                                                                                                                              <  2:452782/274‑1 (MQ=255)
                                                                                                         aaGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGGTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTGACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAAcc                                                                                                                                                         <  2:11071/275‑1 (MQ=255)
                                                                                                           gATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGcc                                                                                                                                                                                                                >  4:297880/1‑218 (MQ=255)
                                                                                                               cGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCa                                                                                                                                                                                                            <  2:297881/218‑1 (MQ=255)
                                                                                                                ggTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGAGCTTCAGTCTTGAGATAATGGGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc                                                                                                                                                                            <  2:9269/249‑1 (MQ=255)
                                                                                                                              ttCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaa                                                                                                                                                                                                                                                                  >  7:251370/1‑149 (MQ=255)
                                                                                                                              ttCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAACGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaa                                                                                                                                                                                                                                                                  >  7:615157/1‑149 (MQ=255)
                                                                                                                               tCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                                                                                                                                                                                                 <  8:251370/149‑2 (MQ=255)
                                                                                                                               tCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAAt                                                                                                                                                                                                                                                                 <  8:641577/149‑1 (MQ=255)
                                                                                                                                cGAAGCCAGTACACGCTGATGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagt                                                                                                                                                                                                                                                                <  8:615157/149‑3 (MQ=255)
                                                                                                                                    gCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCg                                                                                                                                                                                                                                                             <  5:662922/148‑1 (MQ=255)
                                                                                                                                       aGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtttt                                                                                                                                                                                                                                                         <  7:627844/149‑1 (MQ=255)
                                                                                                                                        gTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCCATGGCGACAACCAAAATg                                                                                                                          >  4:552247/1‑275 (MQ=255)
                                                                                                                                         tACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                                                                                                                                                                                        <  6:445821/148‑1 (MQ=255)
                                                                                                                                             cGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGaa                                                                                                                     <  2:315675/275‑1 (MQ=255)
                                                                                                                                              gCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcg                                                                                                                                                                                                                                                  >  8:660904/1‑149 (MQ=255)
                                                                                                                                                  aaGATCTTCAGTCTTGAGATAATGTGCCAGCATTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAACTCGGCAATGTCGTTTTCACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTTGCGCTTCAACGCAACGGGCGTACCTGGCAATATTACGAAACCCAAGAAAAATTGCGATTGCCAAACCGCGAATGCTTAACCGta                                                                                                                >  4:600050/1‑274 (MQ=255)
                                                                                                                                                   aGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                                                                                                                                                                                                 <  6:76342/129‑2 (MQ=255)
                                                                                                                                                   aGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCa                                                                                                                                                                                                                                                                   >  2:759746/1‑127 (MQ=255)
                                                                                                                                                   aGATCTTCACTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                                                                                                                                                                                                 >  5:76342/1‑128 (MQ=255)
                                                                                                                                                       cTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGt                                                                                                                                                                                                                                                               <  4:759744/127‑1 (MQ=255)
                                                                                                                                                          cAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGa                                                                                                                                                                                                                                       >  5:691585/1‑148 (MQ=255)
                                                                                                                                                           aGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGAc                                                                                                                                                                                                                                      >  5:361471/1‑148 (MQ=255)
                                                                                                                                                              cTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGt                                                                                                                                                                                                                                   >  6:563781/1‑148 (MQ=255)
                                                                                                                                                                  agaTAATGAGCCAGCACTTCCGCCTGTGGCTGTTCGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGCCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                                                                                                                                                                                              >  8:794998/1‑149 (MQ=255)
                                                                                                                                                                   gaTAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCTGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                                                                                                                                                                                               <  8:523941/147‑1 (MQ=255)
                                                                                                                                                                   gaTAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                                                                                                                                                                                               >  7:523941/1‑147 (MQ=255)
                                                                                                                                                                         gAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcat                                                                                                                                                                                                                        >  5:667295/1‑148 (MQ=255)
                                                                                                                                                                          aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcatc                                                                                                                                                                                                                       >  6:413394/1‑148 (MQ=255)
                                                                                                                                                                           gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                                                                                                                                                                                                      <  5:147924/148‑1 (MQ=255)
                                                                                                                                                                           gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCaa                                                                                                                                                                                                                     <  7:7661/149‑1 (MQ=255)
                                                                                                                                                                           gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTaa                                                                                                                                                                                      <  2:221325/180‑1 (MQ=255)
                                                                                                                                                                           gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGAt                                                                                        <  2:699176/274‑1 (MQ=255)
                                                                                                                                                                            caagaaCTTCCGCCTGTGGCTGTTGGCGTAATGTCATGCGTAAAAAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACCCCCTTACACTTGCCGAATCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCg                                                                                     <  2:313789/270‑1 (MQ=255)
                                                                                                                                                                             cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                                                                                                                                                                                                   >  7:477016/1‑149 (MQ=255)
                                                                                                                                                                              aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGt                                                                                                                                                                                                                                                            >  7:376780/1‑107 (MQ=255)
                                                                                                                                                                              aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGt                                                                                                                                                                                                                                                            <  8:376780/107‑1 (MQ=255)
                                                                                                                                                                              aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                                                                                                                                                                                                   >  5:40967/1‑148 (MQ=255)
                                                                                                                                                                                   ttCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccccc                                                                                                                                                                                                             <  8:638817/149‑1 (MQ=255)
                                                                                                                                                                                     ccGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtt                                                                                                                                                                                                                                                           >  8:645718/1‑101 (MQ=255)
                                                                                                                                                                                      cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTCTTACGAAAACCAAGAAAAATGGCACGGGCGGAACCGGAAATCCTGAACGGTCCGAGCGGAAATTTCGTGCCAAGCTGCGTGTACCAg                                                                             >  4:89801/1‑274 (MQ=255)
                                                                                                                                                                                            tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                                                                                                                                                                                     <  6:244137/148‑1 (MQ=255)
                                                                                                                                                                                                   ttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCTTGAGGCAATCATCAACTGCCCCCAAAGTTATCGCGCGTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAAt                                                               <  1:138812/275‑1 (MQ=255)
                                                                                                                                                                                                      cGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                                                                                                                                                                                           <  6:40967/148‑1 (MQ=255)
                                                                                                                                                                                                        tAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACg                                                                                                                                                                                         >  6:749313/1‑148 (MQ=255)
                                                                                                                                                                                                          aTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa                                                                                                                                                                                       >  6:14623/1‑148 (MQ=255)
                                                                                                                                                                                                             tCATGCGTAAACAGTGGTCGAAAGCGATGTGCAGCCGGGCAAACGCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGGTTTCGCGCTTCAACGTAACgg                                                                                                                                                                                   >  8:853115/1‑149 (MQ=255)
                                                                                                                                                                                                             tCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACg                                                                                                                                                                                    <  6:125143/148‑1 (MQ=255)
                                                                                                                                                                                                                 gCGTAAACAGCGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCg                                                                                                                                                                                >  6:287406/1‑148 (MQ=255)
                                                                                                                                                                                                                  cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt                                                                                                                                                                               >  6:539439/1‑148 (MQ=255)
                                                                                                                                                                                                                     aaaCAGCGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc                                                                                                                                                                            <  5:287406/148‑1 (MQ=255)
                                                                                                                                                                                                                       aCAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGt                                                                                                                                                                                                                                   >  2:373612/1‑91 (MQ=255)
                                                                                                                                                                                                                           tGGGCGAAAGCGATTTGCAGCCGGGCAAACTCACTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                                                                                                                                                                                               <  4:373610/91‑1 (MQ=255)
                                                                                                                                                                                                                            gggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGGAACGACGGTGAGGCAATCATCAACTGCCCCCAAAGTTATCGCGCTTCAACGTAACGGGCGTACCCGGCAATATTACGAAAACCAAGAAAAATGGAGAAGGCCACACCGAGAATGCTGAACGATGCGAGGGTGAACTTCATGCCCAGATGAGTTTACCCTGGCAGAATGAAAATAACAACGATAGaaaaacgggcccg                                       >  4:611838/1‑265 (MQ=255)
                                                                                                                                                                                                                                gAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAATGACCGTCAGGCAATCATCAACAGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCCGCATTATTAc                                                                                                                                                                <  8:586491/149‑1 (MQ=255)
                                                                                                                                                                                                                                gAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAg                                                          <  2:611840/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCattatt                                                                                                                                                                  <  6:438592/142‑1 (MQ=255)
                                                                                                                                                                                                                                     cGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCattatt                                                                                                                                                                  >  5:438592/1‑142 (MQ=255)
                                                                                                                                                                                                                                     cGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAg                             <  2:281827/275‑1 (MQ=255)
                                                                                                                                                                                                                                     agATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGaaa                                                                                                                                                            >  6:725268/2‑148 (MQ=255)
                                                                                                                                                                                                                                      gATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTCACCTCGCGCAGTAACGACCGTGAGGCAATCATCACCTGCCCCCAAAGTTTTCTCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAAc                                                                                                                                                          <  8:601698/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATCACGAAAACCaa                                                                                                                                                       <  7:794998/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ccGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGaaaaa                                                                                                                                                 >  8:678812/1‑149 (MQ=255)
                                                                                                                                                                                                                                                cGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAat                                                                                                                                                >  7:68412/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 gggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAat                                                                                                                                                <  6:68956/148‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAat                                                                                                                                                <  7:620029/147‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcg                                                                                                                                            <  6:691585/148‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGTCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGa            <  2:701186/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         tCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatgg                                                                                                                                        >  4:250438/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCattat                                                                                                                                                                   >  5:446131/1‑120 (MQ=255)
                                                                                                                                                                                                                                                          cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatggc                                                                                                                                       >  6:44182/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             ttACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAAc                                                                                                                                                          >  8:451315/1‑126 (MQ=255)
                                                                                                                                                                                                                                                              tACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGacac                                                                                                                                  <  8:211559/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatggc                                                                                                                                       <  8:444794/139‑1 (MQ=255)
                                                                                                                                                                                                                                                               gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatggc                                                                                                                                       >  7:444794/5‑143 (MQ=255)
                                                                                                                                                                                                                                                                 acTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAg  <  2:552249/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ttGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGNCATTATTACGAAAACCAAGAAAANNGGCGATGGNGACAcc                                                                                                                                 >  3:52002/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                   ttGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCgaga                                                                                                                             >  8:613838/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCg                                                                                                                                                                                >  2:687025/1‑96 (MQ=255)
                                                                                                                                                                                                                                                                      ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatg                                                                                                                                         <  5:155607/134‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatg                                                                                                                                         >  6:155607/1‑134 (MQ=255)
                                                                                                                                                                                                                                                                       cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCgaga                                                                                                                             <  1:52003/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGc                                                                                                                         <  7:660904/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaTCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc                                                                                                                                                                            <  4:687023/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGa                                                                                                                      <  5:749313/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGGGAATTTGATGc                                                                                              >  3:108673/1‑172 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
AGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGCTGAACGGTGCGAGGGTGAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAG  >  NC_000913/1608439‑1608969

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