Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,529,856 T→C 5.5% K107K (AAA→AAG)  ligA ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,529,8560TC 5.5% 198.3 / 2.3 75K107K (AAA→AAG) ligADNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (41/28);  new base C (2/2);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATGGCGCGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACG  >  NC_000913/2529712‑2529996
                                                                                                                                                |                                                                                                                                            
aaaggCGCGAATAGTACGGACATTCGACGTGATATCTTCCCCGGTGGTGCCATCGCAACGGGTCGCGGCACTGACTAAACCGCCATTTTCATACATAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGAGACTTtct                                                                                                                                          <  2:935302/146‑1 (MQ=255)
    gcgcGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCCTCTCAACGGGTCGCGGCACTGACTAAAACGCCAATTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCgtt                                                                                                                                      >  1:710487/1‑149 (MQ=255)
    gcgcGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCgtt                                                                                                                                      <  2:352567/149‑1 (MQ=255)
      gcgAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCgttgt                                                                                                                                    >  1:625933/1‑149 (MQ=255)
        gattAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACCTTCTCGTTGttt                                                                                                                                  >  1:779175/4‑149 (MQ=255)
        gAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACCTTCTCGTTGttt                                                                                                                                  >  1:779173/1‑149 (MQ=255)
         aaTAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGCCGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGtttt                                                                                                                                 >  2:371401/1‑149 (MQ=255)
          aTAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCTCTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGttttt                                                                                                                                >  1:413441/1‑149 (MQ=255)
                 gCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACgg                                                                                                                         <  2:161178/149‑1 (MQ=255)
                   acaTTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTc                                                                                                                       >  1:190864/1‑149 (MQ=255)
                   acaTTAGACGTGATATCTTCCACGGTGGTGCCATCGCCACGGTTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTc                                                                                                                       >  1:611491/1‑149 (MQ=255)
                     aTTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCt                                                                                                                     <  1:331429/149‑1 (MQ=255)
                         gACGTGATATCTTCCCCGGTGGAGCCATCGCCACGGGTCGCGGCACTGACTAAAACTCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGCCGGTCCTGcac                                                                                                                 >  1:534435/1‑149 (MQ=255)
                          aCGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCCGCTTCAGCTCACAGCACCAGGTGAATTTCTCGTTGTTTTTCAGCCGGTCCTGcacc                                                                                                                >  2:919820/1‑148 (MQ=255)
                            gTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACCGCACCAGGTGACTATCTCGTTCTTTTTCAGACGGTCCTGcacaag                                                                                                              >  2:254534/1‑147 (MQ=255)
                              gATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGtt                                                                                                            >  2:395524/1‑149 (MQ=255)
                                 atCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTttgt                                                                                                         >  2:64654/1‑149 (MQ=255)
                                  tCTTCCCCGGGGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCAAACAGAATCCTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTATTTCAGACGGTCCTGCACACGTttgtt                                                                                                        <  2:563897/149‑1 (MQ=255)
                                    ttCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCCCCAGGTGACTGTCTCGTTGTTTTTCAGCCGGTCCTGCACACGTTTGTTGa                                                                                                      >  1:463718/1‑149 (MQ=255)
                                    ttCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGa                                                                                                      >  1:271153/1‑149 (MQ=255)
                                        ccGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGc                                                                                                  >  2:637490/1‑149 (MQ=255)
                                           gtggtgCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGaaagaaag                                                                                               >  2:690765/1‑149 (MQ=255)
                                               tgCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAg                                                                                           <  1:168525/149‑1 (MQ=255)
                                                gCCATCGCCACGGGTCGCGGCACTGACTAAACCGCCATTTTCATACAGAATACTGACGGCAAGCCCATCCAGCTTCAGCTCACAGCCCCAGGTGACTTTCTCGGTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGc                                                                                          >  1:441690/1‑149 (MQ=255)
                                                gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGATTTTCAGACGGTCCTGCACACGTTTGTAGAAAGCAAGaaatc                                                                                          >  1:430386/1‑147 (MQ=255)
                                                 ccATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTGTGTTGAAAGCAAGAAAGCt                                                                                         >  1:52057/1‑149 (MQ=255)
                                                  cATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCtt                                                                                        <  1:280293/149‑1 (MQ=255)
                                                   aTCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCttt                                                                                       >  1:858082/1‑149 (MQ=255)
                                                   aTCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCttt                                                                                       >  1:1019270/1‑149 (MQ=255)
                                                    tCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGGGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTttc                                                                                      >  2:816237/1‑149 (MQ=255)
                                                         aCGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGTCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTtcat                                                                                 <  2:561159/149‑1 (MQ=255)
                                                          cGGGTCGCGGCACTGACTAAACCGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCTCCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTtcatc                                                                                <  1:754557/149‑1 (MQ=255)
                                                                   gCACTGACTAAAAGGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGtt                                                                       <  1:931874/149‑1 (MQ=255)
                                                                      cTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATc                                                                    <  2:902307/149‑1 (MQ=255)
                                                                       tGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCCGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGATTTTCAGACGGTCCTGCACACTTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGATATcc                                                                   >  2:338244/1‑149 (MQ=255)
                                                                       tGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCCCCAGGTGACGATCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATcc                                                                   >  2:945486/1‑149 (MQ=255)
                                                                         aCTAAAACGCCATTTTCATACAGAATACTGACGGCAAGCCCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAATGCTTTCTTCATCAAAAACGTTATCCAg                                                                 >  2:215598/1‑149 (MQ=255)
                                                                         aCTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGATAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAg                                                                 <  2:833599/149‑1 (MQ=255)
                                                                         aCTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAg                                                                 >  1:282000/1‑149 (MQ=255)
                                                                          cTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGt                                                                <  2:451239/149‑1 (MQ=255)
                                                                             aaaCGCCATTTTCATACAGAATACTGACGGCAAGACAATCCAGCTTCAGCTAACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTACTGCACACGTTTGTTGAAAGAAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGaa                                                             >  2:142445/1‑148 (MQ=255)
                                                                               aCGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAg                                                           <  1:690765/149‑1 (MQ=255)
                                                                               aCGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAg                                                           <  2:271153/149‑1 (MQ=255)
                                                                                agCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCATGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGc                                                          <  1:1022194/148‑1 (MQ=255)
                                                                                  gcATTTTCATACAGAATACTGACGGCATGACCATCCAGCTTCAGCTCACAGCACCAGGTGTCTTTCTCGTTGTTTTTCAGACGGTCCTTCACACGTTTTTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCAt                                                        <  2:250225/148‑1 (MQ=255)
                                                                                       ttCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACATCATTGGTa                                                   >  1:645639/1‑149 (MQ=255)
                                                                                        tcctACAGAATACAGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGCTGACATTCTCGTGGTTTTTCAGACGGACCTGCACACGTTTGTTTAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTAc                                                  <  1:574325/146‑1 (MQ=255)
                                                                                        tCATACTGTATACTGACTGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGTAAGCTAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTAc                                                  <  1:556641/149‑1 (MQ=255)
                                                                                            aCAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCa                                              <  1:847979/149‑1 (MQ=255)
                                                                                              aGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATg                                            <  1:216801/149‑1 (MQ=255)
                                                                                              aGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATg                                            <  2:645341/149‑1 (MQ=255)
                                                                                              aGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACCTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATg                                            <  2:779175/149‑1 (MQ=255)
                                                                                              aGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACCTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATg                                            <  2:779173/149‑1 (MQ=255)
                                                                                                 aTACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACAGGTTTTTTGAAAGCAAGAAAGCTTTCTTCATCAATAACGTTATCCAGTGACAGCATTGGTACTTCATGGCg                                         <  2:534435/149‑1 (MQ=255)
                                                                                                    cTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTAt                                      >  2:366526/1‑149 (MQ=255)
                                                                                                    cTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTGTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGAGTAt                                      >  2:772632/1‑149 (MQ=255)
                                                                                                      gACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGCCTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCt                                    >  2:828748/1‑149 (MQ=255)
                                                                                                      gACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCt                                    >  1:836317/1‑149 (MQ=255)
                                                                                                       aCGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTg                                   >  1:986129/1‑149 (MQ=255)
                                                                                                        cGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCTAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTgg                                  <  1:254534/149‑1 (MQ=255)
                                                                                                        cGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTgg                                  <  2:190864/149‑1 (MQ=255)
                                                                                                         ggCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGc                                 <  1:371401/149‑1 (MQ=255)
                                                                                                         ggCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGc                                 >  2:426933/1‑149 (MQ=255)
                                                                                                          gCAAGACCATCCAGCTTCAGCTCACAGCACCAGGGGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGTAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCt                                <  2:161763/149‑1 (MQ=255)
                                                                                                                 cATCCAGCTTCAGCTCACAGCACCAGGGGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAgcc                         <  2:463718/149‑1 (MQ=255)
                                                                                                                    ccAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAgccgcc                      >  1:847967/1‑149 (MQ=255)
                                                                                                                        cTTCAGCTCACAGCACCAGGTGACTTACTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAATCAAGAATGCTTTCTTCATCAAAAACGATATCCACTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAgcg                  <  1:477152/149‑1 (MQ=255)
                                                                                                                         ttCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAgcgg                 >  1:558356/1‑149 (MQ=255)
                                                                                                                         ttCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAgcgg                 >  2:711522/1‑149 (MQ=255)
                                                                                                                          tCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCCgcggc                >  1:566190/1‑149 (MQ=255)
                                                                                                                             gCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGCCAGCATTGGTACTTCATGGCGTATCTGGCTGAACGCCGCCAGCGGCGCa             >  2:1016541/1‑149 (MQ=255)
                                                                                                                                 acaGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc         >  2:136633/1‑149 (MQ=255)
                                                                                                                                 acaGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc         >  1:974321/1‑149 (MQ=255)
                                                                                                                                      aCCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTaca    >  1:290297/1‑149 (MQ=255)
                                                                                                                                        cAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACg  >  2:301758/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                            
AATGGCGCGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGCTTCAGCTCACAGCACCAGGTGACTTTCTCGTTGTTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACG  >  NC_000913/2529712‑2529996

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