Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 271,822 G→A pseudogene (59/427 nt) insN → IS911A regulator fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913271,8220GA98.4% 178.9 / ‑4.3 63pseudogene (59/427 nt)insNIS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/30);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

GAGAATACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTG  >  NC_000913/271581‑272060
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
atgccgcgaaagctatggatgttggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCtgt                                                                                                                                                                                                                                       <  1:245‑M1/70‑1 (MQ=255)
         aagctatggatgttggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTatcctgtggtcacact                                                                                                                                                                                                                              <  2:192610‑M1/79‑8 (MQ=255)
            ctatggatgttggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTg                                                                                                                                                                                                                           <  1:70426‑M1/82‑1 (MQ=255)
                        gcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCAGTATCCTGTGGTCACACTCTGCCATGTGTTCgg                                                                                                                                                                                                               >  1:314525‑M1/158‑251 (MQ=255)
                           tttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGt                                                                                                                                                                                                            <  2:670982‑M1/97‑1 (MQ=255)
                            ttccacaatgacaagatgcgtcaaacaactgcttgatgatcgtcaggacaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagcttcaacgcattgaaatggataatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGTCTCCCTGTACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGt                                                                                                                                                                                                            >  2:183173‑M1/154‑250 (MQ=255)
                             tccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGtt                                                                                                                                                                                                           >  2:113640‑M1/153‑250 (MQ=255)
                               cacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaattaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCAt                                                                                                                                                                                                        >  1:598867‑M1/151‑251 (MQ=255)
                                acaatgacaagatgggtcaaacaactgcgggatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATc                                                                                                                                                                                                       >  2:157206‑M1/150‑251 (MQ=255)
                                      acaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgcggaaaaagctacaacgcattgaaatggcgaatgaaatattaaaaaggctACCGCGCTCTTGCTGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCgcagc                                                                                                                                                                                                  <  1:85815‑M1/107‑1 (MQ=255)
                                      acaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCgcagc                                                                                                                                                                                                  <  1:493600‑M1/107‑1 (MQ=255)
                                       caagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCgcagc                                                                                                                                                                                                  <  1:382304‑M1/107‑1 (MQ=255)
                                             ggggcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACa                                                                                                                                                                                           >  2:351215‑M1/137‑250 (MQ=255)
                                               gtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCtacagat                                                                                                                                                                                        >  1:35991‑M1/135‑248 (MQ=255)
                                                   aacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTg                                                                                                                                                                                    >  1:587390‑M1/131‑251 (MQ=255)
                                                   aacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaaacgaaatacgtaagcggaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTg                                                                                                                                                                                    >  2:493600‑M1/131‑251 (MQ=255)
                                                    acacctgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGAGGGCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTgg                                                                                                                                                                                   >  2:350598‑M1/130‑251 (MQ=255)
                                                       actgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGaa                                                                                                                                                                                 <  1:520742‑M1/124‑1 (MQ=255)
                                                          gcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGaaaa                                                                                                                                                                               >  1:670982‑M1/124‑249 (MQ=255)
                                                                       cagggcaaaacaccaaaagcctctccgattacacccgaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgcaatggagaatgaaatattaaaaaggctACCGCGCTCTTGCTGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACAATCTGCATTGTGTTCGGGGATCATCGCAGCAGCTACAGTTACTGGAAAAACCGTCAtgaaaacc                                                                                                                                                                <  2:98645‑M1/141‑3 (MQ=255)
                                                                            caaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcagac                                                                                                                                                           <  1:2382‑M1/146‑1 (MQ=255)
                                                                               aacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAACCcagacggc                                                                                                                                                        <  2:160079‑M1/149‑1 (MQ=255)
                                                                                   ccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcagacggcagac                                                                                                                                                    >  2:340009‑M1/99‑251 (MQ=255)
                                                                                    caaaagcctctcctataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcagacggcagacg                                                                                                                                                   >  2:2382‑M1/98‑251 (MQ=255)
                                                                                     aaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcagacggcagacgg                                                                                                                                                  >  1:334266‑M1/97‑251 (MQ=255)
                                                                                        agcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCt                                                                                                                                               >  1:208152‑M1/94‑251 (MQ=255)
                                                                                             ctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTAtt                                                                                                                                          <  1:411672‑M1/163‑1 (MQ=255)
                                                                                               ccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTAtt                                                                                                                                          >  1:260615‑M1/87‑249 (MQ=255)
                                                                                                  ataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattacaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGc                                                                                                                                      <  1:227835‑M1/167‑1 (MQ=255)
                                                                                                  ataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCtgtattccgc                                                                                                                                      <  1:628708‑M1/167‑5 (MQ=255)
                                                                                                         cagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGt                                                                                                                              >  1:272061‑M1/77‑251 (MQ=255)
                                                                                                                    gaaataaataagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGACACACTCTGCCATGTGTTCGGGGTTCATCACCGCAGCTACAGATACTGGAAAACCCATCCTGAAAAACCAGACGGCAGACGGGCTGTattacgcagtcagtgaattgaactc                                                                                                                    <  2:631024‑M1/185‑13 (MQ=255)
                                                                                                                      aatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGgtacttgagctacat                                                                                                                 >  2:382304‑M1/64‑244 (MQ=255)
                                                                                                                             cAAAACACCGAAAGAGATAATGGAAAGGGGTGTTGCAGTGACAGAGTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCGCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTGCGGGGTTCATCGCAGCAGCTACAAACACTGGAAAGACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATc                                                                                                           <  1:292170/250‑1 (MQ=39)
                                                                                                                             aagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATc                                                                                                           <  1:5366‑M1/194‑1 (MQ=255)
                                                                                                                             aagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCa                                                                                                          >  1:160079‑M1/57‑251 (MQ=255)
                                                                                                                                        aaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggttcg                                                                                               <  1:198662‑M1/206‑5 (MQ=255)
                                                                                                                                          agctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGACATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGAAGACGGGCTGTATTACGCAGTCCGGTACTTGAGCTACatggcctcagccccggttcggc                                                                                             <  2:495355‑M1/208‑11 (MQ=255)
                                                                                                                                               caacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAg                                                                                        >  1:224230‑M1/39‑251 (MQ=255)
                                                                                                                                                 acgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGc                                                                                       <  1:286888‑M1/214‑1 (MQ=255)
                                                                                                                                                          aatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACggg                                                                                                                                                 >  1:233165‑M1/28‑183 (MQ=255)
                                                                                                                                                             ggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGcc                                                                          <  1:147007‑M1/227‑1 (MQ=255)
                                                                                                                                                                gaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGAAACTGGAAAAACCGTCCTGAAAAACCCGACGGCAGACGGGCTGAATTACGCAGTAAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATcgccgc                                                                        <  2:573284‑M1/229‑3 (MQ=255)
                                                                                                                                                                   tgaaatattaaaaaggctACAGAGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATg                                                                    <  2:494678‑M1/233‑1 (MQ=255)
                                                                                                                                                                     aaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATCATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTCTTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGc                                                                  <  2:336895‑M1/235‑1 (MQ=255)
                                                                                                                                                                           taaaaatgctTCCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCGCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGTTCACACTCTGCCATTTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTGCGGCCGGAGCAATAAGCATCGCCACAATGGCAAccc                                                             >  2:268440‑M1/11‑250 (MQ=255)
                                                                                                                                                                             aaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCaacccgg                                                           <  1:540371‑M1/242‑3 (MQ=255)
                                                                                                                                                                                aggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCggagag                                                       >  1:336895‑M1/6‑251 (MQ=255)
                                                                                                                                                                                 ggctaCCGCGCTCTTGATGTCCGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCATAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTTGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGAACTTGCGCTACATAGCATCAGACACGGTTCGGCCGGAGCAAGACGCATCGCCACAATGGCCACCCggaga                                                        >  2:643164/5‑249 (MQ=18)
                                                                                                                                                                                 ggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTTTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCggaga                                                        >  1:643164‑M1/5‑249 (MQ=255)
                                                                                                                                                                                     tccGCGCTCTTGATGTCAGTCTCCCTGAACAGTTCACGATAATCGGGAAACGCAGATCGCATTTACCTTTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAggcta                                                   <  2:552735/249‑4 (MQ=17)
                                                                                                                                                                                               tGATGTCAGACTCCCGGAACAGTTCTCGAGAAGCGGGAAACTCAGAGCGCATTATCCTGTGGTCCCACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTTAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCTGGAGAGGCTACCAGATggga                                         <  2:568070/250‑2 (MQ=17)
                                                                                                                                                                                                gctGTCAGACTCCCTGAACAGTGCTCGATAAACGGGAAACTCTGAGCGCCTTCTACTGTGGTCACTCTCTGGCAAGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTTGAAAAACCGTCCTGAAAAACCAGACGGCTGACGGGCTGAATTACGCAGTCAGGTACTTGAGCTACATCGCATCAGCACCGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATgggat                                        <  2:182460/248‑3 (MQ=9)
                                                                                                                                                                                                            acgctctgccgaactaTAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCt                             <  2:255338/234‑1 (MQ=14)
                                                                                                                                                                                                                   aGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTc                    >  1:370654/1‑251 (MQ=14)
                                                                                                                                                                                                                    gTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTc                    >  1:167215/1‑250 (MQ=14)
                                                                                                                                                                                                                         cGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGaa               >  1:442682/1‑250 (MQ=14)
                                                                                                                                                                                                                           aTAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCt                             >  1:255338/1‑234 (MQ=14)
                                                                                                                                                                                                                             aaTCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAagag           <  2:411672/250‑1 (MQ=14)
                                                                                                                                                                                                                                cGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTgg       >  1:649254/1‑251 (MQ=14)
                                                                                                                                                                                                                                  ggAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTg                                    <  1:271801/220‑1 (MQ=14)
                                                                                                                                                                                                                                  ggAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTgggg     >  1:389421/1‑251 (MQ=11)
                                                                                                                                                                                                                                     aaCTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGGTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTGCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACAAGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTg  <  1:4290/251‑1 (MQ=11)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
GAGAATACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTG  >  NC_000913/271581‑272060

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: