Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,587,370 C→T T143I (ACT→ATT) 
*248* (TAG→TAA) 		 yhhA →
ugpQ ←		 DUF2756 family protein
glycerophosphodiester phosphodiesterase, cytosolic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,3700CT100.0% 228.4 / NA 71T143I (ACT→ATT) 
*248* (TAG→TAA) 		yhhA
ugpQ		DUF2756 family protein
glycerophosphodiester phosphodiesterase, cytosolic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/36);  total (35/36)

ACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCT  >  NC_000913/3587241‑3587502
                                                                                                                                 |                                                                                                                                    
aCGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGcc                                                                                                                             >  1:745467/1‑139 (MQ=255)
aCGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGcc                                                                                                                             >  1:1005983/1‑139 (MQ=255)
   gcggGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          <  2:1750766/139‑1 (MQ=255)
        aTGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGtt                                                                                                                     >  2:985275/1‑139 (MQ=255)
            tAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                 <  2:1366170/139‑1 (MQ=255)
            tAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                 <  2:1005983/139‑1 (MQ=255)
               aCAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCa                                                                                                              <  1:431165/139‑1 (MQ=255)
               aCAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCa                                                                                                              <  2:745467/139‑1 (MQ=255)
                  gCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  <  2:1481605/116‑1 (MQ=255)
                  gCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  >  1:1481605/1‑116 (MQ=255)
                  gCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATc                                                                                                           >  2:771526/1‑139 (MQ=255)
                   cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          >  2:1422637/1‑123 (MQ=255)
                   cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCAACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          <  1:1422637/123‑1 (MQ=255)
                     aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACg                                                                                                        >  1:1670150/1‑139 (MQ=255)
                     aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACg                                                                                                        >  1:791681/1‑139 (MQ=255)
                         tAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCaa                                                                                                    >  2:1792448/1‑139 (MQ=255)
                          aaTCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  2:1670150/139‑1 (MQ=255)
                            tCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCg                                                                                                 <  1:771526/139‑1 (MQ=255)
                               ggATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTTAAAATTATTGGGCCGGAAAGTTCGGACCAATCACGTCAATCGCAt                                                                                              <  1:1037424/139‑1 (MQ=255)
                                gATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  2:180420/133‑1 (MQ=255)
                                gATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   >  1:180420/1‑133 (MQ=255)
                                gATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATc                                                                                             >  1:262419/1‑139 (MQ=255)
                                 aTAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCg                                                                                            >  1:245345/1‑139 (MQ=255)
                                     ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTAc                                                                                        <  2:1386953/139‑1 (MQ=255)
                                     ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACAT‑‑CGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTAc                                                                                        <  2:913470/137‑1 (MQ=255)
                                     ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACAT‑‑CGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTAc                                                                                        >  1:913470/1‑137 (MQ=255)
                                       cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          >  1:706593/1‑103 (MQ=255)
                                       cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          <  2:706593/103‑1 (MQ=255)
                                       cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  2:1671813/126‑1 (MQ=255)
                                       cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   >  1:1671813/1‑126 (MQ=255)
                                       cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAg                                                                                      >  2:292418/1‑139 (MQ=255)
                                          tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCa                                                                                               <  1:231423/127‑1 (MQ=255)
                                          tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                   <  1:162560/139‑1 (MQ=255)
                                          tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGACGTAAAGTTCGGACCAATCACGTCAATCGCa                                                                                               >  2:231423/1‑127 (MQ=255)
                                               cagcagCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCCCGTCAAt                                                                                                   >  1:380664/1‑118 (MQ=255)
                                               cagcagCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  2:380664/118‑1 (MQ=255)
                                                       tatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAcccc                                                                      <  1:2904/139‑1 (MQ=255)
                                                            tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  >  1:177576/1‑74 (MQ=255)
                                                            tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  <  2:177576/74‑1 (MQ=255)
                                                            tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                 >  2:1141391/1‑139 (MQ=255)
                                                                 gagaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                            >  1:679092/1‑139 (MQ=255)
                                                                 gagaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                            <  2:679092/139‑1 (MQ=255)
                                                                     ggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAActct                                                        >  1:1297165/1‑139 (MQ=255)
                                                                             ggCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                                <  1:985275/139‑1 (MQ=255)
                                                                             ggCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                                <  1:1141391/139‑1 (MQ=255)
                                                                              gCGACATGTTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  >  1:542600/1‑56 (MQ=255)
                                                                              gCGACATGTTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  <  2:542600/56‑1 (MQ=255)
                                                                               cGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCt                                              <  1:292418/139‑1 (MQ=255)
                                                                                gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTg                                             <  2:245345/139‑1 (MQ=255)
                                                                                 aCATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTgg                                            >  1:979537/1‑139 (MQ=255)
                                                                                           tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  2:262419/139‑1 (MQ=255)
                                                                                               cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  2:970081/70‑1 (MQ=255)
                                                                                               cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   >  1:970081/1‑70 (MQ=255)
                                                                                               cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  >  1:491584/1‑135 (MQ=255)
                                                                                               cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  2:491584/135‑1 (MQ=255)
                                                                                                ccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGta                             >  1:205790/1‑139 (MQ=255)
                                                                                                ccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGta                             >  1:523669/1‑139 (MQ=255)
                                                                                                   gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  >  2:664831/1‑131 (MQ=255)
                                                                                                   gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  1:664831/131‑1 (MQ=255)
                                                                                                         cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGa                    >  1:543970/1‑139 (MQ=255)
                                                                                                            aGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATACAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATc                 >  1:588072/1‑139 (MQ=255)
                                                                                                              ccTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATc                                                                                                           >  1:242614/1‑47 (MQ=255)
                                                                                                              ccTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATc                                                                                                           <  2:242614/47‑1 (MQ=255)
                                                                                                              ccTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCg               >  1:141883/1‑139 (MQ=255)
                                                                                                              ccTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCg               <  2:791681/139‑1 (MQ=255)
                                                                                                                   cATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCg               >  1:1638541/1‑134 (MQ=255)
                                                                                                                   cATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCg               <  2:1638541/134‑1 (MQ=255)
                                                                                                                      tCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg       <  2:931046/139‑1 (MQ=255)
                                                                                                                      tCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg       <  2:1297165/139‑1 (MQ=255)
                                                                                                                        cGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   >  1:735665/1‑45 (MQ=255)
                                                                                                                        cGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  2:735665/45‑1 (MQ=255)
                                                                                                                        cGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg     >  1:95938/1‑139 (MQ=255)
                                                                                                                           tGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCt  <  2:141883/139‑1 (MQ=255)
                                                                                                                            gAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGc                                                                                                >  2:1154158/1‑44 (MQ=255)
                                                                                                                            gAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGc                                                                                                <  1:1154158/44‑1 (MQ=255)
                                                                                                                              aaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   >  1:994356/1‑39 (MQ=255)
                                                                                                                              aaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  2:994356/39‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
ACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCT  >  NC_000913/3587241‑3587502

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