Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 832,662 T→C M97V (ATG→GTG)  ybiA ← DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913832,6620TC100.0% 216.7 / NA 61M97V (ATG→GTG) ybiADUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/32);  total (29/32)

GTCGATTCTTGCCCTTACCATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTT  >  NC_000913/832518‑832808
                                                                                                                                                |                                                                                                                                                  
gTCGATTCTTGCCCTTACCATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTg                                                                                                                                              <  2:190347/151‑1 (MQ=255)
  cGATTCTTGCCCTTACCATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGtt                                                                                                                                            <  1:365849/151‑1 (MQ=255)
   gATTCTTGCCCTTACCATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTc                                                                                                                                           >  2:184712/1‑151 (MQ=255)
    aTTCTTGCCCTTACCATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCt                                                                                                                                          >  1:491163/1‑151 (MQ=255)
    aTTCTTGCCCTTACCATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCt                                                                                                                                          >  2:389525/1‑151 (MQ=255)
                  cATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCa                                                                                                                            <  2:156749/151‑1 (MQ=255)
                  cATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCa                                                                                                                            <  2:362347/151‑1 (MQ=255)
                      accaccGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGttt                                                                                                                        >  2:753813/1‑151 (MQ=255)
                        caccGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGttttt                                                                                                                      <  2:499737/151‑1 (MQ=255)
                          ccGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTAc                                                                                                                    >  1:200130/1‑151 (MQ=255)
                                ccccAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCagag                                                                                                              <  2:125252/151‑1 (MQ=255)
                                ccccAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCa                                                                                                                 <  1:758608/148‑1 (MQ=255)
                                 cccAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAgg                                                                                                             <  2:748619/151‑1 (MQ=255)
                                  ccAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGc                                                                                                            >  1:191295/1‑151 (MQ=255)
                                     gTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTa                                                                                                         <  1:584451/151‑1 (MQ=255)
                                         gCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAAc                                                                                                     >  1:153859/1‑151 (MQ=255)
                                          cATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACg                                                                                                    <  2:162298/151‑1 (MQ=255)
                                          cATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACg                                                                                                    <  1:158857/151‑1 (MQ=255)
                                             cGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATc                                                                                                 >  2:58072/1‑151 (MQ=255)
                                              gTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATcg                                                                                                <  2:217493/151‑1 (MQ=255)
                                                                 aGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTgcggcg                                                                              >  2:655831/1‑150 (MQ=255)
                                                                         gggCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAggg                                                                     >  1:105608/1‑151 (MQ=255)
                                                                           gcgcGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAgggg                                                                    >  1:521304/1‑150 (MQ=255)
                                                                             gcgGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGga                                                                   >  1:209211/1‑149 (MQ=255)
                                                                               gggtGCCAGCAAGAGCGCACGCAGTTCTGCATTCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGCGCGGCGACCCTAGGGgaaga                                                                >  1:22327/4‑150 (MQ=255)
                                                                                   gCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAAcc                                                            <  1:316901/150‑1 (MQ=255)
                                                                                           gagCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGAtc                                                   <  1:701257/151‑1 (MQ=255)
                                                                                            agCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGAtc                                                   >  1:78760/1‑150 (MQ=255)
                                                                                               gcacgcaGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCtt                                               <  2:748476/151‑1 (MQ=255)
                                                                                                  cgcaGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTcgc                                            >  2:345825/1‑151 (MQ=255)
                                                                                                  cgcaGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTcgc                                            <  1:58041/151‑1 (MQ=255)
                                                                                                  cgcaGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTcgc                                            <  1:185783/151‑1 (MQ=255)
                                                                                                   gcaGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTcgcg                                           <  1:188338/151‑1 (MQ=255)
                                                                                                     aGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGt                                         >  2:368187/1‑151 (MQ=255)
                                                                                                      gTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGt                                         >  2:472343/1‑150 (MQ=255)
                                                                                                           gCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTAtt                                      <  2:563241/148‑1 (MQ=255)
                                                                                                           gCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTAtt                                      <  1:56093/148‑1 (MQ=255)
                                                                                                            cATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtc                                  >  1:380530/1‑151 (MQ=255)
                                                                                                             aTGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtc                                  >  2:462765/1‑150 (MQ=255)
                                                                                                               gCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtcatc                               <  1:32537/151‑1 (MQ=255)
                                                                                                               gCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTtcat                                <  2:463968/150‑1 (MQ=255)
                                                                                                                  gTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCaa                             <  2:758196/150‑1 (MQ=255)
                                                                                                                  gTTCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAg                            <  2:744488/151‑1 (MQ=255)
                                                                                                                    tCGAATTTAGGGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCGAGTTTTTACGCAGAGGCTTAGAACGATCGCGACCCATGCGTGCGGCGCCCATATGGGAAGAAACCCGACGGATCCATTCGCGGTATTTTTCAGCaa                             >  1:641078/1‑148 (MQ=255)
                                                                                                                    tCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAgg                           <  2:162119/150‑1 (MQ=255)
                                                                                                                    tCGAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAgg                           >  1:573091/1‑150 (MQ=255)
                                                                                                                      gAATTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAt                        >  2:535944/1‑151 (MQ=255)
                                                                                                                       aaTTTAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAtt                       >  1:567480/1‑151 (MQ=255)
                                                                                                                           tAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAAttttt                    <  2:528947/150‑1 (MQ=255)
                                                                                                                           tAGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTtg                   >  2:630950/1‑151 (MQ=255)
                                                                                                                            aGCGCGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTtgt                  >  2:169707/1‑151 (MQ=255)
                                                                                                                             gcgcGAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTtgtg                 <  2:699758/151‑1 (MQ=255)
                                                                                                                                 gAAGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTg             <  2:713483/151‑1 (MQ=255)
                                                                                                                                   aGTGCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGaa           <  1:753738/151‑1 (MQ=255)
                                                                                                                                      gCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAAt        >  2:283643/1‑151 (MQ=255)
                                                                                                                                      gCTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAAt        >  1:370986/1‑151 (MQ=255)
                                                                                                                                       cTTTTCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAAt        >  2:172920/1‑150 (MQ=255)
                                                                                                                                         tttCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATa       <  1:244436/149‑1 (MQ=255)
                                                                                                                                         tttCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAgt     >  1:499979/1‑151 (MQ=255)
                                                                                                                                          ttCGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAgt     <  2:545190/150‑1 (MQ=255)
                                                                                                                                            cGCACCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGtt  <  2:159292/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GTCGATTCTTGCCCTTACCATGACCACCGTCTCCCCAGTAAGCATCGTTTTCCGTATGCTCAACCAGTTTTGCGGGCGCGGTTGCCAGCAAGAGCGCACGCAGTTCTGCATGCTGTTCGAATTTAGCGCGAAGTGCTTTTCGCATCACCTGTTCTTTGACCGACTCCCAGTTTTTACGCAGAGGCTTAGAACGATCGCGCCCCATGCGTGCGGCGACCATAGGGGAAGAAACCCGACGGATCTCTTCGCGGTATTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTT  >  NC_000913/832518‑832808

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