Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,759,478 C→T 5.9% S349N (AGT→AAT)  sufS ← cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,759,4780CT 5.9% 238.1 / 4.0 68S349N (AGT→AAT) sufScysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (2/2);  total (34/34)
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.98e-01

AGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGG  >  NC_000913/1759343‑1759616
                                                                                                                                       |                                                                                                                                          
aGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGcc                                                                                                                                        <  1:635752/140‑1 (MQ=255)
  tggATCCACTTCTTCATGGTTGTTTTACATGGCCAGCGACGCCCGACACCTCGCAGGGACGTGGTTCTAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCATTGCCGTAATTATCGAGATAACTGCCaa                                                                                                                                      <  2:1657685/138‑1 (MQ=255)
  tcGATACACTTCTTCATGGGTGTTTTACATGGCAAGCGACGCCCGACACATCGCAGGGACGTTTTAATAGGACATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCaa                                                                                                                                      <  2:1657677/139‑1 (MQ=255)
  aCGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCaa                                                                                                                                      <  2:1571744/140‑1 (MQ=255)
  aCGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATAGCAGGGACGTTGTAATAGGCCATCAATGTCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCaa                                                                                                                                      <  2:1657589/140‑1 (MQ=255)
     aTCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAAcat                                                                                                                                   >  2:292246/1‑140 (MQ=255)
     aTCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAAcat                                                                                                                                   >  1:1509520/1‑140 (MQ=255)
       ctaCTTGTTCTTGGGTGTTATACATGGCCAGCGACGCCCGACTCATAGCAGGGACGTTGTATTAGGCCATCAATGGCAGTGCGCAGTGATGTCCGGTAGGCACAGCACGGCCGTAATTATCGAGAAAACTGCCAAcatca                                                                                                                                 <  2:1681848/138‑1 (MQ=255)
            tcttcATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGc                                                                                                                            <  1:1094751/140‑1 (MQ=255)
              ttcATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCgt                                                                                                                          <  1:1177127/140‑1 (MQ=255)
              ttcATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCAGTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCgt                                                                                                                          <  2:397983/140‑1 (MQ=255)
                     gtgtTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGttt                                                                                                                   >  1:854077/1‑140 (MQ=255)
                        ttATACATGGCCAGCGACGCCCGACACATCGCAGGGACTTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTAcc                                                                                                                >  1:834338/1‑140 (MQ=255)
                              aTGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAGTGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGAtt                                                                                                          >  1:1586870/1‑140 (MQ=255)
                              aTGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGAtt                                                                                                          >  1:722341/1‑140 (MQ=255)
                                ggCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTaa                                                                                                        >  1:1394404/1‑140 (MQ=255)
                                ggCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTaa                                                                                                        >  2:1482473/1‑140 (MQ=255)
                                ggCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTaa                                                                                                        >  1:1527691/1‑140 (MQ=255)
                                ggCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTaa                                                                                                        >  2:805816/1‑140 (MQ=255)
                                ggCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTaa                                                                                                        >  2:586722/1‑140 (MQ=255)
                                  ccAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTaaaa                                                                                                      >  1:1412947/1‑140 (MQ=255)
                                                        aGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATTGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGtttt                                                                                <  1:1286610/140‑1 (MQ=255)
                                                        aGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGtttt                                                                                <  1:829817/140‑1 (MQ=255)
                                                          ggACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTtgt                                                                              <  2:324740/140‑1 (MQ=255)
                                                          ggACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTtgt                                                                              <  1:1021592/140‑1 (MQ=255)
                                                          ggACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTtgt                                                                              <  1:1277941/140‑1 (MQ=255)
                                                          ggACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTtgt                                                                              <  1:1633083/140‑1 (MQ=255)
                                                          ggACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTtgt                                                                              <  1:310200/140‑1 (MQ=255)
                                                          ggACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTtgt                                                                              <  1:1633187/140‑1 (MQ=255)
                                                               ttGTAATAGGCCATCAATGTCATAGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACAGACAACATCAAAGGCGTGGTGTTTAAAGAGAAAAAAAGCAATAACGACAAGACTGTAATGAGGGcc                                                                         >  2:795968/1‑140 (MQ=255)
                                                                tGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCa                                                                        >  2:1198134/1‑140 (MQ=255)
                                                                    aTAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATaga                                                                    <  1:410888/140‑1 (MQ=255)
                                                                        gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGt                                                                >  1:476396/1‑140 (MQ=255)
                                                                         ccATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTg                                                               <  1:978734/140‑1 (MQ=255)
                                                                         ccATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTg                                                               <  1:586625/140‑1 (MQ=255)
                                                                            tCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTgaga                                                            >  1:825156/1‑140 (MQ=255)
                                                                            tCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTgaga                                                            >  1:739581/1‑140 (MQ=255)
                                                                            tCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTgaga                                                            >  1:1078189/1‑140 (MQ=255)
                                                                            tCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTgaga                                                            >  2:914186/1‑140 (MQ=255)
                                                                            tCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACAGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTgaga                                                            >  1:1106776/1‑140 (MQ=255)
                                                                              aaTGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATc                                                          <  2:1580091/140‑1 (MQ=255)
                                                                              aaTGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATc                                                          <  2:1580118/140‑1 (MQ=255)
                                                                              aaTGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATc                                                          <  2:833231/140‑1 (MQ=255)
                                                                                   cATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAATTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTa                                                     >  2:1331305/1‑140 (MQ=255)
                                                                                   cATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAATTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTa                                                     >  2:1332245/1‑140 (MQ=255)
                                                                                        cgcAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCTAGAAAACTGCCAACATCATTGGCGTGGTGTTTAGCGAGATTAAAAGCAATAACGACAAGCCTGTTTTGTGGGCCTTATAGAGTGAGATCCGGTACAGAt                                                >  2:633385/1‑140 (MQ=255)
                                                                                           aGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTcc                                             <  1:1137652/140‑1 (MQ=255)
                                                                                           aGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTcc                                             <  1:1137541/140‑1 (MQ=255)
                                                                                           aGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTcc                                             <  2:547177/140‑1 (MQ=255)
                                                                                                  tCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCtgtg                                      >  1:997247/1‑140 (MQ=255)
                                                                                                  tCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCtgtg                                      >  1:362997/1‑140 (MQ=255)
                                                                                                  tCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTAGTGGGCCATAGAGAGTGAGATCCGGAACAGATTCCAGCtgag                                      >  1:613005/1‑138 (MQ=255)
                                                                                                        aCGCACAGCAATGCCGTTATTATCGAGAAAACTGCCAACATCATAGGAGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCATCTGTGATAgcg                                <  2:1625211/140‑1 (MQ=255)
                                                                                                              aGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAAt                          <  2:533478/140‑1 (MQ=255)
                                                                                                                   tGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATc                     >  1:1536177/1‑140 (MQ=255)
                                                                                                                       gTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGAt                 <  1:1028101/140‑1 (MQ=255)
                                                                                                                        tAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGAtt                <  1:939887/140‑1 (MQ=255)
                                                                                                                        tAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGAtt                >  2:967224/1‑140 (MQ=255)
                                                                                                                        tAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGACTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGAtt                <  1:940707/140‑1 (MQ=255)
                                                                                                                          aTTATCGAGAAAACTGCCAACATCCTAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCt              >  2:977687/1‑140 (MQ=255)
                                                                                                                          aTTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCACGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCt              >  2:111506/1‑140 (MQ=255)
                                                                                                                          aTTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCt              >  1:972093/1‑140 (MQ=255)
                                                                                                                            tATCGAGAAAATTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGt            <  1:1332245/140‑1 (MQ=255)
                                                                                                                            tATCGAGAAAATTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGt            <  1:1331305/140‑1 (MQ=255)
                                                                                                                               cGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCa         >  2:1127889/1‑140 (MQ=255)
                                                                                                                                  gaAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATAc      >  1:166191/1‑140 (MQ=255)
                                                                                                                                  gaAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATAc      >  1:162879/1‑140 (MQ=255)
                                                                                                                                      aCTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCgg  <  1:24488/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
AGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGG  >  NC_000913/1759343‑1759616

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: