Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 364,170 2 bp→AA coding (2135‑2136/3075 nt) lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913364,1700GA95.0% 117.4 / ‑5.5 40A712V (GCC→GTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/29);  minor base C (2/0);  total (11/29)
Fisher's exact test for biased strand distribution p-value = 7.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01
*NC_000913364,1710GA100.0% 122.1 / NA 40A712V (GCC→GTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/29);  total (11/29)

ATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGG  >  NC_000913/363957‑364408
                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                              
agcgcgCAGCGGCGGCAGCAGTTGTGTGTTATCGCCAATACACATGTGTGAAAGAAAGCGTGACTGGCGGTTAAATTGCCAACGCGTATGACCCAGCTCGAGGCAAAAATCCATTTCGCTGGTGGTGAGATGCGGTAGGGCGTGGGACGGGGCGGGGAGGGTTACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTAcgc                                                                                                                                                                                                            <  1:269349/249‑1 (MQ=255)
    cgcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTAc                                                                                                                                                                                                        <  2:145487/251‑1 (MQ=255)
     gcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAATGCCAACCCCTATTAACCAGATCGATGCAAAAATTCCTTTTGCTGGGGGGCAGATGCGGGGTGGGGTGGGACGCGGCGGGGGGCGGCACACTGAGGGTTTTCCCCCA‑ACGCCACCGCTTCCCAGCGCTGGTGTGCCCCACCTCTTACCATGCGGTCGCGGTCgggttcccttagcctacc                                                                                                                                                                                                       >  1:104564/1‑235 (MQ=255)
          gtcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTCCGCGTACTGTGAg                                                                                                                                                                                                  <  2:159108/249‑1 (MQ=255)
                   atttgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTg                                                                                                                                                                                         <  1:235106/249‑1 (MQ=255)
                          ttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGgc                                                                                                                                                                                   <  1:64842/250‑1 (MQ=255)
                                    aaTCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCGCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGGCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCt                                                                                                                                                                        <  2:53306/251‑1 (MQ=255)
                                               gtgAAAGAAAGCCTTGCTGGCGGTTAAATTGCCAAAGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCGGGTTGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGAGGTAGTTc                                                                                                                                                              <  2:146966/250‑1 (MQ=255)
                                                    agaaagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAg                                                                                                                                                        <  1:141572/251‑1 (MQ=255)
                                                       aagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑CTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGAGGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCGGCGGTAGTTCAGGCAGTtc                                                                                                                                                     <  1:38534/251‑1 (MQ=255)
                                                                      ttAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCGCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TAc                                                                                                                                      <  1:73748/250‑1 (MQ=255)
                                                                      ttAAATTGCCAACGCATATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGGGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACCCTGAGG‑GTTTTCGCCCGACGCCACTGATGCCAAGGGCTGATGTTCCCAACTTCTGACCCTGCGGTCCCGTTCTGTTGCACCACGCGTAATGTGAGCCTGGGTTGACCGGCGCTCTCCGCCAGCGGGTGTTCGGGCCGTTCAATCGACTGtt‑ttac                                                                                                                                     >  2:150838/1‑248 (MQ=255)
                                                                        aaaTTGCCAACGCTGATTACCCAGCTCGATGCAAAAAGCCATTGCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTTAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTCACCATGCGGTCTCGTTCGGTTCCACTACGCGGACTGTGAGCCAGAGTTGCCCGTCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTG‑TACCt                                                                                                                                    <  2:260607/250‑1 (MQ=255)
                                                                             gCCAACGCTTATTACCCAGCTTGATGCAAAAATCCAGTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTGCCGCCAGACGCCACTTCTGCCAGGCGCTGATGTTCCCAACTGCTGACCATGCGGTCGCGTTCGGTTGCACTACGCTTACTGTGAGCCAGAGTTGCCCTGCGCGCTCCGGCGGCGGTAGGTCAGCCAGTTCAATCAACTGGT‑TACCGTGTgg                                                                                                                               <  2:261817/250‑1 (MQ=255)
                                                                              ccAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAg                                                                                                                             <  2:62955/251‑1 (MQ=255)
                                                                               cAACGCTTATTACCCAGCTCGATGCAAAACTCCATTTAGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTTCCAGGCGCTGATGTGCCCAACTTCTGACCATGAGGTCGGGTTCGGTTGCACTACGCGGTCTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCCGTTCAATCAACTGTT‑TTCCTAGTggtgc                                                                                                                            >  2:44016/1‑248 (MQ=255)
                                                                                      tattaCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACTCCCCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCTTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAAGTGTT‑TACCTTGTGGAGCGACAt                                                                                                                       <  2:181658/249‑1 (MQ=255)
                                                                                       attaCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCa                                                                                                                    <  2:224680/251‑1 (MQ=255)
                                                                                        ttaCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTAAGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCCGTTCAATCAACTGTT‑TTCCTTGTGGAGCGACATcc                                                                                                                     >  2:251196/1‑249 (MQ=255)
                                                                                                      gCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCg                                                                                                      <  2:253198/250‑1 (MQ=255)
                                                                                                        aaaaaTCCATTTCGGTGGGGGTCAGATGGGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGGCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGTGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTTTGGAGCGACATCCAGAGGCACTTCACCGCtt                                                                                                   <  2:175172/251‑1 (MQ=255)
                                                                                                         aaaaTCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TTCCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTg                                                                                                  <  1:162148/251‑1 (MQ=255)
                                                                                                             tCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCCGCTCTCTCCGGCTGCGGTCGTTCAGGCAGTTCAATCCACTGTT‑TTCCTTGTTGAGCGACATCCAGAGGCACTTCACCGCTTGCCac                                                                                              >  1:59397/1‑250 (MQ=255)
                                                                                                              ccATTTCGGTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGc                                                                                             <  2:59397/251‑1 (MQ=255)
                                                                                                              ccATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGTACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAg                                                                                              <  1:126863/250‑1 (MQ=255)
                                                                                                                    cGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCCCACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCCGTTCAATCAACTGTT‑TAACTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCtt                                                                                        >  1:166060/1‑250 (MQ=255)
                                                                                                                        ggtggtCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCTCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTAcc                                                                                     <  2:118971/249‑1 (MQ=255)
                                                                                                                                         tGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATcc                                                                    >  1:182038/1‑249 (MQ=255)
                                                                                                                                                  aCGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACGTCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGc                                                         <  1:26478/251‑1 (MQ=255)
                                                                                                                                                   cgcgGCGGGGAGCGTCACACTGAGG‑GTTTCCGCCAGACGCACCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGGCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACAGCCAGAGGCACTGCACCGCTTGCCAGCGGCTTACCATCCACCGCCTCCATCCAGTGCAGGAGCt                                                        <  2:166060/251‑1 (MQ=255)
                                                                                                                                                      ggcggGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGt                                                     >  1:9179/1‑251 (MQ=255)
                                                                                                                                                         ggggAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTAt                                                  <  1:197392/251‑1 (MQ=255)
                                                                                                                                                               cGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTAt                                            <  2:182038/251‑1 (MQ=255)
                                                                                                                                                               cGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTAt                                            >  1:210920/1‑251 (MQ=255)
                                                                                                                                                                                tCCGCCCGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTAtt                             <  1:154836/250‑1 (MQ=255)
                                                                                                                                                                                   gCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCt                         <  1:19773/251‑1 (MQ=255)
                                                                                                                                                                                             aCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCACGTTCGGTTGCACTACGCGTACTGTGAGCCCGAGTTGCCCGGCGCTCTCCGGCTGCGGTCGTTCAGGCGGTTCGATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGGGCTCGTTATCGCTATGACGGAACAGGCATTCGCTGGTCACTTCg               >  2:257516/1‑251 (MQ=255)
                                                                                                                                                                                                gctgcCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATg            <  1:251196/251‑1 (MQ=255)
                                                                                                                                                                                                         cgcTGATGTGCACAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAAC‑GTTGTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGcccc   <  2:148251/251‑2 (MQ=255)
                                                                                                                                                                                                          gcTGATGTGCCCCACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCgg  >  2:230401/1‑251 (MQ=255)
                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                              
ATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGG‑TTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTT‑TACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGG  >  NC_000913/363957‑364408

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