Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913271,8170GA48.1% 13.2 / 61.9 55pseudogene (54/427 nt)insNIS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base A (14/12);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGG  >  NC_000913/271588‑272055
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
gaaagctatggatgttggcctttacacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatgtagaatgaaatattaaaaaggctACCGTGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATcctgtggtcaca                                                                                                                                                                                                                           >  1:186047‑M1/175‑246 (MQ=255)
 aaagctatggatgttggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagccgctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTAtcctgtggtcacac                                                                                                                                                                                                                          >  1:164738‑M1/174‑245 (MQ=255)
     gggCTAATTCGGCAGTACTTTCCTAAAAAGACATGGCTTACCCAATATACTCAACATGAACGAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTg                                                                                                                                                                                                                      <  2:120859/251‑1 (MQ=255)
        gggatgttggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaatgctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACtctgaca                                                                                                                                                                                                                   >  2:236411‑M1/167‑248 (MQ=255)
         ggagttttgcctttccacaataacaagatgtgtcaaacaactgcgtgatgagcgtcaggtcaaaacaccaaaatcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCTATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTATGCCAt                                                                                                                                                                                                                  >  2:314148‑M1/166‑251 (MQ=255)
                  aGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCggg                                                                                                                                                                                                         <  1:231616/251‑1 (MQ=255)
                     aCTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGt                                                                                                                                                                                                       <  1:52831/250‑1 (MQ=255)
                       ccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCa                                                                                                                                                                                                    >  1:28445‑M1/152‑251 (MQ=255)
                                     aTGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTAc                                                                                                                                                                                       >  1:98392/1‑250 (MQ=255)
                                       gTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACaa                                                                                                                                                                                     <  1:153674/250‑1 (MQ=255)
                                                    taCTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAAcc                                                                                                                                                                       >  1:349356/1‑251 (MQ=255)
                                                            aTGAACTAGCTCTGGTTGCTGCTCAGCTAAAAAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACCAGCGCTCTTGATGTCAGACTCCCTGAACAGATCTCGAGACTCGGGAAACTCAGAGCGCGTTATCATTTGGCCACTCTCTGCCATGTGTTAGTGGTTCCTCGCACCGACTACAAATCCTGGACAAACCTTCCtta                                                                                                                                                                >  2:314477/1‑248 (MQ=18)
                                                                         ggTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCcaga                                                                                                                                                       <  1:287003/246‑1 (MQ=255)
                                                                         ggTTGCTACTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCcagacggca                                                                                                                                                  <  1:312704/251‑1 (MQ=255)
                                                                               tgctCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTa                                                                                                                                                                                                                                        <  1:54073/159‑1 (MQ=255)
                                                                               tgctCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACAGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTa                                                                                                                                                                                                                                        >  2:54073/1‑159 (MQ=255)
                                                                               tgctCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGAAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTa                                                                                                                                                                                                                                        >  2:35019/1‑158 (MQ=255)
                                                                               tgctCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGAAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTa                                                                                                                                                                                                                                        <  1:35019/158‑1 (MQ=255)
                                                                                     tctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattggactggagactgaaatattaaaaaggctACCGCGCTCTTGATGTCAGCCTCCCTGAACAGTTCTCGATAATCGTGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGCTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTAt                                                                                                                                      <  1:179167‑M1/162‑1 (MQ=255)
                                                                                     gCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTAt                                                                                                                                      >  1:267216/1‑251 (MQ=255)
                                                                                      cTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTAtt                                                                                                                                     >  1:269551/1‑251 (MQ=255)
                                                                                       tAAACACCAGACCGAGAAATACACTTAATTTCAAAACACTGAAAGAAATAAGTGAAAGGGTGGTTGCATTGACAGATTGAAGCTACAACCGTGCTCTTGATGTCAGATTCCCTTAACAGTTCTCAATAAACGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTTAAAAGCCAGACGGCAGACGGGCTGTAtt                                                                                                                                     <  2:300476/250‑1 (MQ=255)
                                                                                          acaacaGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGc                                                                                                                                 >  1:86772/1‑251 (MQ=255)
                                                                                             aacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAg                                                                                                                               >  1:88109‑M1/82‑250 (MQ=255)
                                                                                                 ccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgccttgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAgg                                                                                                                          <  1:313635‑M1/174‑1 (MQ=255)
                                                                                                       caaatcgaaatacgtaagatgaggaaaaagctacaacgcattgaaatggcgaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGGGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTACTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTg                                                                                                                    <  1:154738‑M1/180‑1 (MQ=255)
                                                                                                        aaGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGa                                                                                                                   <  1:282080/251‑1 (MQ=255)
                                                                                                          atcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTtgagc                                                                                                                 >  1:316444‑M1/69‑249 (MQ=255)
                                                                                                          atcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTtgagc                                                                                                                 <  1:237092‑M1/183‑3 (MQ=255)
                                                                                                                 tacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATgg                                                                                                          <  1:80496‑M1/190‑1 (MQ=255)
                                                                                                                     taagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATc                                                                                                      <  1:216110‑M1/194‑1 (MQ=255)
                                                                                                                      cAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATc                                                                                                      >  1:190713/1‑250 (MQ=255)
                                                                                                                                aaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCcacggttc                                                                                           >  1:294588‑M1/47‑248 (MQ=255)
                                                                                                                                    agaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGCCTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACTGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGcc                                                                                       <  1:36712/251‑1 (MQ=255)
                                                                                                                                                   gggtTGCATTGACAGATTGAATCTACAACCGCGTTCTTGATGTCAGGCTCCCTGAACAGTTCTCGATAATCGGGAATCTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCTGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATc                                                                        <  1:366867/249‑1 (MQ=255)
                                                                                                                                                      ttGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGcc                                                                     >  1:180769/1‑251 (MQ=255)
                                                                                                                                                         gaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCcaca                                                                  >  1:59982‑M1/22‑251 (MQ=255)
                                                                                                                                                          aatgaaatattaaaaaggctaccacgcTATTGATGTCAGACTCCCTGAACAGTTCTCGATAAACGGGAAACTCAGCGCGAATATTACTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCGGCTACAGATCCTGGAAAAACCGTCCTGAAACACCAGACGGCAGACGGGATGAATTGCGAAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGTCAGGAGCAGAAAGAATCGCcaca                                                                  >  2:275308/25‑250 (MQ=12)
                                                                                                                                                            tgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCAACGGTTCGGCCGGAGAAAGAAGCATCGCCACAATg                                                               <  1:350045‑M1/233‑1 (MQ=255)
                                                                                                                                                               aatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCa                                                            >  1:25614‑M1/16‑251 (MQ=255)
                                                                                                                                                                  aTTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCGCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAAcc                                                         <  1:252634/251‑1 (MQ=255)
                                                                                                                                                                   ttGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAAccc                                                        >  1:75415/1‑251 (MQ=255)
                                                                                                                                                                       aaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCaacccgg                                                      <  1:79929‑M1/242‑3 (MQ=255)
                                                                                                                                                                             tACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAgaggct                                               <  1:184237‑M1/249‑3 (MQ=255)
                                                                                                                                                                              aCCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTAc                                             >  1:381293/1‑251 (MQ=18)
                                                                                                                                                                               ccGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTAGATAATCGGGAAACTCAGAGCGCATTATCCGGTGGTCACACTCTGACATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTAcc                                            <  2:79929/251‑1 (MQ=18)
                                                                                                                                                                                  cgcTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTAcc                                            >  1:284313/1‑248 (MQ=18)
                                                                                                                                                                                     tctTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATgg                                      >  1:371688/1‑251 (MQ=255)
                                                                                                                                                                                         gATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATgggacg                                  <  1:377654/251‑4 (MQ=17)
                                                                                                                                                                                           ggtcatgctcgttGAACAGTTCTCGAAAAACGGGAAATTCATAGCGCATTTTCCTGCGGTCACACTCGGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATAATGGAAAAACCTTCCTGAAAAACCAGGCGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGTAAGAAGCAGCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCt                                <  2:179167/239‑1 (MQ=2)
                                                                                                                                                                                                        gAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGCTGGGACGCTGGCTTGCTGGc                     >  1:217417/1‑249 (MQ=17)
                                                                                                                                                                                                             gTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCGGATGGGGCGCTGGCTTGCTGGCAGACTc               >  2:163310/1‑250 (MQ=255)
                                                                                                                                                                                                                  cGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGaa          >  1:166910/1‑250 (MQ=255)
                                                                                                                                                                                                                         cGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTgg  <  1:298163/251‑1 (MQ=255)
                                                                                                                                                                                                                          gtgAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCaga                                                     <  1:352679/197‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
CAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGG  >  NC_000913/271588‑272055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: