Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,459 A→G 8.3% S103G (AGC→GGC)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,4590AG 8.3% 793.6 / 41.0 255S103G (AGC→GGC) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (126/106);  new base G (9/13);  total (135/120)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 260 total aligned reads displayed.

ATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGG  >  NC_000913/3991203‑3991736
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                     
            tCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCct                                                                                                                                                                                                                                                         <  2:124042/275‑1 (MQ=255)
                  gcTTGCGTCTATGGGGCCTGCATTCAAACAGGTCTAAAGGTGACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGTTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGAc                                                                                                                                                                                                                                                    <  4:562104/274‑1 (MQ=255)
                            aTGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTg                                                                                                                                                                                                                                         <  2:230996/275‑1 (MQ=255)
                                                ggTTTACGGTCTACTGCCGACATTGTTGCCCTTACACCATCCGCTACTGCGGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTgg                                                                                                                                                                                                                      <  4:589936/274‑1 (MQ=255)
                                                           tACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGaa                                                                                                                                                                                                          >  2:410634/1‑275 (MQ=255)
                                                           tACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTCACACGCCTGATGAAACTCAACGCTACTACCTGAACGAGCTTGAACTGTATCGTGGAATCTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCTTCTGCCCTGGATATCTGGGTCTGTCATCAACCCTGGCCCGATTGCCaa                                                                                                                                                                                                          >  4:185712/1‑275 (MQ=255)
                                                                     aTTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAAt                                                                                                                                                                                                >  2:571293/1‑275 (MQ=255)
                                                                           gCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTAc                                                                                                                                                                                          <  2:185712/275‑1 (MQ=255)
                                                                               tATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGCAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACaa                                                                                                                                                                                        >  2:653018/1‑273 (MQ=255)
                                                                                     cATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGTaa                                                                                                                                                                                                                                                                   <  4:46272/192‑1 (MQ=255)
                                                                                       tCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCctc                                                                                                                                                                                                                                                        >  2:339858/1‑201 (MQ=255)
                                                                                           cTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGAc                                                                                                                                                                                                                                                    <  4:339858/201‑1 (MQ=255)
                                                                                                   ggattACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTCCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTg                                                                                                                                                                                                                                                      >  4:768488/4‑191 (MQ=255)
                                                                                                              gtGGAAACTTTTCAAAAGGCATTGTGCTGTACNCGCCGGGCGAAACTAAACGCCACTCNTGTAGCGAGCTTGAAATGATTCATGGAATGTCAGGAAAGGATCCGCGGACAGGTGAGCTTCAAATTACTGGTGTAGACACNATGGGCAGCACCTCGGCCGTAGGGGAAAGTTGGTCCTGTGACCTGGATATCTGGGGCTGTCAGAAATCCTTGCTCGATAGCGAAGGGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCcgcg                                                                                                                                                    <  3:126647/277‑1 (MQ=255)
                                                                                                                ggTAACGTTCCCAAAGGCATTTGCCTTTACCCGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGAGCCGGCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAc                                                                                                                                                                                                                                                                                   <  6:555833/149‑1 (MQ=255)
                                                                                                                   aaCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTc                                                                                                                                                                                                                                                                                <  6:66843/149‑1 (MQ=255)
                                                                                                                    aCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTCAATGTAGCCTGCTGGAAAACTGGGCCGCCTCCCt                                                                                                                                                  >  2:157210/1‑274 (MQ=255)
                                                                                                                       ttCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc                                                                                                                                                                                                                                                                             >  7:953811/1‑148 (MQ=255)
                                                                                                                       ttCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc                                                                                                                                                                                                                                                                             >  7:325711/1‑148 (MQ=255)
                                                                                                                          ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTTAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGt                                                                                                                                                                                                                                                                          <  8:344697/148‑1 (MQ=255)
                                                                                                                          ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCCCCTCGTCCGTa                                                                                                                                                                                                                                                                         <  6:534255/149‑1 (MQ=255)
                                                                                                                          ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGt                                                                                                                                                                                                                                                                          <  8:950546/148‑1 (MQ=255)
                                                                                                                          ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGt                                                                                                                                                                                                                                                                          <  7:282772/148‑1 (MQ=255)
                                                                                                                          ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTa                                                                                                                                                                                                                                                                         >  6:502497/1‑149 (MQ=255)
                                                                                                                           taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAg                                                                                                                                                                                                                                                                        >  6:233359/2‑149 (MQ=255)
                                                                                                                           cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTa                                                                                                                                                                                                                                                                         <  8:242018/148‑1 (MQ=255)
                                                                                                                            aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTgg                                                                                                                                                                     >  4:207891/1‑247 (MQ=255)
                                                                                                                            aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAgg                                                                                                                                                                                                                                                                       >  5:228406/1‑149 (MQ=255)
                                                                                                                            aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCAGGGGCGGCACCTCGTCCGTAg                                                                                                                                                                                                                                                                        <  7:622407/148‑1 (MQ=255)
                                                                                                                             aaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAgg                                                                                                                                                                                                                                                                       >  8:87264/1‑148 (MQ=255)
                                                                                                                               ggCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCa                                                                                                                                                                                                                                                                    >  6:1001434/1‑149 (MQ=255)
                                                                                                                                gCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCa                                                                                                                                                                                                                                                                    >  7:1046677/1‑148 (MQ=255)
                                                                                                                                 cATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCaaa                                                                                                                                                                                                                                                                  <  5:644419/149‑1 (MQ=255)
                                                                                                                                 cATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCaa                                                                                                                                                                                                                                                                   >  8:765471/1‑148 (MQ=255)
                                                                                                                                   tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTGTCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCaa                                                                                                                                                                                                                                                                   <  7:401856/146‑1 (MQ=255)
                                                                                                                                   tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCaa                                                                                                                                                                                                                                                                   >  8:401856/1‑146 (MQ=255)
                                                                                                                                        cTTTATACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAATGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTc                                                                                                                                                                                                                                                           <  5:303958/149‑1 (MQ=255)
                                                                                                                                        cTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGCGCTTGAACTGTATCGTGCAATGTCAGTACAGGATCCGCCGAAAGGTGAGGTTCCAAATACTGGTGTATAGACCATGGGCAGCACCTCGACCGTAGGGCAACGTTGTTc                                                                                                                                                                                                                                                           <  6:714689/149‑1 (MQ=255)
                                                                                                                                        cTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTc                                                                                                                                                                                                                                                           >  5:878702/1‑149 (MQ=255)
                                                                                                                                        cTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTc                                                                                                                                                                                                                                                           <  5:1031921/149‑1 (MQ=255)
                                                                                                                                         tttACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTc                                                                                                                                                                                                                                                           >  7:712645/1‑148 (MQ=255)
                                                                                                                                                 ccTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCt                                                                                                                                                                                                                                                  <  6:317394/149‑1 (MQ=255)
                                                                                                                                                     atgaAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGAc                                                                                                                                                                                                                                                    >  5:496565/1‑143 (MQ=255)
                                                                                                                                                       gaAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGatat                                                                                                                                                                                                                                            >  5:774451/1‑149 (MQ=255)
                                                                                                                                                        aaaCTCAACGCCACTACCTGAACGAGCTTGAACTGAATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGatat                                                                                                                                                                                                                                            >  8:1001278/1‑148 (MQ=255)
                                                                                                                                                           cTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTg                                                                                                                                                                                                                                         <  8:1046677/148‑1 (MQ=255)
                                                                                                                                                               aCGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTc                                                                                                                                                                                                                                     <  7:383265/148‑1 (MQ=255)
                                                                                                                                                               aCGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCt                                                                                                                                                                                                                                    <  5:590135/149‑1 (MQ=255)
                                                                                                                                                                  ccACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTc                                                                                                   <  2:276074/275‑1 (MQ=255)
                                                                                                                                                                   cACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAgttgtt                                                                                                                                                                                                                                                            <  7:915837/121‑1 (MQ=255)
                                                                                                                                                                   cACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAgttgtt                                                                                                                                                                                                                                                            >  8:915837/1‑121 (MQ=255)
                                                                                                                                                                   cACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGtca                                                                                                                                                                                                                                >  6:429265/1‑149 (MQ=255)
                                                                                                                                                                     cTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCACGATAGCGAAGAGCGCCAATTGCTa                                                                                                                                                                                           <  1:68336/184‑1 (MQ=255)
                                                                                                                                                                      tACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAAGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc                                                                                                                                                                                                                                                                             >  6:131545/1‑101 (MQ=255)
                                                                                                                                                                         cTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATc                                                                                                                                                                                                                          >  5:623271/1‑149 (MQ=255)
                                                                                                                                                                             aaCGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGCAAGGTTAGCTTCCAATTACTGGTGTATACATCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTATGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGTAGAGCGCCAATTGCTACAACGTAAATGTAACCTGCTGCAACACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCTGGATTGATGAAACCCCCTTCCGTCATAATGAAAg                                                                                         >  4:293247/1‑275 (MQ=255)
                                                                                                                                                                             aaCGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCt                                                                                                                                                                                                                        <  7:166043/148‑1 (MQ=255)
                                                                                                                                                                             aaCGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTg                                                                                                                                                                                                                       >  5:425071/1‑149 (MQ=255)
                                                                                                                                                                                  gCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATc                                                                                                                                                                                                                          >  6:604156/1‑141 (MQ=255)
                                                                                                                                                                                   cTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGctt                                                                                                                                >  2:727141/1‑230 (MQ=255)
                                                                                                                                                                                   cTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAg                                                                                   >  4:527385/1‑275 (MQ=255)
                                                                                                                                                                                    ttGAAGTGGATCGTTGAATGGAAGTGAAGGATCCGCCGAAAGGTGAGCTTCCAATTACGGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGAAAAGTTGTTCCTCTGAACTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTgg                                                                              <  1:127738/279‑1 (MQ=255)
                                                                                                                                                                                     tGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATa                                                                                                                                                                                                               <  5:1050830/149‑1 (MQ=255)
                                                                                                                                                                                          tGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGaa                                                                                                                                                                                                          >  6:493240/1‑149 (MQ=255)
                                                                                                                                                                                            tATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCct                                                                                                                                                                                                                                                         >  1:33608/1‑100 (MQ=255)
                                                                                                                                                                                             aTCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAagag                                                                                                                                                                                                       <  6:900953/149‑1 (MQ=255)
                                                                                                                                                                                             aTCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAaga                                                                                                                                                                                                        >  8:437806/1‑148 (MQ=255)
                                                                                                                                                                                                 tGGAATGTCAGTACAGGATCCGCCGAAAGGTTAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGATAGCGGCAGCCAGGGGGGCGaa                                                                      >  2:611161/1‑274 (MQ=255)
                                                                                                                                                                                                      tGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCCCCTCCCTGGGTGTGGAATTAAGCTTATTCCTGATGGCTTACAACCGCCTCGGAAAAAATNAAAAANNCAGCGGCGGGGGGGaaaggtggtggc                                                            >  3:17541/1‑269 (MQ=255)
                                                                                                                                                                                                       gTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGc                                                                                                                                                                                             >  5:363579/1‑149 (MQ=255)
                                                                                                                                                                                                        tCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGaa                                                                                                                                                                                                          >  5:907199/1‑135 (MQ=255)
                                                                                                                                                                                                            tACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACa                                                                                                                                                                                         >  7:427766/1‑148 (MQ=255)
                                                                                                                                                                                                              ctggtacccccgAAAGGTGTGTGTCCAATTCCTTGTGGATACACCATGGGGAGCACCGCGTCCGTAGGGCAAAGTTGTTCCTCTGAGCTGGATATCTGGGTCTGTCATCAATCCTTGCTCGATAGCGAAGAGCGCCAATTGCTACAACGGAAATGTAGCCTGCTGGAAAACTGGGCCACCTCGCTGGGTGTGTAAGTCAGCTTCTTCCGGATTGCTGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTTGGGGGCGAAGATTGTGGCTCCAc                                                        <  4:157210/266‑1 (MQ=255)
                                                                                                                                                                                                               aGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACg                                                                                                                                                                                      <  8:888101/148‑1 (MQ=255)
                                                                                                                                                                                                                  gtCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTaaa                                                                                                                                                                                  >  6:913075/2‑149 (MQ=255)
                                                                                                                                                                                                                  aTCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTaa                                                                                                                                                                                   <  7:724514/148‑1 (MQ=255)
                                                                                                                                                                                                                      gccgAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTCTGGCTCCACCCAGCata                                                >  2:469548/1‑275 (MQ=255)
                                                                                                                                                                                                                               tGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTgg                                                                                                                                                                     >  5:164787/1‑149 (MQ=255)
                                                                                                                                                                                                                                gAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTcc                                                                                                                                                                                                                                                          <  6:212725/63‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCTTCCAATTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGaa                                                                                                                                                                   >  6:357382/1‑149 (MQ=255)
                                                                                                                                                                                                                                     tCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAAt                                                                                                                                                                                                                           <  7:939804/89‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTTACTGGTGTATACACCATGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTgg                                                                                                                                                             >  7:554496/1‑148 (MQ=255)
                                                                                                                                                                                                                                        aaTTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTgg                                                                                                                                                             <  8:67636/148‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGgc                                                                                                                                                           <  6:970506/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGgc                                                                                                                                                           <  5:954934/149‑1 (MQ=255)
                                                                                                                                                                                                                                          ttACTAGGGTATACACCATGGGCAGACTCTGGTCCGTAGGGAAAAGTTGTTCCTCTGACCTGGAAACCTGGGTCTGTCATCAGTCCTGGCTGGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAAttt                            <  4:25002/275‑1 (MQ=255)
                                                                                                                                                                                                                                             atGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAg                                                                                                                                                                             >  7:237386/2‑127 (MQ=255)
                                                                                                                                                                                                                                               gggttttatacCCTGGGCGTCACCCGGTCCGTAGGGCAAAGTGGTTCCCCTGACCTGGGTATATGGGTCTTGCATCCATCCTGGCTCGATAGCGAAAAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTCATGAAAACCGCTTCCGGCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATc                        <  4:611161/265‑1 (MQ=255)
                                                                                                                                                                                                                                                        cccATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTGCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGgtg                                                                                                                                             <  8:230746/147‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccATGGGCGGCACCTCGTCCGAAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTATGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTCCAACGTAGATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTgtgtgt                                                                                                                                            >  8:660090/1‑148 (MQ=255)
                                                                                                                                                                                                                                                         ccATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTg                                                                                                                                                 <  8:562577/143‑1 (MQ=255)
                                                                                                                                                                                                                                                          cATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATtgt                                                                >  2:689461/1‑223 (MQ=255)
                                                                                                                                                                                                                                                          cATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACGTCTGCTTGCCGAATTTTATCGTACCGCCGttc            >  4:253232/1‑273 (MQ=255)
                                                                                                                                                                                                                                                           aTGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGAGTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCTTATACTTCTGCTTGACGAATTTTATCGTACCGCCGTACg           >  2:366056/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            tGGGCGGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAg                                                                                         >  4:85602/1‑196 (MQ=255)
                                                                                                                                                                                                                                                            tGGGCATCACGTCTTACGTAGGGCAAAGTTGTTCCGCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGtctc       <  1:14121/278‑1 (MQ=255)
                                                                                                                                                                                                                                                            tGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGaa                                                                                                                                        >  5:525211/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             gggCAGCACCTCGTCCGTAGGGCAAAGTTGTTNNNCGTACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATGGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCg      <  3:217433/278‑1 (MQ=255)
                                                                                                                                                                                                                                                             gggCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTNGGTGTGGAAGTCAGCTTCTTCCGGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGGNGATGGTGGCGCCACCCAGCAGATACTGCTGNNTGACGAAATTTACCGTAACGCCGTGCGtaccng     >  3:96920/1‑274 (MQ=255)
                                                                                                                                                                                                                                                             gggCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATAGTCCCGCCGTGCGTCTCg      >  1:217433/1‑278 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                     
ATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGG  >  NC_000913/3991203‑3991736

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