Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 344,991 G→A *272* (TAG→TAA yahL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913344,9910GA100.0% 237.1 / NA 73*272* (TAG→TAAyahLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/33);  total (40/33)

CTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGG  >  NC_000913/344860‑345125
                                                                                                                                   |                                                                                                                                      
ctctctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGaa                                                                                                                                 <  2:89923/139‑1 (MQ=255)
  ctctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAg                                                                                                                               >  1:150543/1‑139 (MQ=255)
   tctGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAgg                                                                                                                              >  2:896698/1‑139 (MQ=255)
         cTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCa                                                                                                                        <  1:190424/139‑1 (MQ=255)
          tGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCAc                                                                                                                       >  2:402549/1‑139 (MQ=255)
            gACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCt                                                                                                                     >  1:177760/1‑139 (MQ=255)
              cGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGc                                                                                                                   <  1:574316/139‑1 (MQ=255)
                               atgatgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCt                                                                                                  >  1:103165/1‑139 (MQ=255)
                               atgatgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCt                                                                                                  <  2:856657/139‑1 (MQ=255)
                                tgatgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTg                                                                                                 >  1:1019063/1‑139 (MQ=255)
                                  atgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTgag                                                                                               >  2:785398/1‑139 (MQ=255)
                                     cgcgATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGt                                                                                            >  2:106/1‑139 (MQ=255)
                                     cgcgATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGt                                                                                            <  2:177760/139‑1 (MQ=255)
                                      gcgATTAAATGAATGGATATCTAATAGAGTGCAGGAGTTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                                                                                           <  1:673367/139‑1 (MQ=255)
                                      gcgATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                                                                                           <  1:106/139‑1 (MQ=255)
                                                   tgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAg                                                                              <  1:896698/139‑1 (MQ=255)
                                                           tAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCAt                                                                      >  2:477399/1‑139 (MQ=255)
                                                              tAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTcc                                                                                                      >  1:782129/1‑104 (MQ=255)
                                                              tAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTcc                                                                                                      <  2:782129/104‑1 (MQ=255)
                                                                gagTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCat                                                                 >  1:492423/1‑139 (MQ=255)
                                                                gagTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCat                                                                 <  2:492423/139‑1 (MQ=255)
                                                                 agTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                                                                                           <  1:509764/112‑1 (MQ=255)
                                                                 agTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                                                                                           >  2:509764/1‑112 (MQ=255)
                                                                   tacAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTa                                                                                                          >  1:502510/3‑95 (MQ=255)
                                                                   tacAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTa                                                                                                          <  2:502510/93‑1 (MQ=255)
                                                                   tGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctg                                                              >  2:39281/1‑139 (MQ=255)
                                                                   tGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctg                                                              >  2:673707/1‑139 (MQ=255)
                                                                   tGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAACCGCAATGCAGGTAATCATTTCatctg                                                              >  2:682972/1‑139 (MQ=255)
                                                                    gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                                                                                           >  1:997062/1‑109 (MQ=255)
                                                                    gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                                                                                           <  2:997062/109‑1 (MQ=255)
                                                                    gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgt                                                             >  1:882779/1‑139 (MQ=255)
                                                                    gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgt                                                             >  1:920645/1‑139 (MQ=255)
                                                                    gCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgt                                                             >  2:476697/1‑139 (MQ=255)
                                                                     cAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTTAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCat                                                                 <  2:219520/134‑1 (MQ=255)
                                                                     cAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTTAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCat                                                                 >  1:219520/1‑134 (MQ=255)
                                                                     cAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgta                                                            <  1:673707/139‑1 (MQ=255)
                                                                      aGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtat                                                           <  1:476697/139‑1 (MQ=255)
                                                                             aGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgtatatctgta                                                    <  2:882779/139‑1 (MQ=255)
                                                                                 gACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGcc                                                <  1:39281/139‑1 (MQ=255)
                                                                                 gACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGcc                                                <  2:920645/139‑1 (MQ=255)
                                                                                 gACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGcc                                                <  2:150543/139‑1 (MQ=255)
                                                                                                gTATATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGc                                         <  1:696196/131‑1 (MQ=255)
                                                                                                gTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCCGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGCGAATCGCAATGCAGGTAATCATTTCATCTGTTTATCTGTATGCCAAc                                             >  2:414031/1‑127 (MQ=255)
                                                                                                gTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAAc                                             <  1:414031/127‑1 (MQ=255)
                                                                                                gTAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGc                                         >  2:696196/1‑131 (MQ=255)
                                                                                                 tAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                                                                   <  2:618883/104‑1 (MQ=255)
                                                                                                 tAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                                                                   >  1:618883/1‑104 (MQ=255)
                                                                                                 tAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgta                                                            >  1:952722/1‑111 (MQ=255)
                                                                                                 tAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCatctgta                                                            <  2:952722/111‑1 (MQ=255)
                                                                                                 tAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCAtct                                >  2:705141/1‑139 (MQ=255)
                                                                                                 tAAATAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCAtct                                >  2:209893/1‑139 (MQ=255)
                                                                                                   aaTAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAgttgt                                                                                                            >  2:555158/1‑61 (MQ=255)
                                                                                                   aaTAAGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAgttgt                                                                                                            <  1:555158/61‑1 (MQ=255)
                                                                                                      aaGCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCtttt                           >  2:378609/1‑139 (MQ=255)
                                                                                                        gCACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCtttttt                         >  1:992512/1‑139 (MQ=255)
                                                                                                         cACGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTc                        >  1:339600/1‑139 (MQ=255)
                                                                                                           cGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                                                                   <  1:327065/94‑1 (MQ=255)
                                                                                                           cGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTc                                                                   >  2:327065/1‑94 (MQ=255)
                                                                                                           cGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAAc                                             >  1:316899/1‑116 (MQ=255)
                                                                                                           cGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAAc                                             <  2:316899/116‑1 (MQ=255)
                                                                                                           cGTTAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGc                      >  1:682204/1‑139 (MQ=255)
                                                                                                             ttAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCa                                                                               >  1:1063084/1‑80 (MQ=255)
                                                                                                             ttAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCa                                                                               <  2:1063084/80‑1 (MQ=255)
                                                                                                             ttAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCg                    <  1:378609/139‑1 (MQ=255)
                                                                                                             ttAAAATCCGCTACACACTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCg                    >  1:266842/1‑139 (MQ=255)
                                                                                                              tAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGcc                                                <  2:630974/110‑1 (MQ=255)
                                                                                                              tAAAATCCGCTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGcc                                                >  1:630974/1‑110 (MQ=255)
                                                                                                                       cTACACTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGGAATCATTTCATATGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTttgtt          >  2:46879/1‑139 (MQ=255)
                                                                                                                          cacTCTTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTtgtg       <  2:1019063/139‑1 (MQ=255)
                                                                                                                            ctctTTAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGAt     <  1:705141/139‑1 (MQ=255)
                                                                                                                               tttAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCaa                                                                                   >  1:356931/1‑58 (MQ=255)
                                                                                                                               tttAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCaa                                                                                   <  2:356931/58‑1 (MQ=255)
                                                                                                                               tttAATTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCgg  >  2:523427/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
CTCTCTGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAATCATTTCATCTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGG  >  NC_000913/344860‑345125

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: