Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,886,722 G→A G122S (GGC→AGC)  b3708 → tryptophanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,886,7220GA100.0% 1021.9 / NA 277G122S (GGC→AGC) b3708tryptophanase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (140/137);  total (140/137)

Only 100 of 284 total aligned reads displayed.

GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCG  >  NC_000913/3886632‑3886810
                                                                                          |                                                                                        
ggTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGc                                                                                        <  1:3750786/93‑1 (MQ=255)
 gTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGcc                                                                                       >  1:1523161/1‑93 (MQ=255)
  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                      <  1:110463/93‑1 (MQ=255)
  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                      <  2:405337/93‑1 (MQ=255)
  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                      <  1:549960/93‑1 (MQ=255)
  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCAGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGccg                                                                                      >  1:164783/1‑92 (MQ=255)
   tATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTg                                                                                     >  1:793731/1‑93 (MQ=255)
    aTCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTgg                                                                                    >  1:46540/1‑93 (MQ=255)
    aTCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTgg                                                                                    >  2:4093498/1‑93 (MQ=255)
    aTCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGACTgg                                                                                    >  1:2721505/1‑93 (MQ=255)
     tCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGa                                                                                   <  2:4441075/93‑1 (MQ=255)
      cAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGAt                                                                                  <  1:1841768/93‑1 (MQ=255)
       aaTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATc                                                                                 <  2:265443/93‑1 (MQ=255)
       aaTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATc                                                                                 <  1:2297758/93‑1 (MQ=255)
         tACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAATCCTGGATCgc                                                                               >  1:4411041/1‑93 (MQ=255)
          acacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgca                                                                              >  2:3010763/1‑93 (MQ=255)
           cacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcag                                                                             >  1:291233/1‑93 (MQ=255)
            acCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagc                                                                            >  1:718881/1‑93 (MQ=255)
            acCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagc                                                                            >  1:4110054/1‑93 (MQ=255)
            acCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagc                                                                            >  2:3838370/1‑93 (MQ=255)
             ccATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagca                                                                           >  1:3180112/1‑93 (MQ=255)
              cATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaa                                                                          >  1:944392/1‑93 (MQ=255)
              cATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaa                                                                          >  2:3410124/1‑93 (MQ=255)
               aTTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGAACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaa                                                                         >  1:2191768/1‑93 (MQ=255)
                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                        >  1:3453284/1‑93 (MQ=255)
                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                        <  1:868670/93‑1 (MQ=255)
                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                        <  1:594528/93‑1 (MQ=255)
                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                        >  1:2860854/1‑93 (MQ=255)
                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                        <  2:3581672/93‑1 (MQ=255)
                 tCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAt                                                                       <  1:1488658/93‑1 (MQ=255)
                 tCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAt                                                                       <  2:2125138/93‑1 (MQ=255)
                  ccGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtg                                                                      <  2:11430/93‑1 (MQ=255)
                   cGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtgg                                                                     <  1:632395/93‑1 (MQ=255)
                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                    >  1:388953/1‑93 (MQ=255)
                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                    <  1:754987/93‑1 (MQ=255)
                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                    <  2:585459/93‑1 (MQ=255)
                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                    >  1:904796/1‑93 (MQ=255)
                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                    <  2:2161324/93‑1 (MQ=255)
                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                    <  2:4180294/93‑1 (MQ=255)
                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                   >  2:1471194/1‑93 (MQ=255)
                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                   >  2:3723036/1‑93 (MQ=255)
                      cTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtgg                                                                  <  1:1744922/93‑1 (MQ=255)
                      cTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtgg                                                                  <  1:4217400/93‑1 (MQ=255)
                      cTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtgg                                                                  <  2:2491281/93‑1 (MQ=255)
                       tCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGc                                                                 >  1:1934874/1‑93 (MQ=255)
                       tCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGc                                                                 <  1:3318823/93‑1 (MQ=255)
                           cAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtc                                                             <  1:1043842/93‑1 (MQ=255)
                           cAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtc                                                             >  2:811562/1‑93 (MQ=255)
                            aGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtct                                                            >  1:848876/1‑93 (MQ=255)
                            aGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGGGGCGGtct                                                            >  1:4198903/1‑93 (MQ=255)
                             gggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                                           >  1:1302344/1‑93 (MQ=255)
                             gggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                                           >  1:2350660/1‑93 (MQ=255)
                             gggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGAGTtctc                                                           >  2:3892036/1‑93 (MQ=255)
                              ggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTctct                                                          <  1:3486842/93‑1 (MQ=255)
                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                                       <  1:4171658/93‑1 (MQ=255)
                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                                       <  2:2823908/93‑1 (MQ=255)
                                    ggCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAt                                                    <  2:4369053/93‑1 (MQ=255)
                                    ggCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAt                                                    <  2:4237586/93‑1 (MQ=255)
                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTGTAACTAtt                                                   >  2:1053958/1‑93 (MQ=255)
                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                                   <  1:4039695/93‑1 (MQ=255)
                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                                   <  2:3238939/93‑1 (MQ=255)
                                          gagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCttt                                              >  2:1277920/1‑93 (MQ=255)
                                          gagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCttt                                              >  2:3195049/1‑93 (MQ=255)
                                          gagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCttt                                              >  2:1395556/1‑93 (MQ=255)
                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                            <  1:1727600/93‑1 (MQ=255)
                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                            >  1:3601315/1‑93 (MQ=255)
                                              aaaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa                                          >  2:620844/1‑93 (MQ=255)
                                               aaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                                         <  1:2461547/93‑1 (MQ=255)
                                               aaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                                         >  1:1351027/1‑93 (MQ=255)
                                                aTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAcc                                        >  1:668226/1‑93 (MQ=255)
                                                 tCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                                       <  2:1883836/93‑1 (MQ=255)
                                                 tCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGAGCGCAGCAAGATGGTGGCGTTGTCTAACTATTTCTTTGATACCa                                       >  2:1509420/1‑93 (MQ=255)
                                                  cTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                                      >  1:4583246/1‑93 (MQ=255)
                                                   tataTTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACg                                     >  2:1651516/1‑93 (MQ=255)
                                                       ttCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAgg                                 >  1:3824723/1‑93 (MQ=255)
                                                       ttCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAgg                                 <  2:2396027/93‑1 (MQ=255)
                                                          cGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGcc                              >  2:332544/1‑93 (MQ=255)
                                                            gTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCAt                            >  2:1994184/1‑93 (MQ=255)
                                                            gTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCAt                            >  1:3964819/1‑93 (MQ=255)
                                                             tACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTGACTATTTCTTTGATACCACGCAGGGCCATa                           <  2:2900056/93‑1 (MQ=255)
                                                             tACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGAGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATa                           <  2:1783552/93‑1 (MQ=255)
                                                              aCTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGGGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAg                          >  1:3374001/1‑93 (MQ=255)
                                                               cTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc                         <  2:3444086/93‑1 (MQ=255)
                                                               cTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc                         <  1:2473615/93‑1 (MQ=255)
                                                               cTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGCGCCATAGc                         >  2:4335889/1‑93 (MQ=255)
                                                                 gATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGTGCCATAGCCa                       >  2:704099/1‑93 (MQ=255)
                                                                 gATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa                       >  2:3784436/1‑93 (MQ=255)
                                                                 gATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa                       <  2:4555625/93‑1 (MQ=255)
                                                                  aTTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAg                      >  1:4338371/1‑93 (MQ=255)
                                                                   ttAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGa                     >  1:2555916/1‑93 (MQ=255)
                                                                    tAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGAt                    >  1:3118941/1‑93 (MQ=255)
                                                                     aaaaaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATc                   <  2:1731705/93‑1 (MQ=255)
                                                                        aaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAAc                >  2:2281475/1‑93 (MQ=255)
                                                                         aaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACg               <  1:1179310/93‑1 (MQ=255)
                                                                                gCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTAc        >  2:4631082/1‑93 (MQ=255)
                                                                                   ggAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGTGCCATAGCCAGATCAACGGCTGTACCGt     <  1:2125264/93‑1 (MQ=255)
                                                                                   ggAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGt     >  2:1209210/1‑93 (MQ=255)
                                                                                    gAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTg    <  1:1968964/93‑1 (MQ=255)
                                                                                     aaaaaaGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGc   >  2:1329112/1‑93 (MQ=255)
                                                                                      aaaaAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCg  >  1:1234821/1‑93 (MQ=255)
                                                                                          |                                                                                        
GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCG  >  NC_000913/3886632‑3886810

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