Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,557,286 G→A 6.5% T265M (ACG→ATG)  rtcB ← RNA‑splicing ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,557,2860GA 6.5% 190.2 / 2.5 61T265M (ACG→ATG) rtcBRNA‑splicing ligase
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

AATTCCATATTGACCAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCG  >  NC_000913/3557147‑3557433
                                                                                                                                           |                                                                                                                                                   
aatTCCATATTGACCAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACTGTTTTCTGCGt                                                                                                                                           <  1:243034/150‑1 (MQ=255)
         ttGACCAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCATTTTCTGCGTAATGCTCTg                                                                                                                                  >  1:491139/1‑150 (MQ=255)
         ttGACCAGCACGCGCAGACACCGCCCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTg                                                                                                                                  >  2:621952/1‑150 (MQ=255)
           gACCAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCa                                                                                                                                >  2:970485/1‑150 (MQ=255)
           gACCAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCa                                                                                                                                >  1:742165/1‑150 (MQ=255)
           gACCAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCa                                                                                                                                >  1:617868/1‑150 (MQ=255)
              cAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATg                                                                                                                             >  1:1123264/1‑150 (MQ=255)
                 cACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCg                                                                                                                          >  2:1224607/1‑150 (MQ=255)
                 cACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCg                                                                                                                          >  1:536282/1‑150 (MQ=255)
                    gcgcAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTa                                                                                                                       >  2:742061/1‑150 (MQ=255)
                      gcAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAcc                                                                                                                     >  1:451477/1‑150 (MQ=255)
                           cacaGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGAGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGtt                                                                                                                >  2:1124038/1‑150 (MQ=255)
                            acCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGttt                                                                                                               <  2:536282/150‑1 (MQ=255)
                               gcgcCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc                                                                                                            <  2:219496/150‑1 (MQ=255)
                                  ccTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCcatc                                                                                                         >  2:1334427/1‑150 (MQ=255)
                                   cTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCcatca                                                                                                        <  2:358401/150‑1 (MQ=255)
                                    tttACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCcatcat                                                                                                       >  1:962270/1‑150 (MQ=255)
                                      tACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATcatcg                                                                                                     <  1:958006/150‑1 (MQ=255)
                                      tACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATcatcg                                                                                                     <  1:81592/150‑1 (MQ=255)
                                      tACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCATTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATcatcg                                                                                                     <  1:1120475/150‑1 (MQ=255)
                                       aCGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATcatcgc                                                                                                    >  1:945191/1‑150 (MQ=255)
                                                 tAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTa                                                                                          <  2:1338319/150‑1 (MQ=255)
                                                      ctctTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCt                                                                                     <  2:412013/150‑1 (MQ=255)
                                                        ctTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTa                                                                                          >  1:1321851/1‑143 (MQ=255)
                                                        ctTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTa                                                                                          <  2:1321851/143‑1 (MQ=255)
                                                        ctTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTgg                                                                                   >  2:674434/1‑150 (MQ=255)
                                                        ctTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTgg                                                                                   >  2:1039025/1‑150 (MQ=255)
                                                        ctTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCATTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTgg                                                                                   >  2:757901/1‑150 (MQ=255)
                                                        ctTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCATTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTgg                                                                                   >  2:1083959/1‑150 (MQ=255)
                                                         ttCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGc                                                                                  <  2:212351/150‑1 (MQ=255)
                                                           cACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCca                                                                                                           >  1:904963/1‑123 (MQ=255)
                                                           cACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCca                                                                                                           <  2:904963/123‑1 (MQ=255)
                                                                     cTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTgcgc                                                                      <  1:923256/150‑1 (MQ=255)
                                                                        ttCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCa                                                                   <  2:954696/150‑1 (MQ=255)
                                                                         tCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAg                                                                  <  1:361265/150‑1 (MQ=255)
                                                                         tCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAg                                                                  >  1:1371827/1‑150 (MQ=255)
                                                                         tCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCATTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAg                                                                  <  1:1083959/150‑1 (MQ=255)
                                                                          cTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGc                                                                                  <  1:1079083/133‑1 (MQ=255)
                                                                          cTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGc                                                                                  >  2:1079083/1‑133 (MQ=255)
                                                                          cTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAgg                                                                 <  2:617868/150‑1 (MQ=255)
                                                                            ttttGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGcc                                                               >  1:1268927/1‑150 (MQ=255)
                                                                               tGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACg                                                            >  2:323584/1‑150 (MQ=255)
                                                                                 cacaTAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGc                                                          <  2:1302338/150‑1 (MQ=255)
                                                                                  acaTAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTGCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCt                                                         <  1:742061/150‑1 (MQ=255)
                                                                                      aGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCATTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCa                                                     <  1:757901/150‑1 (MQ=255)
                                                                                        ttGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAgg                                                   <  1:1039025/150‑1 (MQ=255)
                                                                                         tgtgGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGt                                                  <  1:621952/150‑1 (MQ=255)
                                                                                           tggtgACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCATTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTaa                                                <  2:491139/150‑1 (MQ=255)
                                                                                                  aGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCaa                                         <  2:962270/150‑1 (MQ=255)
                                                                                                   gTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCaaa                                        <  1:749726/150‑1 (MQ=255)
                                                                                                        tctcTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTAtt                                   >  1:217113/1‑150 (MQ=255)
                                                                                                         ctctTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTc                                  <  2:805841/150‑1 (MQ=255)
                                                                                                               cATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTAcc                            >  1:1277276/1‑150 (MQ=255)
                                                                                                                   gCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTcc                        <  1:1014122/150‑1 (MQ=255)
                                                                                                                       gCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATaa                    <  2:1123264/150‑1 (MQ=255)
                                                                                                                       gCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATaa                    >  2:853539/1‑150 (MQ=255)
                                                                                                                             gtgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACg              <  1:323584/150‑1 (MQ=255)
                                                                                                                             gtgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGATTTCCATCATCGAATCGCGGTTAAGCCTGGCAAAAAGATGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACg              >  2:707694/1‑150 (MQ=255)
                                                                                                                                    tCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATc       <  1:1334427/150‑1 (MQ=255)
                                                                                                                                      tGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCAc     <  1:970485/150‑1 (MQ=255)
                                                                                                                                         ccGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGCCCACGGCTTTCAGGTATTCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCAcgcg  >  2:751828/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
AATTCCATATTGACCAGCACGCGCAGACACCGCGCCTTTACGCGTCACGTAGATCTCTTCACCAAAGTGCTGTTCTTTTTGCACATAGTTGTGGTGACAGTTGATCTCTTCCATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCATCATCGCATCGCGGTTAAGGCTGGCAAAAAGCTGCGCCCAGGCCACGGCTTTCAGGTAATCATCAAAGTATTCCGTACCTTCCATAAAGTACGCCAGATCACGCG  >  NC_000913/3557147‑3557433

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