Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,483,536 T→C intergenic (‑17/‑250) pepA ← / → yjgP aminopeptidase A/I/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,483,5360TC100.0% 673.8 / NA 208intergenic (‑17/‑250)pepA/yjgPaminopeptidase A/I/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (90/118);  total (90/118)

Only 100 of 211 total aligned reads displayed.

GGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAACG  >  NC_000913/4483446‑4483625
                                                                                          |                                                                                         
ggTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCt                                                                                         <  1:2204884/93‑1 (MQ=255)
    cGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTaaa                                                                                     <  1:542889/93‑1 (MQ=255)
     gAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAg                                                                                    >  2:3732347/1‑93 (MQ=255)
      aaGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGa                                                                                   <  1:4008227/93‑1 (MQ=255)
      aaGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGa                                                                                   <  1:3380656/93‑1 (MQ=255)
       aGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGAc                                                                                  <  1:4039340/93‑1 (MQ=255)
       aGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGAc                                                                                  <  2:3855446/93‑1 (MQ=255)
       aGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGAc                                                                                  >  1:2432626/1‑93 (MQ=255)
        gACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACa                                                                                 >  1:2546623/1‑93 (MQ=255)
        gACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACa                                                                                 >  1:4481739/1‑93 (MQ=255)
        gACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACa                                                                                 >  2:1026530/1‑93 (MQ=255)
        gACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACa                                                                                 >  1:2946569/1‑93 (MQ=255)
        gACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACa                                                                                 >  1:2817028/1‑93 (MQ=255)
        gACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACa                                                                                 >  2:2801245/1‑93 (MQ=255)
          cGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAc                                                                               <  1:923238/93‑1 (MQ=255)
            cccACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcgg                                                                             <  2:1638662/93‑1 (MQ=255)
             ccACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggc                                                                            <  2:1955402/93‑1 (MQ=255)
             ccACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggc                                                                            <  1:1013947/93‑1 (MQ=255)
             ccACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggc                                                                            >  1:168973/1‑93 (MQ=255)
               acgacgATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggcgg                                                                          >  1:3025435/1‑93 (MQ=255)
                cgacgaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggcggg                                                                         >  2:1060662/1‑92 (MQ=255)
                cgacgaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGt                                                                         >  1:2978961/1‑93 (MQ=255)
                cgacgaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGt                                                                         <  1:3435738/93‑1 (MQ=255)
                cgacgaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGt                                                                         <  2:4531568/93‑1 (MQ=255)
                  acgaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTgg                                                                       <  1:3995839/93‑1 (MQ=255)
                    gaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCt                                                                     >  2:3958646/1‑93 (MQ=255)
                     aTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTa                                                                    <  2:2640399/93‑1 (MQ=255)
                     aTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTa                                                                    >  2:4493307/1‑93 (MQ=255)
                       gCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACa                                                                  <  1:145987/93‑1 (MQ=255)
                         aGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGa                                                                <  1:3240899/93‑1 (MQ=255)
                         aGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGa                                                                <  2:1755666/93‑1 (MQ=255)
                         aGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGa                                                                <  2:4257538/93‑1 (MQ=255)
                         aGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGa                                                                <  2:2603987/93‑1 (MQ=255)
                          ggCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACaa                                                                                >  2:88854/1‑76 (MQ=255)
                          ggCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACaa                                                                                <  1:88854/76‑1 (MQ=255)
                             gcTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGa                                                            <  1:810631/92‑1 (MQ=255)
                             aCTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGa                                                            <  1:2186389/93‑1 (MQ=255)
                             aCTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGa                                                            <  2:3336392/93‑1 (MQ=255)
                               tCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAAt                                                          >  1:3659674/1‑93 (MQ=255)
                               tCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAAt                                                          <  2:2191428/93‑1 (MQ=255)
                                 cGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTg                                                        >  1:1329959/1‑93 (MQ=255)
                                  gCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGc                                                       <  2:1764961/93‑1 (MQ=255)
                                   cTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCa                                                      >  2:2404077/1‑93 (MQ=255)
                                     gTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAg                                                    >  1:4275738/1‑93 (MQ=255)
                                     gTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAg                                                    >  1:2258870/1‑93 (MQ=255)
                                     gTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAg                                                    >  1:3982410/1‑93 (MQ=255)
                                      tttCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGc                                                   <  2:273968/93‑1 (MQ=255)
                                       ttCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCt                                                  <  1:1129123/93‑1 (MQ=255)
                                        tctcCGGGCTACCGCTTTTTACACTAAGCTCCATGCACTACGCTCCTGAACCTTAAAGACAGCGGCGGTGGCTACAGATAGAATTGCAAGCtt                                                 <  1:2297344/93‑1 (MQ=255)
                                        tctcCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCtt                                                 >  2:1413723/1‑93 (MQ=255)
                                         ctcCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCCCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCttt                                                <  2:2216594/93‑1 (MQ=255)
                                          tcCGGTCTACCGCTTTTTACACTAAACTCCAGGCACTACGCTCTTTTATCATAAAGACAACGGCGTTGGCTATAGATAGAATAGCAAGCCTTc                                               >  1:4685577/1‑93 (MQ=255)
                                          tcCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTc                                               >  2:4505893/1‑93 (MQ=255)
                                              ggCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTaa                                           >  1:4449837/1‑93 (MQ=255)
                                                cTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACt                                         <  1:1627212/93‑1 (MQ=255)
                                                cTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACt                                         <  2:2395302/93‑1 (MQ=255)
                                                  aCCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCa                                       <  1:4535704/93‑1 (MQ=255)
                                                  aCCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCa                                       >  2:1929370/1‑93 (MQ=255)
                                                   ccGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCAt                                      <  2:442429/93‑1 (MQ=255)
                                                    cGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATg                                     <  1:2327837/93‑1 (MQ=255)
                                                     gCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGt                                    >  2:3667495/1‑93 (MQ=255)
                                                      cTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTc                                   >  2:1383925/1‑93 (MQ=255)
                                                      cTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTc                                   >  1:3273324/1‑93 (MQ=255)
                                                       tttttACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTcc                                  >  2:3526262/1‑93 (MQ=255)
                                                        ttttACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCg                                 >  2:2417892/1‑93 (MQ=255)
                                                         tttACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGc                                <  2:1718507/93‑1 (MQ=255)
                                                          ttACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCt                               <  2:2825713/93‑1 (MQ=255)
                                                            acacTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGt                             <  2:406068/93‑1 (MQ=255)
                                                            acacTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGt                             <  2:3041251/93‑1 (MQ=255)
                                                            acacTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGt                             <  2:2977571/93‑1 (MQ=255)
                                                              acTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTg                           >  1:2637949/1‑93 (MQ=255)
                                                              acTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTg                           <  2:3701232/93‑1 (MQ=255)
                                                               cTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGc                          <  2:2257407/93‑1 (MQ=255)
                                                               cTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTAGTAACTCATGTCCGCTGTTGc                          >  1:767962/1‑93 (MQ=255)
                                                                 aaaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGa                        <  1:2020981/93‑1 (MQ=255)
                                                                 aaaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGa                        <  2:341227/93‑1 (MQ=255)
                                                                 aaaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGa                        <  2:1282271/93‑1 (MQ=255)
                                                                  aaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGAt                       >  1:1848796/1‑93 (MQ=255)
                                                                  aaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGAt                       <  2:2864621/93‑1 (MQ=255)
                                                                  aaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGAt                       <  2:578756/93‑1 (MQ=255)
                                                                    cTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGa                     <  1:1786207/93‑1 (MQ=255)
                                                                    cTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGa                     <  2:4328077/93‑1 (MQ=255)
                                                                     tCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGAc                    <  2:833031/93‑1 (MQ=255)
                                                                      ccATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGACTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACt                   <  1:24824/93‑1 (MQ=255)
                                                                       cATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACtt                  >  1:1102918/1‑93 (MQ=255)
                                                                        aTGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTc                 >  1:4603179/1‑93 (MQ=255)
                                                                           cACTACGCCCCTGAACCTTAAAGACAACGGCGTTGTCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCgtg              <  2:2881033/93‑1 (MQ=255)
                                                                            aCTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCgtgt             >  1:949503/1‑93 (MQ=255)
                                                                             cTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCAGCTGTTGCGATGACTTCGTGtt            >  1:3784531/1‑93 (MQ=255)
                                                                              tACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTa           >  2:2527410/1‑93 (MQ=255)
                                                                              tACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTa           >  2:518100/1‑93 (MQ=255)
                                                                              tACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTa           >  1:4155092/1‑93 (MQ=255)
                                                                               aCGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTaa          <  2:1983203/93‑1 (MQ=255)
                                                                               aCGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTaa          <  1:3716268/93‑1 (MQ=255)
                                                                                  cTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCt       >  2:4685321/1‑93 (MQ=255)
                                                                                   tCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCtt      <  1:383672/93‑1 (MQ=255)
                                                                                   tCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCtt      <  1:4046956/93‑1 (MQ=255)
                                                                                    ccTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTa     <  2:1149853/93‑1 (MQ=255)
                                                                                      tGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAAc   >  2:4000803/1‑93 (MQ=255)
                                                                                       gAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAACg  >  2:2148249/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAACG  >  NC_000913/4483446‑4483625

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