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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913271,8170GA95.0% 164.6 / ‑1.3 60pseudogene (54/427 nt)insNIS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (34/23);  total (36/24)
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.39e-01

GGGAACAAATGAGAATACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGG  >  NC_000913/271571‑272061
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
acggtggcagatgccgcgaaagctatggatgttggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAgcgcgc                                                                                                                                                                                                                                                  <  2:53756‑M1/60‑1 (MQ=255)
                cgaaagctatggatgttggcctttccacaatgacaagatgggtcaaacacctgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtacgctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGACCAGTTCTCGATAATCAGGAAAATCAAAGCGCATtatcctgtgctcac                                                                                                                                                                                                                                  <  2:275076‑M1/76‑6 (MQ=255)
                        atggatgttggactgtccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaacgcctctccgctaacaccagaacaaatcgaaatacgtaagctgaggaaaaagatacaacgcattgaaatggagaatgaaaaattaaaaaggcaAACGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTAGATATTCGGGAAACTCAGAGCgcattatcatgtggtcaacctctgcc                                                                                                                                                                                                                          <  2:403175‑M1/84‑13 (MQ=255)
                            atgttggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGTACACACTCTGCCAtgt                                                                                                                                                                                                                      <  1:279449‑M1/88‑1 (MQ=255)
                                tggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGtt                                                                                                                                                                                                                   >  1:317729‑M1/160‑250 (MQ=255)
                                 ggcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggcaACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCCGGATAATCGGGAAACTCAAAGCGCATTATCCTGTGGTCACACTCTGCCATGtgttgg                                                                                                                                                                                                                 <  1:270384‑M1/93‑3 (MQ=255)
                                  gcctttccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCgg                                                                                                                                                                                                                >  1:245938‑M1/158‑251 (MQ=255)
                                   cctttccacaatgacaagatgggtcacacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCggg                                                                                                                                                                                                               <  2:417687‑M1/95‑1 (MQ=255)
                                       tccacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGtt                                                                                                                                                                                                            >  1:293144‑M1/153‑250 (MQ=255)
                                         cacaatgacaagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCAt                                                                                                                                                                                                         >  2:51665‑M1/151‑251 (MQ=255)
                                             atgacaagattggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCgca                                                                                                                                                                                                     >  1:101230‑M1/147‑251 (MQ=255)
                                               gacaagatgggtcaaacaactgcgtgacgagcgtcagggcaaaacaccaaaagtctctccgataacaccagaacaaagcgaaatacggaagctgaggaaaaagctacaactcatttaaatggagaatgaaatattaaaaaggctACCGCGTTCTTGATGTCAGACTCCCTGAACAGTCCTCGATAAGCGCGAAACTCTGAGCGCATTATCCTGTGGTCACCCTCTGCCTTGTGTTCGGGGTtcatctcag                                                                                                                                                                                                    >  2:174889‑M1/145‑246 (MQ=255)
                                                           caaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATAc                                                                                                                                                                                       >  1:449792‑M1/133‑251 (MQ=255)
                                                              acaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagcgacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTgg                                                                                                                                                                                    >  2:314314‑M1/130‑251 (MQ=255)
                                                                aactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCagcagctacagataatggc                                                                                                                                                                                   <  2:280190‑M1/123‑10 (MQ=255)
                                                                aactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGaa                                                                                                                                                                                  >  1:280190‑M1/128‑251 (MQ=255)
                                                                       tgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgacatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCCCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCGGCTACAGATACTGGAAAAACcgg                                                                                                                                                                           <  2:427201‑M1/131‑2 (MQ=255)
                                                                        tCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCCCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTc                                                                                                                                                                          >  1:251442/1‑251 (MQ=255)
                                                                           gagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCACGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCAGGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTg                                                                                                                                                                       >  1:417687‑M1/117‑251 (MQ=255)
                                                                             gcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGaa                                                                                                                                                                     >  1:437738‑M1/115‑251 (MQ=255)
                                                                                  agggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCctgaaaaacc                                                                                                                                                                >  1:345417‑M1/110‑248 (MQ=255)
                                                                                         aactccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcatggaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcagacggc                                                                                                                                                         >  2:229638‑M1/103‑251 (MQ=255)
                                                                                               aaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctaaaaagcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcaaacggcagacgg                                                                                                                                                   <  1:342424‑M1/155‑1 (MQ=255)
                                                                                                         ccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTa                                                                                                                                          >  1:53756‑M1/87‑250 (MQ=255)
                                                                                                                 accagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTc                                                                                                                                   >  2:60738‑M1/79‑249 (MQ=255)
                                                                                                                 accagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTc                                                                                                                                   <  1:60738‑M1/171‑1 (MQ=255)
                                                                                                                    ggaacaaatcgaaatacgtacgctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCGGACTCCCTGAACAGTTCTCGATATTCGGGAAACTCAGAGAGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGCAAAACCGTCCTGAAGAACCAGCCGGCAGACGGGctgaattacgagatcaggt                                                                                                                               <  2:95146‑M1/175‑10 (MQ=255)
                                                                                                                           atcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctAACGCGCTCTTGATGTCAGACTCCTTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGTATTACCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGTAGTCAGGTACTtgagc                                                                                                                       >  1:95146‑M1/69‑249 (MQ=255)
                                                                                                                             cgagatacgtaagatgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACACCCAGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCATTATCCTGTGGTAACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGCCGTGCTGCATTACGCAGTCAGGAACttgagct                                                                                                                      <  2:75898‑M1/184‑4 (MQ=255)
                                                                                                                                     gtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCAt                                                                                                             <  1:376622‑M1/193‑1 (MQ=255)
                                                                                                                                     gtaagctgaggaaaaagctacaacgcattgaaagggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCTATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTGCGGGGTTCATCGCAGCAGCTACAGATTCTGGAAAAACCGTCCTTAAAAACCAGACGGCAGACGTGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCAt                                                                                                             >  1:463123‑M1/59‑251 (MQ=255)
                                                                                                                                      taagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATc                                                                                                            <  1:366332‑M1/194‑1 (MQ=255)
                                                                                                                                               gaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCATGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGCGCTACATGGCATCAGCCAcggt                                                                                                   <  1:304605‑M1/203‑2 (MQ=255)
                                                                                                                                                  aaagctacaacgcattgacatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAACAACCAGACGGCAGCCGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggatcg                                                                                                <  2:189247‑M1/206‑5 (MQ=255)
                                                                                                                                                           acgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTAGCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCa                                                                                       >  2:386843‑M1/37‑251 (MQ=255)
                                                                                                                                                            cgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAACCCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCa                                                                                       <  2:248114‑M1/215‑1 (MQ=255)
                                                                                                                                                                        gagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGcc                                                                           <  1:149895‑M1/227‑1 (MQ=255)
                                                                                                                                                                        gagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCca                                                                          >  1:115820‑M1/24‑251 (MQ=255)
                                                                                                                                                                          cATTGACAGATTGACTCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAACGCCCGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCAACGcccc                                                                         >  2:48899/1‑248 (MQ=255)
                                                                                                                                                                                aatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcaga                                                                                                                                                             <  1:358349‑M1/145‑1 (MQ=255)
                                                                                                                                                                                aatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACcaga                                                                                                                                                             >  2:358349‑M1/16‑160 (MQ=255)
                                                                                                                                                                                   attaaaaaggctACCGCGCTCTTGATGTCAGTCTCCCTGAACAGTTCGCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAAcc                                                               >  2:366332‑M1/13‑251 (MQ=255)
                                                                                                                                                                                        aaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAACAACCGTCCTGAAACACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGAACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCgccccaatggaaacacgg                                                            <  2:216701‑M1/242‑9 (MQ=255)
                                                                                                                                                                                         aaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCggaga                                                         >  1:291973‑M1/7‑251 (MQ=255)
                                                                                                                                                                                          aggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCAtcagccacggttcggccg                                                                                            >  2:157561‑M1/6‑207 (MQ=255)
                                                                                                                                                                                          aggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCAtcagccacggttcggccg                                                                                            <  1:157561‑M1/210‑9 (MQ=255)
                                                                                                                                                                                              tACCGCGCTCTTGATGTCAGACTCCATGAACAGTTCTCGATAATCGGGAAAATCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGGTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAAACAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTAAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGGAACCCGGAgaggct                                                     <  1:351208‑M1/249‑3 (MQ=255)
                                                                                                                                                                                                ccGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGACACTCAGAGCGAATTATCCTGTAGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCTTGAAAAAACAGACGGCAGACGGGCTGTATTCCGCAGTAAGGTACTTGCACTACATGGCATAAGCCACGGTTTGGCCGGAGCAAGAAGCATCGCCACACGGGCAACCCagagcggcttcc                                                  >  2:60439/1‑243 (MQ=255)
                                                                                                                                                                                                ccGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCa                                                                                                           <  2:228563/194‑1 (MQ=25)
                                                                                                                                                                                                ccGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCa                                                                                                           >  1:228563/1‑194 (MQ=25)
                                                                                                                                                                                                 cgcgCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTAcc                                                  >  1:156408/1‑250 (MQ=18)
                                                                                                                                                                                                     ctctTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATg                                             >  2:79961/1‑251 (MQ=18)
                                                                                                                                                                                                          gATGTCAGACTCCCTGAACAGTGCCCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATgggacg                                        <  1:427201/251‑4 (MQ=17)
                                                                                                                                                                                                                aGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCt                                  >  1:66770/1‑251 (MQ=17)
                                                                                                                                                                                                                           aCAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAggc                       >  1:495618/1‑249 (MQ=14)
                                                                                                                                                                                                                               ttCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCAACAATGGCACACCGGAGAGGCTACCAGATGAGACGCTTGCTTACTGGCAggcaca                    >  2:131604/1‑245 (MQ=14)
                                                                                                                                                                                                                                 ctcGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATg                  >  1:425928/1‑250 (MQ=14)
                                                                                                                                                                                                                                            ggAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCGGGCTCATGAAAGAGCTgggg      >  2:280533/1‑251 (MQ=11)
                                                                                                                                                                                                                                              aaaCTCAGAGCTCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCt    <  2:156408/251‑1 (MQ=11)
                                                                                                                                                                                                                                               aaCTCAGAGCGAATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAACACCGTACTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGAATCAGCCACGGTTCGGCCGGAGAAAGAAGCATCGCCACAATGGCAACCCGGAGCGGCTACCAGATGGGACGCTGGCTTGCTGGACGGCCCATGAAAGAGCTgggggtg   >  2:459821/1‑248 (MQ=11)
                                                                                                                                                                                                                                                aCTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCGCATGAAAGAGCTGGGGCTgg  >  1:143196/1‑251 (MQ=11)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GGGAACAAATGAGAATACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGG  >  NC_000913/271571‑272061

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